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Résumé. — Le spectre d’absorption de l’iodure mercurique rouge est dichroïque. Aux très 
basses températures, on observe un spectre formé de raies polarisées précédant une absorption 
continue de plus courte longueur d’onde. La lumière réfléchie donne lieu au phénomène de rayons 
restant et manquant au voisinage des raies d’absorption. 

On expose, dans ce mémoire, les données expérimentales quantitatives sur l’absorption et la 
réflexion. On en déduit l'intensité d’oscillateur des raies que l’on compare à la théorie d’Elliott des 
spectres de première classe. 

La raie d’absorption ordinaire principale a un comportement exceptionnel. 


Abstract. — The absorption spectrum of red mercuric iodide is diochric. At very low 
temperatures, the spectrum consists of polarised lines preceding a continuum of shorter wave- 
lengths. Near the absorption lines, a maximum and a minimum of reflected light is observed. 

In this paper experimental quantitative results on the absorption and the reflexion spectra of 
Hgl, are presented. We deduce from these data the oscillator strength of the lines which are 
compared to those of exciton spectra of the first class in Elliott’s theory. 

The principal ordinary absorption line shows peculiar properties which are discussed. 


Introduction. — Le spectre d’absorption de 
Hgl, rouge est dichroïque [1]. À 4 °K par exemple, 
on observe, avec les cristaux perpendiculaires à 
-Vaxe, une raie ordinaire intense située à 5 320 À, 
suivie du côté des courtes longueurs d’onde par une 
raie moins intense à 5 250 À précédant un fond 
d’absorption continue (fig. 1). À la même tempé- 
rature, l’absorption extraordinaire de cristaux 
parallèles à l’axe se compose d’une raie à 5 300 À 
parfois accompagnée de raies satellites et d’un fond 
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d’absorption continu débutant par un bord vers 
5 200 À (fig. 1). Dans ces mêmes cristaux paral- 
lèles à l’axe, le spectre ordinaire montre générale- 
ment une raie intense à 5 250 À précédant une 
absorption continue de plus courte longueur d’onde 
(fig. 1). Toutefois, la raie 5 320 À apparaît très 
faiblement dans quelques échantillons. Ce compor- 
tement exceptionnel n’est pas conforme aux théo- 
ries connues actuellement. 

Nous nous proposons, dans ce mémoire, d’expo- 
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ser les données expérimentales sur l’absorption et 
la réflexion de ces cristaux. On a déterminé l’inten- 
sité d’oscillateur des raies et on a comparé les 
valeurs obtenues à celles déduites de la théorie 
d’Elliott pour les spectres de première classe. 


Mesure des coefficients d’absorption et de ré- 
flexion. — Les coefficients d'absorption et de ré- 
flexion ont été mesurés à l’aide du montage repré- 
senté sur la figure 2. 

L’échantillon était éclairé en lumière monochro- 
matique à l’aide d’un monochromateur Zeiss à 
trois prismes. La lumière transmise et la lumière 
réfléchie étaient envoyées sur des photomultipli- 
cateurs R. C. A. L'étude photométrique des absorp- 
tions ordinaire et extraordinaire a nécessité l’emploi 
d’un polaroïde qui était interposé sur le faisceau 
lumineux incident à la sortie du monochromateur. 

Sur le porte-échantillon, on avait fixé, près de 
l’échantillon, un petit miroir dont le pouvoir réflec- 
teur était 0,9. Les mesures s’effectuaient alors en 
trois temps. 

a) On a d’abord mesuré simultanément les inten- 
sités transmises ZP et réfléchie Z® par le cristal, 
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recueillies dans les photomultiplicateurs indépen- 
donts 1 et 2. Il a été nécessaire ensuite de mesurer, 
indépendamment pour chaque photomultiplica- 
teur, les intensités incidentes 7, correspondant aux 
mêmes longueurs d’onde. 

b) A cet effet, on a d’abord soulevé légèrement 
le cryostat, de façon à obtenir la valeur de l’inten- 
sité incidente 12 dans le photomultiplicateur 1. 
La lumière traversait alors un support en verre non 
recouvert de Hgl.. 

c) On a soulevé ensuite le cryostat, de sorte que 
l’image de la fente du monochromateur soit formée 
sur le miroir et se projette sur le photomultipli- 
cateur 2. On a ainsi mesuré l’intensité 2P. 

Ces mesures ont été faites pour toute la région 
spectrale étudiée. 

En admettant que l’on puisse appliquer la loi de 
Lambert, le coefficient d'absorption X est donné 
par : 

14 


RUES [os + log (1 — n)| le (1) 


e désignant l’épaisseur de l’échantillon. 
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L’épaisseur a été mesurée, soit à l’aide d’un com- 
parateur, soit par une méthode interférentielle. La 
mesure interférométrique est la plus précise, mais 
elle n’est pas la plus générale car les échantillons 
doivent avoir des épaisseurs convenables et des 
surfaces régulières. On observe, dans ce cas, par 
transmission ou par réflexion un spectre cannelé de 
lames minces. L’indice de l’iodure rouge étant 
connu pour plusieurs longueurs d’onde, il est aisé 
de déduire l’épaisseur de l’échantillon de la mesure 
de la position des cannelures. On a dû utiliser, en 
général, la première méthode. L’erreur commise 
sur l’épaisseur était alors de l’ordre de + 15 
dans nos expériences. 

Par ailleurs, l’erreur dans la détermination de la 


densité optique au voisinage du maximum d’ab- 
sorption pouvait atteindre 5 %. Finalement, l’er- 
reur commise sur X n’excède probablement”pas 
+ 15 % dans le voisinage de Kax et doit être 
moins grande dans les autres régions du spectre. 
L’intensité d’oscillateur f des raies se déduit alors 
directement de la formule de KRAwETz [2]. 


es Cte f (y dv. (2) 


On n’a pas tenu compte de la variation de 
indice de réfraction. L'erreur commise est petite à 
cause de la forme antisymétrique de la courbe de 
dispersion de part et d’autre de la raie d'absorption. 
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Résultats expérimentaux, — Nous allons décrire 
les résultats expérimentaux obtenus à 290 0K 
77 0K et 4 0K. 

10 A 290 °K : On n’observe pas de raie (fig. 3). 
Les coefficients de réflexion ordinaire et extra- 
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ordinaire, pratiquement constants entre 16 000 et 
19 000 cm1, sont respectivement de l’ordre de 
R, = 0,23 et R; = 0,14 (H 0,02). 

20 À 77 0K : L’épaisseur des échantillons était 
telle qu’il n’a pas été possible d'observer la raie 
ordinaire. Les échantillons d’épaisseur inférieure 
au micron sont souvent troués, de sorte que, par 
endroits, la lumière traversait l’échantillon sans 
être absorbée. Dans ces conditions, seule la raie 
extraordinaire a pu être observée (fig. 4). 
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Le facteur f, par cellule élémentaire, de cette 
rale est de l’ordre de : 


19 12.105 + 4105. 
Sa demi-largeur a pour valeur : 
AVS, 185 cm—1. 


Le coefficient d’absorption maximum, déduction 
faite du fond d’absorption continue, est donné par : 


Kx (raie seule) < 0,8.108 cm—1 


R® reste constant entre 48 000 et 20 000 em—1 
et a pour valeur : 
RQ 0,14 + 0,02. 


Au contraire, le pouvoir réflecteur des échantil- 
lons perpendiculaires à l’axe présente des ano- 
malies au voisinage de la raie d’absorption. On 
observe un maximum de réflexion à 18 670 cm1 
et un minimum à 18 800 cm1. 

Ces accidents de réflexion sont dus à la disper- 
sion anormale au voisinage de la raie d’absorption 
et ont été respectivement appelés «bandes de 
rayons restants » et «bandes de rayons man- 
quants » par S. Nikitine et R. Reiss [3]. Maïs, ces 
phénomènes sont plus nets à 40K. 


39 A 40X. 


a) CRISTAUX PERPENDICULAIRES A L’AXE. — 
On observe la raie ordinaire principale (O0. P.) à 
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5 320 À. Le facteur f, par cellule élémentaire, de 
cette raie est de l’ordre de (fig. 5) : 


OP 192 104 + 4.104 
Sa demi-largeur a pour valeur : 
AVS: 160 cm—1 
et le coefficient d’absorption maximum 
aie (raie seule) æ 6.10% cmt. 


Par réflexion, on observe au voisinage de la raie 
O0. P. (fig. 6) des bandes de rayons restant et man- 
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coefficient de reflexion À 
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quant. Comme l’indice d'absorption reste toujours 
petit (*), on peut utiliser la formule 


Ref) (3) 


qui permet de tracer la courbe de dispersion au 
voisinage de la raie. En appliquant les formules de 
Kramers : 


; L Ne? PRESS &? 
n2— k2 = n2, + 27 : Coë = DE 5 
m  (oÿ— &°?)? + yo (a) 
Æ. 
ne kr Ne? vof 
; m  (oÿ— w?)}? + y? w? i 


(où l’on désigne par y le coefficient d’amortisse- 
ment, k l’indice d'absorption, N le nombre de cel- 
lules élémentaires par cm’, w, la pulsation pour la 
raie d'absorption, m la masse de l’électron, c la 
vitesse de la lumière et f l'intensité d’oscillateur), 
on retrouve approximativement la courbe de dis- 
persion en posant : 


20 LD 
RÉ 
Dans ce calcul, on a supposé nr légèrement 
variable (fig. 7). Cette valeur de f n’est pas en bon 
(*) Au voisinage du maximum de la raie d'absorption, 
l’indice d'absorption k est de l’ordre de 0,025. Cette valeur 
est nettement plus petite que (n —1)?. Au voisinage du 


rayon manquant pour lequel ce terme devient très petit, 
k est également beaucoup plus petit que 0,025. 
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accord avec la valeur mesurée par absorption. 
Mais elle est probablement fortement entachée 
d'erreurs, car il n’est pas certain que l’on puisse 
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appliquer sans restrictions la formule (3). En outre, 


il n’est pas possible d’éliminer les contributions à 
l'indice dues au fond continu et à la raie 5 250 À. 


b) CRISTAUX PARALLÈLES A L’AXE. — Spectre 
extraordinaire : Quand le vecteur électrique est 
parallèle à l’axe, on observe une raie d’absorption 
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extraordinaire principale (E. P.) à 5 296 À dont le 
facteur f a pour valeur (fig. 8) : 


FT 10.105 + 3.105, 
Sa largeur est de : 
AVÈR- 85 cm1 
et son coeflicient d'absorption maximum de : 
Knax (raie seule) æ 8.102 cm—1. 


Le coefficient de réflexion ne présente pas d’ano- 
malies au voisinage de la raie et conserve une 
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valeur pratiquement .constante entre 17 500 et 
19 500 FE : 


R 0,14 + 0,02. 


Spectre ordinaire : Quand le vecteur électrique 
est perpendiculaire à l’axe, on devrait retrouver le 
spectre ordinaire décrit ci-dessus. En fait, on ob- 
serbe une raie très intense (0!) vers 5250 À pré- 
sentant les mêmes caractéristiques que la raie 
O. P. à 5 320 À (fig. 9): 
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Nous n'avons pas trouvé d’explication satis- 
faisante de ce phénomène. Nous nous proposons de 
poursuivre des expériences à ce sujet. La raie 
5 320 À qui a été observée avec quelques échan- 
tillons s’est trouvée être trop faible pourles mesures 
du facteur f avec les échantillons utilisés à cet effet. 

Au contraire, avec le même cristal, on observe 
par réflexion les mêmes phénomènes que pour les 
cristaux perpendiculaires à l’axe. Ceci semblerait 
indiquer que la raie O. P. existe bien dans les 
couches superficielles du cristal. Comme la raie 
5 320 À a été observée dans certains cristaux fins, 
il n’est done pas exclu que la préparation du cristal 
ne joue un certain rôle dans les phénomènes d’ab- 
sorption. 

Une autre hypothèse plausible consisterait à 
admettre l’existence de tensions provoquées par le 
support de verre, qui déplaceraient le spectre ordi- 
naire des cristaux parallèles à l’axe vers les courtes 
longueurs d'onde. Toutefois, le spectre extraordi- 
naire. ne semble pas être affecté par ces tensions 
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comme l’étude des cristaux plus épais, parallèles à 
l'axe, indépendants de tout support, l’a montré. 
Enfin, une dernière possibilité serait celle d’ad- 
mettre que les tensions existant près de l’interface 
verre-cristal modifient l’intensité des raies. 


Classification des spectres. — Le spectre ordi- 
naire présente les caractéristiques d’un spectre 
excitonique de première classe de la théorie d’El- 
liott [4]. Il se compose, en effet, d’une bande princi- 
pale très intense (O. P.) suivie, du côté des courtes 
longueurs d’onde, d’une bande moins intense (0. 
S.). D’après cette théorie, on doit pouvoir classer 
les raies dans une série hydrogénoïde : | 


_ R° 
Vas NES e RM 2: se. CO (5) 
avec 
, _ 
ÊE emo (6) 


où À désigne la constante de Rydberg, u la masse 
réduite de l’exciton, e la constante diélectrique, 
m, la masse de l’électron. 

Les deux raies d’absorption déterminent les para- 


mètres À’ et vs, soit : 
Va = 19 119— (330/n°) 


NAN Oo 


La limite de la série, donnée par ve — 19 119. 
coïncide alors pratiquement avec le bord d’absorp- 
tion. On en déduit : 


Lord © 0,2 M 
et si l’on suppose que Ma © Mtrou : 
mel 0m 
qui est à rapprocher de mur = 0,54 m, [5]. Le 
rayon de la première orbite serait : 
af © 98 À. 


Par ailleurs, l’application de la formule proposée 
par S. Nikitine et H. Haken [6] donnant l’intensité 
d’oscillateur par cellule élémentaire pour les spec- 
tres de première classe : 


= = 7 
fr = fs (7) 


où V, désigne le volume de la cellule élémentaire 
(V, (Hgl,) — 238,46 À5) et fr un facteur calcu- 
lable approximativement 


(5 = 5 


27? MVno 
LE 


*% 
fr grad w, dr, 


où uw, est la fonction de Bloch décrivant les états 
les plus bas de la bande de conduction (4=— 0) et u, 
la fonction de Bloch*décrivant les états les plus 
élevés de la bande de valence (4 — 0)} donne 
pour HeT, : 


fn # 0,35, 
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Cette valeur est plus petite que celles qui ont 
été déterminées pour Cul et CuCl [8] : 


fe (Cul) 1 
fe (CuCl) = 0,9. 


Le spectre extraordinaire pourrait correspondre 
à des transitions défendues. 


Tentatives d'interprétation. — Suivant une sug- 
gestion de S. Nikitine, H. Haken a remarqué que 
le dichroïsme des spectres excitoniques provient 
essentiellement du terme fs (formule 7). Par ail- 
leurs, on peut admettre, dans une première approxi- 
mation (dite de liaison forte), que les fonctions 
d’onde qui interviennent dans f# sont des fonc- 
tions atomiques. Différents cas sont alors possibles. 

En prenant, par exemple, des fonctions ato- 
miques S pour la bande de valence et des fonctions 
atomiques P pour la bande de conduction, fB est 
différent de zéro et la transition est permise. Au 
contraire, si la bande de conduction est composée 
de fonctions atomiques D, la transition est inter- 
dite. 

Ces résultats s’appliquent également aux transi- 
tions excitoniques. 

Il est alors possible d'interpréter qualitative- 
ment les spectres d’absorption de Hgl, si l’on 
admet que l’état fondamental est du type S et la 
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bande de conduction du type P. A cause de l’aniso- 
tropie, les orbitales P, et P;, auront des énergies 
différentes, l’axe du cristal étant dirigé suivant Oz. 
On voit, en outre, que la transition vers les orbi- 
tales P, et P;, ne sera permise que pour les pola- 
risations parallèle et perpendiculaire à l’axe, res- 
pectivement. La transition $ — P, correspondrait 
alors au spectre extraordinaire et la transition 
S — P;, correspondrait au spectre ordinaire. Cette 
discussion s’applique également aux transitions 
excitoniques. 

Il est possible d’évaluer expérimentalement la 
différence d’énergie entre les bandes de conduction 
ordinaire et extraordinaire : 


E(P;— P;,) 565 cm? 


Toutefois, l’anisotropie très forte du cristal pro- 
voque probablement une très forte anisotropie des 
masses effectives et il est probable que dans un tel 
réseau cristallin, le modèle hydrogénoïde et les for- 
mules d’Elliott ne s’appliquent plus. De ce fait, 
la comparaison des intensités d’oscillateur des raies 
ordinaire et extraordinaire principales devient très 
problématique. 

Les auteurs remercient le D' Elliott pour les dis- 
cussions fructueuses qu'ils ont eues avec lui. 


Manuscrit reçu le 19 juin 1961. 
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LA DILATATION THERMIQUE DU MILIEU CRISTALLIN 


Par J. LAVAI, 
Collège de France, Paris. 


IV 


Les coefficients d’élasticité, les tensions et la dilatation thermique. 


Résumé. — Toute déformation isotherme qui entraîne une variation de volume fait augmenter 
ou diminuer les tensions thermiques. En conséquence, tout coefficient d’élasticité qui se rapporte 
à une variation de volume renferme une composante égale à la tension thermique. 

À température croissante, quand le milieu cristallin se dilate libre de toute contrainte, les coef- 
ficients d’élasticité diminuent (en grandeur absolue), et cette diminution est déterminée par la dila- 
tation thermique. La diminution calculée paraît en bon accord avec les données expérimentales. 

Enfin, une dernière étude est consacrée à la variation que le coefficient de compressibilité subit 
lorsque le milieu cristallin, maintenu à température constante, est soumis à une pression croissante. 


Abstract. — Any isothermal deformation leading to a variation of volume causes thermal 


strains to increase or decrease. 


Consequently, any elastic coefficient related to a variation of volume 
includes a component equal to the thermal strain. 


At increasing temperature, when the crystalline medium expands, free from all constraints, the 
elastic coefficients decrease (in absolute magnitude) and this decrease is determined by the thermal 


expansion. 


The calculated decrease gives good agreement with the experimental data. 


Finally, a study is made of the variation of the compressibility coefficient when the crystalline 
medium, kept at constant temperature, is subjected to an increasing pressure. 


A. Les tensions thermiques 
et les coefficients d’élasticité. 


Cas des cristaux cubiques formés par les mé- 
taux (1). — Quand le milieu cristallin subit une 
contraction ou une dilatation isotrope à tempéra- 


‘ ture constante, chaque facteur gs. (9) et chaque 


pulsation ©, (67) varient. Maïs tant que la con- 
traction et la dilatation restent minimes, tant 
qu’elles développent des tensions conformes à la loi 
de Hooke, tous les rapports (24) 


Esr = gs: [@°s+ 


restent à peu près constants ; car la croissance ou la 
décroissance relative des produits 


Da, D Ya 


est, dans l’ensemble, à peu près égale à celle des 
dérivées secondes CS, et partant à celle des carrés 
des pulsations 24, (21). D'autre part, l’élévation ou 
l’abaissement relatif des pulsations w4$, étant très 
petit, chaque énergie Ws. (3) n’éprouve pas de 
variation notable, même aux basses températures. 


, 


(1) Les formules numérotées de 1 à 31, les rétérences de 1 
à 5 et les figures 1 et 2 se trouvent dans la première partie ; 
les formules numérotées de 32 à 53, lessréférences de 6 à 10, 
et les figures 3, 4, 5, 6, 7 dans la seconde partie ; les for- 


.mules numérotées de 54 à 69, la référence numérotée 11, la 


figure 8 dans la troisième partie. 


Pour des raisons semblables, une faible contrac- 
tion ou dilatation du milieu cristallin, à tempéra- 
ture constante, ne fait pas varier considérablement, 
non plus, l’énergie moyenne W,, (27-28-29) ni les 
coefficients b (28) et c (29). 

Cela suppose la dilatation ou la contraction iso- 
trope ; en général 


Avfv = À tes. B—1,2,3 (70) 


et les trois dilatations (gg > 0) ou contractions 
(te < 0) élémentaires diffèrent entre elles. Mais, 
même si la dilatation ou la contraction se réduit à 
une seule composante, teg (B constant), en grandeur 
absolue, les variations qu’elle imprime aux dérivées 
secondes et troisièmes de l’énergie potentielle sont 
de même sens ; dans l’ensemble, les modules de ces 
dérivées décroissent s’il y a dilatation et croissent 
s’il y a contraction ; donc les facteurs g4, et les 
carrés &?4 des pulsations augmentent à la fois ou 
diminuent à la fois ; et tant que la contraction ou la 
dilatation élémentaire {gg reste minime, tant qu’elle 
développe des tensions proportionnelles à son am- 
plitude, tel, les rapports Es. = gs. /&?$., les éner- 
gies Ws et, de même, les coefficients b et c et 
l'énergie W,, ne subissent pas de variation notable. 


(1) tes sont les déformations élémentaires : si u est le 
déplacement du point x 
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Une déformation de cisaillement #3 (5 Æ £) 
accroît les distances entre certains atomes et les 
diminuent entre d’autres. En conséquence, elle 
augmente des dérivées secondes C& (1-2) tandis 
qu’elle réduit d’autres dérivées du même ordre. 
Mais les distances accrues sont en même nombre 
que les distances diminuées, et les dérivées secondes 
affaiblies en même nombre que les dérivées secondes 
renforcées. Et si la déformation #5 est uniforme, 
c’est-à-dire de même grandeur [ts] partout — ce 
que je suppose — et si elle reste en outre minime, 
développant des tensions en raison directe des dé- 
placements, dans tous les éléments de la matrice 
de Fourier (20), les croissances et les décroissances 
des dérivées secondes se compensent sensiblement ; 
de la sorte, les pulsations «4, des oscillations fon- 
damentales ne sont pas sensiblement modifiées ; 
les augmentations et les diminutions des dérivées 
troisièmes se compensent de même dans les expres- 
sions des facteurs g$, (9) qui ne subissent donc pas, 
non plus, de variations notables ; et pour la même 
raison, les coefficients b (28) et c (29) ne sont pas 
sensiblement altérés. 

Ainsi, tant que le milieu cristallin reste le siège 
d’une déformation isotherme qui développe des 
forces de rappel obéissant à la loi de Hooke, la 
tension thermique a sensiblement pour expression : 


To=rP ID. (71) 
B étant une constante : 
BL (Q + 0% + ea) 
donc 
RTeldv)e = — Bu? = — Tofv (71 bis) 


et compte tenu de (70) 


ATe = — T® > te. 


1. COEFFICIENTS D’ÉLASTICITÉ ISOTHERME. — 
Dans une étude précédente [5] (1), j’ai rapporté les 
coefficients d’élasticité isotherme, Merkë (ou cg F8) 
aux translations du milieu cristallin et aux dérivées 
secondes de son énergie potentielle. Je suppose le 
milieu libre de toute tension interne. Cette condi- 
tion n’est pas remplie. À toute température, même 
au zéro absolu, le milieu cristallin est le siège de 
tensions thermiques (2). 


(2) Dans cette étude je désigne les coefficients d’élasticité 
isotherme par N4Y-P5, c’est-à-dire N?4Y-P8 — cay.88, 

(?} De nombreux auteurs font état d’une pression interne 
exercée par les électrons de conductibilité dans les cristaux 
métalliques. Ils définissent cette pression conformément 
aux lois qui régissent la statistique de Fermi et de Dirac. 
C’est tenir les électrons de conductibilité pour des électrons 
libres, et c’est envisager la surface du cristal comme la 
paroi rigide d’un récipient fermé. Les électrons de conduc- 
tibilité sont dans le champ de forces des ions positifs. Et 
c’est l’énergie potentielle mutuelle des électrons de conduc- 
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Chez les cristaux cubiques formés par les métaux, 
la dilatation thermique est isotrope. Elle développe 
des tensions 74% (x — 1, 2, 3) qui ont la même 
intensité 3 a 14 = Ta = 1 Et la condition 
pour que le milieu cristallin, libre de toute con- 


frainte, soit en équilibre se traduit par l’égalité [2]. ! 


Ta + T9 = 0: (72) 

Imprimons au milieu cristallin, abandonné à lui- 
même, dont la dilatation thermique s’élève à eg, 
une légère déformation isotherme  (fgs) (°). 
(|tgs| € 1). Les dilatations sont portées à eg + tag. 
Et la tension thermique et la tension développée 
par la dilatation deviennent respectivement 
To + AT® et Ta + AT. Au lieu de la tension 
Ta34 égale et opposée à la tension Tel1z4 (fig. 2), 
prenons en compte la tension T,11,4, exercée par le 
demi-cristal À sur le demi-cristal B, donc de même 


—— 
sens Ox, que la tension Te1144. La condition d’équi- 
libre devient : 


Ta 48 = To" 84 (72 bis) 


et 
ATq1 48 = = co2E8 teg 


ATo!184 = — To!184 = 1pg. 


Ainsi, compte tenu de la condition d’équilibre 
(72 bis) la déformation isotherme (#4) développe 


suivant Ori, une tension 
6 = - (co!:88 + To) tep 
et, en général, des tensions 
Ga — > (Cot*88 + To) te; à, = 1,2, 8. 


La tension 6H a pour expression explicite : 


6 = {ce} + Te) tix + (Co !#?? + Te) (te, + ts) 


(co 1133 = co112?), 


tibilité et des ions positifs qui assure en majeure partie la 
cohésion des atomes qui forment les cristaux métalliques 
(Wigner, Seitz, Fuchs, n21, [13], [14]}). Ainsi, les électrons 
de conductibilité sont des électrons liés ; ils n’exercent point 
de tension interne. D’autre part, les oscillations de plasma 
telles qu’elles sont présentement définies [15], ne déve- 
loppent pas, non plus, de tensions internes, précisément des 
tensions qui soient notables. Donc, à des différences infimes 
près, les tensions internes qui s’exercent dans le milieu 
cristallin sont dues exclusivement à l’agitation thermique 
des atomes. 

Aux basses températures, la dilatation n’obéit plus exac- 
tement à la loi de Grüneisen. Klémens interprète cet écart 
par une dilatation d’origine électronique [16]. Mais, comme 
je l’ai montré plus haut, la loi de Grüneïsen n’est plus 
valable aux basses températures. 

(?) Par (tes) je désigne le tenseur de la déformation, 
soit les 9 déformations élémentaires tes ; 6, 8 = 1, 2, 3. 


| 
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Done, les tensions G% se rapportent aux coef- 
ficients d’élasticité (isotherme) apparents : 


; , 
ca #1 — CO het _ To ; ca? 2 C@ :122 + Te. (73) 


Les déformations de cisaillement t95 (3 7 ) 
n’entraînent nulle variation de volume, partant 
nulle variation sensible de la tension thermique ; 
donc le coefficient cel??? reste inchangé (1) : 


L 
cg "1? = ç12,12, 


Et le coefficient de compressibilité (isotherme) 
apparent a pour expression (2) : 


’11:12 ?1122 
Te — CES To TE (74) 
Cas général. — Ua cristal triclinique est le siège 
de six tensions thermiques distinctes 7'@%* 
(Te*= TT), et ses déformations ont en général 
pour effet de déplacer les uns par rapport aux 
autres les différents réseaux formés par les atomes 
en même position dans la maille [18], [5]. Mais entre 
les limites où la loi de Hooke reste valable, ces dé- 
placements relatifs restent extrêmement petits. [Ils 
r’entraînent pas une variation considérable des dé- 
rivées secondes C%, troisièmes D$Y, etc. (1-2) de 
l’énergie potentielle ; partant ils n’apportent pas 
des modifications importantes aux pulsations ws, 
des oscillations fondamentales, ni aux fac- 
teurs gs, (9) ni aux coefficients b (28) et c (29). Et, 
tant que ces modifications restent négligeables, une 
légère déformation isotherme (#%$) développe dans 
un milieu cristallin triclinique des tensions : 


Bar = À coT-B8 tes + To®t leg. 
FE 


D'où l’on conduit les coefficients d’élasticité (iso- 
therme) apparents : 


cat — co" + To°Y 
mais, si à diffère de £, 


Corse PES (À £ P). 

2. COEFFICIENTS D’ÉLASTICITÉ ADIABATIQUE. — 
Dans un milieu cristallin, libre de toute contrainte, 
à température constante et uniforme, la masse spé- 
cifique est en tous points la même à petite échelle : 
c’est-à-dire, mesurée dans une sphère dont le rayon 
minime, égal à quelques dizaines d’angstrôms, est 
plus petit que le rayon d’action sensible r des forces 
répulsives qui’forment les tensions thermiques [2]. 


(1) Le coefficient c1212 = 1212, de même que le coef- 
ficient c1212 (Voigt) n’est pas affecté par la tension ther- 
mique [5]. | à DURE 

(2) Dans le calcul du coefficient de dilatation, je fais état 
du coefficient de compressibilité T' (8) et non du coefli- 
cient I” (74) ; mais, dans le calcul de l’accroissement qui 
est apporté à la tension thermique par une élévation de la 
température, je tiens compte de l’affaiblissement subi par 
cette tension du fait de la dilatation (33). 
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Parmi les composantes harmoniques de l’agitation 
thermique, les oscillations du genre acoustique, qui 
ont les plus basses fréquences, si elles ne sont pas 
strictement transversales, développent dans le 
milieu cristallin des contractions et des dilatations 
locales qui s’étendent sur des distances plus grandes 
que le rayon d’action sensible r des forces répul- 
sives ; mais ces dilatations et ces contractions sont 
compensées à chaque instant par les variations de 
volume qui sont produites par les autres compo- 
santes harmoniques de l’agitation thermique. De la 
sorte, dans le milieu cristallin abandonné à lui- 
même et à température constante et uniforme, les 


tensions thermiques ont partout la même intensité. 


D'autre part, du zéro absolu aux températures 
moyennes, tant que les oscillations secondaires 
restent négligeables, les oscillations fondamentales 
de l’agitation thermique restent conformes aux 
oscillations déterminées par la matrice de Fou- 
rier [20]. Cela étant, les oscillations fondamentales 
de basses fréquences se rapportent aux coefficients 
d’élasticité dynamique DY [5], déduits direc- 
tement de la matrice de Fourier, ou aux coefficients 
d’élasticité isotherme de la théorie classique cgrfs, 
définis par Voigt et par Max Born [17], [18] (), si 
l’on admet que les forces de cohésion qui unissent 
les atomes engagés dans le milieu cristallin peuvent 
être tenues pour centrales. Sous l’action des ondes 
élastiques, le milieu cristallin diffuse les rayons X, 
les électrons, les neutrons. Et de l’intensité que 
prennent ces diffusions, il est possible de déduire les 
coefficients d’élasticité dynamique D'$ où les coef- 
ÉGents d’élasticité isotherme de Voigt-Max Born 
cg". 

Mais cette uniformité des tensions thermiques 
dans un cristal abandonné à lui-même est rompue 
par une oscillation étrangère, excitée dans le milieu 
cristallin par des impulsions externes ou par l’effet 
piézoélectrique. Comme toute composante harmo- 
nique de l’agitation thermique et suivant la même 
loi, l’oscillation excitée développe des tensions in- 
ternes [2]. Mais si grande que soit son amplitude, 
les tensions qu’elle produit restent minimes, négli- 
geables par rapport aux tensions thermiques, même 
au zéro absolu de la température. C’est par une 
autre action que l’oscillation excitée influe sur les 
tensions thermiques. Si elle n’est pas strictement 
transversale, elle imprime au milieu cristallin des 
contractions et des dilatations locales qui se répè- 
tent périodiquement dans l’espace et le temps et 
qui s'étendent sur des longueurs très grandes par 
rapport au rayon d’action sensible, r, des forces 
répulsives. Ainsi l’oscillation excitée fait osciller à 
grande échelle le volume spécifique du milieu cris- 


(2) Les coefficients c@"P dont je fais état sont ceux que 
Max Born a définis dans le traité intitulé « Atomtheorie des 
Festenzustandes » Max Born suppose les forces de cohésion 
centrales et rapporte de la sorte les coefficients d’élasticité 
de Voigt à la structure atomique du milieu cristallin, 
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tallin ; elle fait donc osciller, de même, les tensions 
thermiques. Ces oscillations des tensions ther- 
miques étant rapides, augmentent ou diminuent 
localement l’énergie thermique, donc élèvent ou 
abaissent localement la température selon qu’elles 
sont produites par une contraction ou par une dila- 
tation du milieu cristallin. En conséquence, l’oscil- 
lation excitée se rapporte aux coefficients d’élas- 
ticité adiabatique apparents (1) : 


£ à 
Des = Dark 2 To: (75) 


si l’on adopte la théorie que jai proposée [5], ou 
aux coefficients 


co vPE = cev8e “e Ta t (75 bis) 


si l’on adopte la théorie classique (Voigt-Max Born) 
[17], [18]. Bien‘entendu si 3 diffère de B, on a : 


1Y8 Ô 
Doùs Re Doi 
cerrre 42 


ou ne ee 

On retombe sur les relations (75-75 bis) par le 
calcul qui donne les équations classiques, relatives 
à la propagation des oscillations acoustiques [20]. 
Si le cristal est sillonné par un train d’ondes élas- 
tiques planes, chacun de ses éléments est soumis à 
une force (oscillante) df ayant pour composantes : 


Ad JET nas. 
di > É (al dx1 dx? dx 


CV 1270. 


Et selon la théorie classique de l’élasticité [20] la 
densité d’énergie Æ se rapporte aux coefficients 
d’élasticité adiabatique apparents «vf (75 bis). 

Ces conclusions concernent une oscillation excitée 
dont la longueur d’onde À est beaucoup plus grande 
que le rayon d’action sensible, r, des forces répul- 
sives qui forment les tensions thermiques [2] (soit à 
une oscillation ayant une fréquence moindre 
que 1010 s1). 

Envisageons une oscillation excitée dont la lon- 
gueur d’onde À diminue continuellement. Quand 
cette longueur d’onde devient inférieure au rayon 
d’action sensible, r, des forces répulsives, les varia- 
tions locales que l’oscillation excitée imprime aux 
tensions thermiques s’atténuent progressivement, 
puis finissent par s’annuler, et les coefficients d’élas- 
ticité auxquels se rapportent l’oscillation excitée, 
égaux primitivement aux coefficients d’élasticité 


(2) Rappelons que les coefficients d’élasticité adiaba- 


tique CTP sont liés aux coefficients d’élasticité isotherme 
c&"P8 par la relation 
xuY APS 
cotY»PÈ = co2Y>P8 +- = = O 
PC 


où X2 A6 sontles coefficients thermiques des tensions To%* 
et T'olÿ, © Ja chaleur spécifique à volume constant, et pla 
masse spécifique [19]. 
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adiabatique apparents, ce%88 (75 bis), tendent vers 
les coefficients d’élasticité isotherme c@°*F8 (72). 

Chez les cristaux cubiques formés parles métaux, 
à la température ordinaire, la tension thermique Te 
n’atteint que le centième, environ, du coefficient de 
compressibilité isotherme T (8). Mais aux tempé- 
ratures élevées, proches du point de fusion. elle 
peut excéder le dixième du même coefficient. 


B. Variation des coefficients d’élasticité 
avec la température. 
Cas des cristaux eubiques formés par les métaux. 


À température croissante, lorsque le milieu cris- 
tallin se dilate libre de toute contrainte externe, la 
dérivée (négative) (69) 


d Log.ci1[d Log.v (?) | 


reste en gros constante ; et, à la température ©, le 
coefficient d’élascitié apparent cg! (73) a pour 
expression approximative : 


’ 
1 
Cg1 — 


(50) = Fe Be) ao) 


où cl est, au zéro absolu, la grandeur du coef- 
ficient ci, et oll un coefficient constant, égal et 
de signe opposé, au module moyen de la dérivée 
d Log.cl/d Log.v, entre le zéro absolu et la tem- 
pérature de fusion. A la même approximation nous 
avons : 


11 15(® 
Co op. ? F2 + To 


deg* HO = — B(@) 11 c'e *° + dTefd0. (76) 


Quand la température tend vers le zéro absolu, le 
coefficient de dilatation cubique B(O) et la dérivée 
dTe/4@ s’annulent, donc la dérivée dc'11/4d@ tombe 
aussi à Zéro. 

Aux températures moyennes, tant que la déri- 
vée dx/d® ((54) et fig. 6)) reste sensiblement cons- 
tante, l’exposant 11 $(O) reste petit, de quelques 
centièmes au dixième seulement, et 


11 
RE Ce 


âT® 


(1— e23(0) +7 on 310) + T6. 


La température s’élevant, le produit 
Fee 
BE) (1— ant (0) +5 qu 540) 
conserve grossièrement la même grandeur, car 


() ce, cl?, c sont les symboles à deux indices des 
coefficients d’élasticité (Voigt) 


cu = cl: çlè = çll22: Ad 623,28 — p12,12, 
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l’accroissement du coefficient B(@) compense à peu 
près la décroissance du facteur 


(1— en 0) +5 qu (0) 


et comme la dérivée (positive) 47% /d0 est encore 
minime par rapport au produit 


89) et at [1 — g100) +3 en 0°), 


la dérivée dc’#t /dO reste à peu près constante. 
Aux températures élevées, proches du point de 


. fusion, la dérivée d7 4/40 tend vers 


* 


(kJuv) [3E(1 — Be). 6) 
+ bKO(2 — Po). O) + ck? O2(3 — Bio).0)] 


Elle atteint, de la sorte, le dixième environ de la 
dérivée 


del [4® = — B(O) pt œite 42, (77) 


Ainsi, quand la température s’élève, dès le zéro 
absolu, le coefficient d’élasticité isotherme appa- 
rent ct diminue d’abord lentement (jusqu’à 
O — 0} (fig. 6)) (branche AB, fig. 9), puis sa dé- 
croissance devient à peu près linéaire (branche BC, 
fig. 9), enfin au voisinage du point de fusion elle 
est moindre que celle du coefficient ckt (77) 
{branche CD, fig. 9); car aux températures élevées 
la tension thermique, T@, n’est plus négligeable par 
rapport aux coefficients d’élasticité isotherme c&. 

Pour les mêmes raisons, nous avons, à la tempé- 
rature ©, de même approximativement 


: 12 —w12 5 
lonercie 0 ET 


Te ner (5 


4 44 —4445(® 
CE = e Hô(0) 


(78) 


oùc,2?,c*#etl',sont,au zéro absolu,respectitement 


les grandeurs des coefficients d’élasticité, c&, c$ 
et du coefficient de compressibilité l'& isotherme. 
Les coefficients ol, o1?, o44et ©, diffèrent entre 


_ eux car, lorsque la température s’élève, les modules 


des dérivées secondes C% (1-2) ne diminuent pas 
uniformément. Remarquons que le coefficient 
d’élasticité apparent c$**est indépendant de la ten- 
sion thermique (78). En conséquence, si la même 
élévation de la température réduisait dans la même 
proportion les modules de toutes les dérivées se- 


. condes, à température croissante, la diminution du 


coefficient d’élasticité apparent cg serait plus 
rapide que celle des coefficients cg1 et c&?, et par- 
tant que celle du coefficient de compressibilité 


To (74) | 


Les expériences qui ont pour objet la variation 


(1) Développant le raisonnement dont j’ai rappelé les 
grands traits, Grüneisen définit la variation du module de 
compressibilité avec la température [10]. J’obtiens un 


. résultat tout différent de celui que trouve Grüneisen. 
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des coefficients d’élasticité avec la température [21], 
[22] donnent des résultats qui paraissent confirmer 
les conclusions précédentes. En particulier, les 
recherches faites par W. C. Overton et J. Gaffney 
sur l’élasticité des cristaux de cuivre [22] révèlent 


AL 


FIG19; 


que les dérivées dc&!/d0, dcé?/dO, dc'#/dO, 
s’annulent bien au zéro absolu de la température. 
Mais les conclusions précédentes supposent que 
la tension thermique est isotrope et partant que la 
dilatation thermique est de même. Elles concernent 
donc seulement l’élasticité des cristaux cubiques. 
Un cristal triclinique est le siège de six tensions 
thermiques T@°* distinctes et, quand latempérature 
s’élève, la diminution que subissent ses coefficients 
d’élasticité est fonction de toutes les composantes 
nes de la dilatation thermique (6, à — 1, 2, 3). 
Cependant, il est à présumer que cette diminution 
est déterminée principalement par l’augmentation 


relative de volume, Ÿ ge, produite par la dila- 
- 8 
tation. 


C. Variation de la compressibilité isotherme 
avec la pression. 
Cas des cristaux eubiques formés par les métaux. 


A température constante et tant que la pression 
exercée sur le milieu cristallin développe des forces 
de rappel sensiblement conformes à la loi de Hooke, 
la tension thermique Te (71) — nous l’avons vu — 
reste sensiblement en raison inverse du volume spé- 
cifique V (58) (ou du volume de la maille élémen- 
taire v). Admettons que cette propriété subsiste 
encore sous les fortes pressions, lorsque la limite de 
l’élasticité est dépassée. D’après l’égalité (71 bis) 
nous avons : : 


(à Log.T'o [à Log.v)o = —1 (79) 


788 


et le rapport 7° /l restant toujours petit par rap- 
port à l’unité, l’expression (74) donne : 


É Te) He) É TE) 
d Log.v/6e  \d Log.v/6 


d Log.T'\ 1 ](Z eo. Té Te. 
 - Log. Se to T2 À j° 


La dérivée (négative) 


F (80) 


(ù Log.1'/d Log.v)e 
est moindre, en grandeur absolue, que la dérivée 
d Log.l'/d Log.v 


prise par rapport à la dilatation thermique du 
cristal libre de toute contrainte extérieure : 


d Log.T É FE) il () 
d Log.v ùd Log.v/e  B(O)T \d8/, 


Et cette différence (négative), 


mon (re). 


croît, en grandeur absolue, quand la température 
s'élève. 


(81) 


Posons : 
d Log.T'/d Log.v = — 6; (d Log.1'/ Log.v)e = — 
(à Log.T”/ Log.v)e = — &. (82) 
Compte tenu de (81), nous avons : 
gp << 


et compte tenu de (80) : 


F ; F Té 
(reg (ets) 


Donc la différence @ -— € croît avec la tempé- 
rature. 
Supposons la dérivée 


(d Log.I"”/à Log.v)9 = — té 


invariable et la température constante. Si p est la 
. pression exercée sur le cristal, 


dIMIALE=NC 
D'où nous déduisons : 
1] = Te + Cp 
2) AVIV = — dpJ(Ts + Cp), 


l étant la grandeur vers laquelle tend le coefficient 
de compressibilité isotherme apparent I” (74) 
lorsque la pression tombe à zéro. 

Soicnt : V, le volume spécifique (58) sous une 
pression nulle, et à la contraction, en grandeur 
absolue E Tr Aie 1) : 


(83) 
(84) 
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L’équation (84) admet la solution 


(VIVo)S = (86) 
c’est-à-dire : 


(86 bis). 


a 


i 
7 6 suivant les puis- 


sances entières et positives de pIT, on exprime Ja 
contraction à (85) par une série qui converge lente-- 
ment, même si le rapport He est minime (et 


diverge dès que le produit € + P. s'élève au-dessus 


T6 
de l’unité), (1) (2) 
2 0p, LONDRES PONS 
Dh (R) + 3! (#) : 41 
AU + QU +920 UE (00008 
TL n! (&) 


RE | 


Si, comparé à l’unité, le rapport p/T% est trés | 
petit, on a sensiblement : 


p.25 LP ÉCRIS 
SAS Lo à (r) 
Selon Bridgman, la contraction 5 (85) mesurée 1 


peut être exprimée : 
Ô — Ap — Bp? 


(88) | 


: 
(89) | 
À et B étant deux constantes ; d’où l’on COTE 1 
par comparaison avec l’égalité précédente (88), 1 a 
relation découverte par Slater à la suite d’un rai- : 
sonnement tout différent [23] 


Ua pe 


A2 


(90) | 


(1) Le module € de la dérivée (d Log.T”/9 Log. v)@ est : 
toujours positif. Ils’élève à quelques unités pour les cristaux : 
cubiques formés par les métaux (à 6 environ pour les cris- : 
taux de cuivre). 

(?) La solution (86) de l'équation (84) peut encore être 
formulée : 


p= lt 5) — 1). 


Et, en développant (1— 3) suivant les puissances posi- è 
tives ‘et entières de à on définit la pression p par une sérié, , 
qui est toujours convergente, et He la convergence est | 


plus rapide (à grandeur égale de Te ae s) RE || 
de la série (87) F5 
DEA Ê | 8 pl et Li (+2) 55 
3! 
J6 +2)... (C+n—1) 


11e Ne 42 


Mais la dérivée A ne reste pas exacte- 
d Log.v}, 


ment constante. Elle varie avec l’état du milieu 


+ cristallin, c’est-à-dire avec l’état auquel elle se rap- 
à porte ; donc elle varie avec la température et le 
: volume spécifique du milieu cristallin (1). Si elle 


n’augmente ou ne diminue, en grandeur absolue, 
que dans une proportion minime — de l’ordre du 
centième — sur toute l’étendue de la dilatation 


Dthermique, du zéro absolu de la température au 


point de fusion (69), elle varie sans doute nota- 
blement sous les fortes contractions que Bridgman 
fait subir aux cristaux métalliques. D’autre part, le 
rapport T@/l diminue quand la contraction à (85) 
“augmente, car d’après les relations (79) et (82), 
nous avons : 


Te Te —(p'— 
LE De PRES À A) ? 15 ? 
L: IPF : 2 ? ie 


Te et lo étant respectivement la tension ther- 
mique et le coefficient de compressibilité sous une 
pression nulle. En conséquence, la dérivée (80) 


(à Log.I”/d Log.v)e = — & 


ne reste pas invariable. De ce fait, les relations (85) 
et (86-87) ne sont qu’approximatives, et l’erreur 


» qu’elles comportent croît probablement avec la 
contraction à (85). Donc, quand la pression devient 


très forte (expériences de Bridgman), la contrac- 


à tion à ne peut être exprimée correctement par une 
» série telle que (88-89) réduite à 2 termes. 


Une formule établie par Einstein [24] implique 
légalité 


U[2 = — (d Log.« fd Log.v) + (1/6) (91) 


» et si l’on admet la relation (90), il vient : 


B[A? = — (d Log.wJd Log.v) + (2/3). 
D’autre part, selon Grüneisen (62) 


d Log.w PV 
“UF Hat 


(92) 


y; (93) 


| de sorte que l’on aurait, compte tenu de (91) : 


C2 = 7% + (1/6), 


} et compte tenu de (92) : 


(BA?) = Y + (213) 


(1) C’est que les dérivées (65) 
(à Log.C48 [do Log.v)o 


1 ne sont pas rigoureusement constantes. Elles varient avec 
| la température et le volume spécifique du milieu cristallin. 


Mais leur variation totale entre le zéro absolu de la tempé- 
rature et le point de fusion, quand le cristal se dilate libre 
de toute contrainte, reste très faible, De même, les dérivées 


(0Ca%8 00) (63) 


sont fonctions de la température et du volume spécifique. 
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Mais l'égalité (91) serait seulement satisfaite 
s’il n'existait pas des tensions thermiques, et si 
les modules de toutes les dérivées secondes C£8 (1-2) 
croissaient uniformément quand la contraction 
5 (85) augmente ; et leur croissance est loin d’être 
uniforme (67) (2). D’autre part la constante y (93) 
de Grüneisen est égale au rapport Ë/u (57), non au 
module de la dérivée moyenne 


d Log. w/d Log. v (62) (#). 


Partant, la relation (94) ne peut être satisfaite 
(sinon l’accord serait fortuit), et encore moins la 
relation (95), surtout si les coefficients À et B (89- 
92) se rapportent à des contractions amples pro- 
duites par de fortes pressions. 


Note sur les tensions thermiques 
et les énergies potentielles triatomiques 


Dans le mémoire intitulé « Les tensions ther- 
miques dans le milieu cristallin » [2], je rapporte 
chaque énergie potentielle «triatomique » w%?7 aux 
translations m + j — p — k, p + k — q — I, 
g +l—m—)j;etu, u?, uf étant les déplacements 
accomplis par les trois atomes en positions m + j, 
p + ket q + j, je développe en série l’énergie w%? 
suivant les puissances entières et positives des neuf 
déplacements relatifs 


mp% 


F0 pa pat pt a* am q* FAC: 
Ujx =U; —Uyx , Ugi = Ur —WÙ , Wjÿ = —u; 


Ph 

La série ainsi obtenue n’est pas déterminée entiè- 
rement. Car les trois translations m +j — p—Kk, 
p+k—q—1,q +1— m—)j, ont une résultante 
nulle, donc deux d’entre elles seulement sont indé- 
pendantes. 

D'autre part, si l’on rapporte chaque énergie tri- 
atomique w%?% à deux translations indépendantes, 
par exemple aux deux translations m + j — p — k, 


er a ANT 
m + j— q — I, les forces de rappel APRES LA 


finies par des expressions dissymétriques, ne se dé- 
duisant pas les unes des autres par la permutation 
circulaire des trois translations m + j, p + Kk, 


g+l. 


Je donne dans cette note de nouvelles expressions 

. M—PD MN—D,MN—A . 
des coefficients DIF AeBDISENN . Ces expressions 
restent symétriques par rapport aux trois trans- 


(?) Einstein suppose que toutes les oscillations ther- 
miques ont la même fréquence, déterminée par les coef- 
ficients d’élasticité (en particulier par le coefficient de com- 
pressibilité) et par la masse spécifique du milieu cristallin. 

1 
sr (Go). 

(*) Pour les cristaux de cuivre, le module de la dérivée 
d Log.w/d Log.v s’élève aux environs de 2,8, et selon 
Grüneisen y = 1,96 [11]. 


En outre, la relation (91) implique 
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lations m + j, p + k, q + I, et sont néanmoins 
entièrement déterminées. 
Afin de simplifier les formules, posons : 


mec) p+k=0, qi. 


Et soient u® ut ut les déplacements absolus des 
atomes a, b et c (en positions a, bet c)et 
ua — ya —— y : URL QU US Ut RUE 
leurs déplacements relatifs. 

Toute énergie potentielle triatomique w‘%* est la 


somme de trois termes : 
gpabe — ga be EL gyb—ca E pyc—ab 


w—t est l’énergie potentielle de l’atome a, dans le 
champ de forces commun aux deux atomes et Le 
wi—ta l’énergie potentielle de l’atome b, dans le 
champ commun aux atomes c et a, pen l’énergie 
potentielle de l’atome c, dans le champ commun 
aux atomes a et b. 

Rapportons l’énergie potentielle #7? aux deux 
translations indépendantes a — b, c — a. Dans la 
série de Taylor qui exprime cette énergie, les termes 
du troisième degré par rapport aux déplacements 


relatifs u%* u* ont pour somme : 
? ? 


1 
S3, a—bc = — ai à CDisUr AE uv avt ua 


es (»,4 se (e 4 
+ DS pes ut aa ac + DE ee uw ue get 


—D" C4 œ ÿ «—G œ Ÿ 
— “DS As 8 y ue us ue un De B y, a—bc u2 nca ue 
Et œ Ÿ C—€, a— œ Y 
— D? 6 Ce 2 UE y ne + D? a+ ao Ua us? ue 


= Dee me u° av u”] 
œ a 


et 
D Ds) D? ma be 
g B Y.a—bc d(az Enr bo) d(ab ÊE bb) d(aY — by) 
A d® wa bc 
De B y, a—be —  D(e — a°) d(cê LE af) d(CY — aY) 
Det RE d? wa be 
SHC dat — ba) D{aê — bB) D(cY — ar) 
“2 D? wa be 
Mere 


dc — ax) d(ai — DE) D(ar — by) 


Rapportons, de même, l’énergie w?—% aux deux 
translations b — c, a — b, donc aux déplacements 


relatifs uw”, ut, l’énergie w°—? aux deux trans- 
lations c — a et b — c, donc aux déplacements 


relatifs u, u, Et soient 8, 5, S3 at, les 
sommes des termes du troisième degré dans les 
séries de Taylor qui expriment respectivement les 
énergies w?—4 et w°—"%?, Posons : 

MT a—b , a—b TT 4—b 

De Y ‘= DE B y.a—bc Le 0, B Y, b—ca 

b—c ,a b—c TD b— 
Da B y = D, B Y.b—en + De By. o—ab 


Détate Die, c—ab + Des: a—bc" (1) 
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Ainsi, l’énergie triatomique w%% est définie par | 
une série de Taylor, où les termes du troisième degré | 


par rapport aux déplacements relatifs, u%%, ui, 
u%* forment une somme (de 21 X 27 termes) : 

Syate = Save + 93,0—ca + 93,e—ab 

= Fo (D, qu poule quuY ne | 


be œ 8 Y 
Se Dire Dr guv* ut uvY s Die ee a qu pc 


EE D 4 b. B uen u* rad + DE b. B uv ue ab 
—D., b= (e4 B Y bX (o] YŸ 
De db, C uv uvc u?° ! Pb, a UT Ua u°? ) 


Lg (DE 4 uv uvc$ uvY 


ess eee dupe voie done es à ns ele ste ee see obels ele pole +5 also nie 


So es ee ets e ess sol oo. le een © see etes sen elele aies Piote san 


se (D. d u°* ue ue 


ole e otre he se.ete ee » 6 er eee © se se ele eee oteles ee een 


eo 0e era te lee je + ha eee » 6/eNe 0 + ee » ee "eo olois se 0/5 'e se ete. s21e)s 


def pyuer)], 


Hormis les trois dérivées — D,#5t,°, — Dies 
eee (1) qui portent individueilement sur deux 
des trois énergies partielles 7%, wi, ge 


COCA IOCICION ICONE | 


toute autre dérivée qui figure dans la somme S3, aber | 
est définie par un symbole comportant deux indices ! 
qui rappellent les deux translations prises pour ! 


variables indépendantes, et partant, selon la con: 


No 42 | 


| 
| 
| 
| 
| 
| 


vention adoptée, qui indique l’énergie Nue | | 


dont elle provient. Par exemple, — 
rapportent aux deux translations a — b et b—6, 
c’est donc une dérivée de l’énergie partielle æ?—t : 
prb bo d° wbca 

BBOY TT Sfar — ba) o(al — 68) d(br — cv): 
Exprimons maintenant la somme S4, ave, en fonc: 


tion des neuf déplacements absolus u%, uw, u”. 
Soit S5, ave, l'expression obtenue. Puis développons 


directement l'énergie w‘* en série suivant les puis- 


sances entières et positives des mêmes dépla- 
cements absolus. Cette opération donne une nou- 
velle expression 85, ax, de la somme que forment 
les termes du troisième degré, par rapport aux dé- 


placements u% ul* u*,. Et, dans les deux sommes 


égales, S3. ave, et S5, ave, les Coefficients d’un même 


produit u‘* ut u* sont égaux. D’où en égalant les 


coefficients de la somme S'ux, aux coefficients | 


correspondants de la somme S” 4e, on obtient, en 
reprenant les translations m + j, p + k, q +1 
(respectivement égales à a, b et c) au lieu des for- 


Pre ns 
mules (10), les coefficients D, 74 #,et Di, FE a) rap 


portés aux énergies partielles we, qpé-da, pete, 


| 
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donc à des dérivées portant toutes sur deux dépla- 
cements relatifs qui sont indépendants 


Mm—D Mm— M—P, P—Q M—D, P—Q 
ji & — ( DESSIN AR Ve CNE Gel 
ARS 8 y 2 De B y DE B y de y 
+ 5 L'£TE + % î fs k 1) 
D D j Di 
œ À œ ê M COS" 
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ERRATUM 


I. Les facteurs de la dilatation. J. Physique Rad., Juillet 
1961, 22, 451. — À la page 457, dans la formule 29, au lieu 


de W lire W. 


Manuscrit reçu le 11 avril 1961. 
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REVUE DES LIVRES 


ABRAGAM (A.), Les principes du magnétisme nucléaire. 
(1 vol. relié de 599 p., 16 X 25 cm, éditions : Clarendon 
Press : Oxford University Press, 1961, en langue anglaise, 
prix 84 shillings.) 

Nul n’était mieux qualifié que M. Abragam, professeur 
au Collège de France, pour traiter ce sujet. D’abord, par 
son importante contribution personnelle. Ensuite (ou 
d’abord, peut être) par son exceptionnelle clarté. Que ce 
livre où l’on explique de quoi on parle, avant le dévelop- 
pement des formules, nous change de la littérature d’au- 
jourd’hui ! L’étude du magnétisme, ou plus précisément du 
paramagnétisme nucléaire, dû aux spins du noyau, est un 
sujet qui date des quinze dernières années. L’étendue de ce 
livre montre la rapidité de son développement. L'étude 
expérimentale est basée sur les phénomènes de résonance 
dans l’absorption des radiations électromagnétiques de 
basses fréquences (de mille mégacyles à deux kilocycles). La 
théorie est basée sur la mécanique quantique, mais non pas 
sur la théorie quantique des champs. Les applications prin- 
cipales concernent la physique nucléaire, l'évaluation des 
divers moments magnétiques nucléaires et la thermodyna- 
mique des très basses températures. Ce livre est donc indis- 
pensable à tout centre de physiciens et à tout chercheur 
qui aura le temps et Je loisir de parcourir et d'étudier ces 
quelques six cents pages. 


Voici un sommaire des matières : 1. Introduction. — 


II. Mouvement des spins libres. — III: Aspects macros- 
copiques du magnétisme nucléaire. — IV. Largeur de raie 
dipolaire dans un réseau solide. — V. Température des 


spins. — VI. Interactions électron-noyau. — VII. Structure 
fine des raies de résonance. Effets quadrupolaires. — 


VIII. Relaxation thermique dans les liquides et les gaz. — 
IX. Relaxation thermique et polarisation dynamique dans 
les solides. — X. La théorie de la largeur des raies spec- 
trales en présence du mouvement des spins. — XI. Struc- 
ture des multiplets de raies de résonance dans les liquides. 
— XII. Effets des champs de radiofréquence intenses dans 
les liquides et les solides. Ce livre comporte 15 planches 
photographiques hors-textes, en plus de nombreux gra- 
phiques. 
J. WINTER. 


Fermi (E.), Notes sur la mécanique quantique. (1 vol. 
broché de 171, p., 14 X 20 em, University of Chicago 
Press, 1961, en langue anglaise, prix Doll. 1,50.) 

On a reproduit les notes écrites qu’utilisa Fermi en 1954 
pour faire un cours de mécanique quantique, un an avant 
sa mort. Mais, si ses notes reproduisent bien les calculs et 
l’essentiel des raisonnements, les explications, les compa- 
raisons, ce qui ait l’originalité d’un cours oral, font défaut. 
Ceux qui pensent retrouver ici la clarté de l’admirabie livre 
« Molecole e Cristalli» (Bologne, 1934, Zanichelli) seront 
déçus. Il eût été plus précieux d’enregistrer le cours oral 
et de le reproduire, plutôt que d'utiliser sans retouche un 
aide-mémoire écrit. Malgré tout, ce livre intéressera tous 
ceux qui font de la mécanique quantique, l’enseignent ou 
l’appliquent. Il contient les bases de la théorie de Dirac, 
qu’on ne trouve généralement que dans des livres encom- 
brants et coûteux. Elles sont exposées d’une manière 
remarquable, et ne serait-ce que pour ces 25 pages finales, 
ce livre vaut d’être acheté. 

J. WINTER. 
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INFLUENCE DU p-Ag DE PRÉPARATION D’UNE ÉMULSION IONOGRAPHIQUE 
SUR LE DIAMÈTRE DES GRAINS 


Par F. SIMON (1), 


Centre de Recherches Nucléaires de Strasbourg, Département de Physique Corpusculaire. 


Résumé. — On expose les résultats obtenus dans la préparation d’émulsions ionographiques 
avec un excès d’ions Ag+ pour certaines conditions de précipitation. Les domaines de p-Ag déter- 
minés sont comparés aux données définies par d’autres travaux. On définit ainsi une bande de 
p-Ag la plus propice à la constitution de grains très fins et bien homogènes. 


Abstract. — We give the results obtained for the préparation of nuclear emulsions with excess 


of silver ions under certain conditions of precipitation. 
We define the most favorable values of p-Ag for obtaining 


red given in data with other works. 
very fine and homogeneous grains. 


Depuis quelques années, la préparation d’émul- 
sions nucléaires à grains ultra fins, c’est-à-dire dont 
le diamètre soit inférieur à 0,1 micron, est l’objet de 
nombreux travaux*à la suite des résultats de 
Perfilov [1] qui a pu produire des grains de bromure 
d’argent de très petites dimensions grâce à des 
conditions particulières de précipitation. Il appa- 
raîit de ces travaux que la présence d’un excès 
d'ions Ag* est indispensable à la formation de 
grains très fins, mais dans ce domaine, il faut tenir 
compte encore de la dispersion autour de la valeur 
moyenne qui dépend pour une grande part du pro- 
cessus de précipitation employé. 

Nous nous proposons d’exposer les résultats de 
nos expériences et de définir une région de p-Ag 
assurant la finesse des grains et leur homogénéité 
en confrontant nos données et celles d’autres 
auteurs. 


Mode de précipitation. — Dans un précédent 
article [2], nous avons détaillé le processus que nous 
employons pour préparer nos émulsions concen- 
trées. Nous avons quelque peu modifié ce procédé 
depuis en assurant d’une part une régulation auto- 
matique de la concentration des ions Ag*, d’autre 
part, en augmentant le volume d’un lot qui a été 
quadruplé. De plus, la durée de toute la précipi- 
tation a été réduite à 15 minutes. La gélatine que 
nous employons est une gélatine (?) de type inerte 
Rousselot dont la viscosité Engler à 10 % est, à 
40 0C, de 209, la teneur en chlorure de 1,150/60 en HCI. 
En fin de précipitation, l’émulsion ne subit aucune 
maturation et est alors ajustée à un p-Ag de 8,5. 

La mesure de la concentration des ions Ag+ est 
réalisée au moyen d’une électrode d’argent et d’une 
électrode de calomel saturé. Celle-ci se trouve 


(1) Adresse actvelle : Faculté des Sciences de Rabat 
(Maroc). Laboratoire de Physique. 

(2) Nous remercions la Cie Centrale Rousselot de nous 
avoir fourni très aimablement cette gélatine et son analyse. 


The values of p-Ag determined are compa- 
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Fic. 1. — Distribution du diamètre des graius de BrAgen 
fonction du p-Ag lors de leur formation. 
Courbe 1 : p-Ag = 3,8 
» : 


» 
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RU 


reliée au milieu à contrôler par un pont de KCI 
saturé et un pont de KNO; N/10. Tout au long 
de la précipitation, l’excès d’ions Ag est main- 
tenu constant. Ceci explique peut-être que nous 
v’ayons pu travailler à des p-Ag très bas sans 
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obtenir au développement de voile très fort contrai- 
rement à ce que Perfilov a pu faire. 

Le domaine des p-Ag dans lequel nous avons 
opéré s'étend de 3,8 à 6,2, ce qui correspond pour 
cette dernière valeur à un excès d’ions Br—, puisque 
le p-Ag isoélectrique est de 5,8 à 50 0C, température 
à laquelle se fait la précipitation du bromure 
d'argent de nos émulsions. De part et d’autre de la 
bande de p-Ag que nous avons étudié, des essais 
de préparation d’émulsions ont cependant été faits, 
mais avec des résultats peu favorables, soit que la 
dimension des grains soit trop forte, ou leur dis- 
persion trop étendue ou le voile trop important. 
Dans le cas où la précipitation s’effectue à un p-Ag 
inférieur à 3,8, nous obtenons une forte montée du 
voile sans bénéfice pour la dimension des grains. 

Si au contraire, les grains de bromure d’argent se 


… Fic. 2, — Microphotographie électronique de grains de bro- 
mure d’argent précipités à un p-Ag de 4,8. — Grossisse- 
ment : 20.000. 


forment dans un milieu dont le p-Ag est supérieur 
à 6,5, nous voyons une augmentation de la dis- 
persion du diamètre des grains apparaître et la 
valeur moyenne de ce diamètre croître pour des 
valeurs croissantes de la concentration en ions 
brome. 

Pour les valeurs de p-Ag comprises entre les 
limites précisées plus haut, nous avons obtenu les 
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meilleurs résultats tant en ce qui concerne le dia- 
mètre des grains que la faiblesse du voile. Les 
courbes de la figure 1 sont les histogrammes faits 
sur cent grains du diamètre mesuré sur des agran- 
dissements à 20 000 de microphotographies électro- 
niques de répliques de l’émulsion très diluée. La 
figure 2 est une partie de l’un de ces agrandis- 
sements obtenu avec une émulsion préparée en 
maintenant son p-Ag égal à 4,8. 

La dépendance entre la concentration en 1ons 
Ag* et le diamètre des grains est surtout marquée 
par rapport à la dispersion qui est minima pour les 
p-Ag les plus proches du p-Ag isoélectrique, le dia- 
mètre le plus faible avec la moindre dispersion 
étant provoqué lorsqu'on se trouve en présence 
d’un excès d’ions Ag*. Nous avons représenté sur 
la figure 3 la variation du diamètre avec le p-Ag ; 
elle montre que le grain est d’autant plus petit que 
la concentration des ions Ag+ est plus élevée. Le 
tableau suivant donne pour différents p-Ag les dia- 
mètres moyens mesurés comme dit plus haut : 


DAS #0 270 GEL ESA 3 
CAP C RNA CR PCR EL ET 


Les diamètres étant exprimés en 10—2? microns. 
Nous avons pensé qu’il pouvait être intéressant de 


Diamètre en 1 * 10 ? 


5 


Fic. 3. — Variation du diamètre des grains en fonction 
du P-Ag. 
Courbe 1 : gélatine Rousselot 12.805. 
Courbe 2 : » Keystone 2191 d’après Oliver. 
Courbe 3 : » Winterthur 5444  » » 


confronter ces données avec celles des différents 
auteurs qui se sont occupés de préparer des émul- 
sions à forte concentration en bromure d’argent et 
à grains très fins. En effet, l'influence des autres 
facteurs entrant dans la formation du grain n’est 
pas négligeable et en particulier le type de gélatine 
utilisée. 
02 
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Parmi les émulsionneurs qui se sont attachés à 
ces problèmes, on doit citer, outre Perfilov, Oliver 
[3], Fourneaux, Narath et Heimann [4], Kubal [5], 
Gauvin et Sebaoun [6], ces derniers ayant étudié 
la précipitation du chlorure d’argent pour des p-Ag 
supérieurs au p-Ag isoélectrique. 

Voici quels sont les diamètres les plus probables 
obtenus par ces auteurs suivant la concentration 
en ions Ag+ de leur milieu de précipitation : 

Oliver avec gélatine Keystone 2191 : 


p-Ag : 6 4,9 k 2,9 
Qu : 0,42 0,06 0,07 . 0,10 
Avec gélatine Wintherthur 5444 : 
P-Ag 6,7 DAS) 3,8 RUE 2 
Ou DÉMO RTTNO PE 0:08 0,06 
Heimann (gélatine D. G. F. 9457) : 

p-Ag : 0,08 0,10 
Ou : 0,08 0,10 


Kubal, sans donner de valeur de p-Ag donne les 
valeurs : 
NC SR 
Ou : 0,19 0,11 


<< 5,8 (si la température est 500 C) 
0,08 


Enfin Perfilov, en faisant varier par bonds le 
p-Ag de 1,56 à 7,8 réussit à préparer des grains de 
diamètre 0,06 microns. 

En rapprochant ces données des résultats que 


N°12. 


nous avons obtenu, on trouve bien confirmée la 
nécessité d'opérer à des p-Ag inférieurs au p-Ag 
isoélectrique quelle que soit la gélatine employée si 
l’on désire préparer des grains très fins de bromure. 
d'argent ; la dispersion des diamètres peut aussi 
être liée à la concentration des ions Ag*, et semble 
être meilleure si l’on est plus près du point isoélec- 
trique. En effet, Kubal, à l’équivalence trouve des 
grains ayant de 0,04 à 0,18 microns de diamètre 
avec un excès d’argent de 0,05 à 0,13 microns, 
tandis qu'avec un excès d’ions bromure les dia- 
mètres sont possibles de 0,06 à 0,32 microns. 

Oliver trouve de même un optimum de la dis- 
persion aux p-Ag 4,0 à 4,9 avec la gélatine Keystone 
ainsi que nous le remarquons dans nos essais ;, 
cependant s’il utilise la gélatine Winterthur, la 
moindre dispersion se produit pour les bas p-Ag 
de 2,2 à 3,5. 

Sans doute d’autres facteurs entrent en ligne de 


compte sans paraître toutefois aussi importants ; | 


en effet, nous avons pu changer les volumes réagis- 


-sants et la durée de précipitation d’un lot sans 


remarquer de différence dans le diamètre des grains 
et sa dispersion. C’est aussi ce qui a lieu entre les 
divers procédés que nous venons de comparer. Un” 
facteur est encore mal connu : la variation du p-Ag 
périodique au cours d’une précipitation, au moins 
en ce qui concerne la dimension des grains. 


Manuscrit reçu le 6 juillet 1961. 
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LES CHAMPS DE DOUBLET DE MASSE 
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Résumé. — Les champs symétriquement self-couplés peuvent être définis pour tous les spins 
vis-à-vis de toutes les huit inversions unitaires. De tels champs sont également réalisables vis-à-vis 
des huit inversions antiunitaires en cas de spins entiers, mais ne le sont que vis-à-vis de Cet PTM 
dans le cas des spins demi-entiers. Chacun de ces champs comporte deux états de masse 
propre différente. Des arguments d’ordre physique excluent la réalité physique de tels champs 
définis vis-à-vis des huit inversions indépendantes de M, mais ceux définis vis-à-vis des huit inver- 
sions comportant M sont bien admissibles au point de vue physique. Cependant des systèmes phy- 
siques particuliers et bien spécifiés peuvent bien être représentés par certains des champs self- 
couplés définis vis-à-vis même des inversions qui ne comportent pas M. Dans le cadre de cette 
théorie, le champ fermionique M-symétriquement self-couplé est proposé pour les doublets-M 
leptoniques, électron-muon et neutrino, et ie champ bosonique C (ou CP)-symétriquement self- 
couplé représente le doublet-C (ou CP) mésonique Æ9-K$. 


Abstract. — Symmetrically self-coupled fields can be defined for all spins with respect to all 
eight unitary inversions. Such fields are also realizable with respect to the eight antiunitary 
inversions in the case of integer spins, but are realizable only with respect to C and PTM in the 
case of half-odd-integer spins. Each such field has two states of different rest-mass. Some 
physical arguments exclude the physical reality of all such fields defined with respect to the eight 
M-independant inversions, while those defined with respect to the eight M-dependent inversions 
are quite admissible in the physical point of view. However, some particular and well specified 
physical systems can quite well be represented by certain of the self-coupled fields which are 
defined even with respect to the M-independent inversions. In the framework of this theory, 
the M-symmetrically self-coupled fermion field is proposed for the leptonie M-doublets,electron- 
muon and neutrino, and the C (or CP)-symmetrically self-coupled boson field represents the meso- 
nic C (or CP)-doublet K?-XK%. 


1. Introduetion. — Dans une étude en vue d’une 
théorie unifiée de l’électron et du muon [1], nous 
avons proposé et examiné un champ de spin 1/2 M- 
symétriquement self-couplé, à savoir, un champ 
spinoriel Ÿ(m) qui est mis en interaction avec son 
«image d’inversion M » d(— m). m désigne ici 
masses propres et M l’opération d’inversion de 
masses propres. En général, des champs de ce 
genre existent bien vis-à-vis de toutes les opérations 
d’inversion du groupe d’inversion 


G={1, P, CT, CPT, M, PM, CTM, CPTM 
C,CP,T, PT, CM, CPM, TM, PTM } 


En d’autre termes, un champ, en général, peut 
être mis de façon compatible en interaction avec 
son image vis-à-vis de chacune des seize inversions. 
De tels champs comportent chacun deux masses 
propres identifiables à deux particules différentes. 
On dira que celles-ci forment un « doublet-G ». 
Comme nous avons expliqué récemment [1], les 

champs que nous étudierons obéissent à des équa- 
tions du type 


Ba dx Ÿ + ma d = m[G]Y. (1) 


Ici m, et m, sont des constantes réelles que nous 
1 AREAS 
supposerons positives toutes les deux. Dans une 


étude générale qui nous occupe maintenant, nous : 


prendrons en considération tous les spins dans ce 
qui suit. Dans un travail antérieur [2], nous avons 


étudié certaines propriétés de symétrie intrinsèques 
vis-à-vis du groupe d’inversion G dans le cadre 
d’une large classe générale des théories des parti- 
cules à spin. Dans le présent exposé, nous nous 
référerons aussi à cette classe générale des théories 
des particules et nous aurons l’occasion de revenir 
à quelques unes de leurs propriétés de symétrie 
intrimsèques. Les matrices 8; obéissent donc aux 
règles non commutatives que nous avons résumées 
dans le travail antérieur mentionné plus haut [2]. 
Dans l’équation (1), [G] 4 désigne le champ engen- 
dré de Ÿ par l’inversion G[3] avec encore un chan- 
gement de nombres quantiques dans ce champ 
engendré, changement impliqué dans la définition 
et dans la caractéristique de l’inversion G[31. Il y 
a donc à distinguer l’opération [ G] de l’opération 
(G) que nous avons définie et étudiée antérieure- 
ment [3] (1). 

À part l’équation de champ (1), notre problème 
des champs symétriquement self-couplés consiste 
à poser encore deux conditions sur l’opération [ G]. 
(z) Le champ engendré de 4, [G] Ÿ, représente 
aussi le même champ bien que non identique à . 
Aîinsi [G] d doit vérifier l'équation 


Ba LG] Ÿ + mlG] 4 = mA[ GI ÿ. (2) 
(:) Les deux opérations [G] et (G) sont identiques dans 


le cas G = M et dans ce cas seulement, car (M) n’engendre 
pas de changement de nombres quantiques. 
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(it) Puisque la présence de [G] 4 dans l’équa- 
tion (1) signifie un certain self-couplage par « image 
d’inversion G », 4 doit apparaître lui aussi comme 
un tel couplage dans l’équation de [G] Ÿ, ce qui 
implique une symétrie par l’échange de 4 et [G] Q. 
Dans ce cas, on doit pouvoir écrire l’équation sui- 
vante : 


Ba d1[G] Ÿ + mA[G] Ÿ = ma 9. (3) 


Des équations (2) et (3), une condition de compa- 
tibilité s'impose dans la définition des champs 
symétriquement self-couplés : 


[GR y = y. (a) 


La question principale est alors la détermimation 
des inversions possibles vis-à-vis desquelles peu- 
vent être définis de façon compatible les champs 
symétriquement self-couplés. 

Des équations (1) et (3), on déduit les champs 
non couplés, 

Ba dde + Ma Ve — 0, 


Bx dx Ve + Mme Lo = 0, (5) 
où l’on définit 
Ma = Mi Me, My — M cn Ma, (6) 
1 FL 
= — (1 G ; = —4—1[G 5 7 
Ÿ ven ele dr ei [G]) à. (? 


Le champ d comporte donc deux masses propres, 
m et M, qui sont identifiables à deux particules 
différentes formant un doublet-G. 

Étant donné que Ÿ et [G] d satisfont, par défi- 
nition, à la même équation (1) et sont mêmes mis 
en interaction, 1l est naturel que tous les deux ont 
des propriétés de transformation identiques. Les 
champs symétriquement self-couplés sont alors 
tous invariants par rapport à toutes les inversions 
et au groupe de Lorentz. 

Dans l’examination des inversions qui permet- 
tent de définir de façon compatible des champs 
symétriquement self-couplés, nous devrons bien 
distinguer les inversions unitaires des inversions 
antiunitaires. À ce sujet, nous les examinerons 
séparément aux paragraphes suivants. 


2. Inversions unitaires. — Désignons par À l’une 
quelconque des huit inversions unitaires : À = 1, 
PECT, CPT;M, PM, CTM, CPTMALEIcRhAMp 
engendré de d par À est alors défini par (voir [3]) 


[A] Ÿ = Ca a 4, 


où l'indice inférieur À dans Ÿ, désigne un change- 
ment de coordonnées (e, r, t, m) et un changement 
de nombres quantiques dans d, changements 
impliqués par la définition de À comme nous 
l'avons prescrit au paragraphe précédent. Les l, 
sont des opérateurs matriciels non singuliers qui 
remplissent les conditions, [2], [3] 


Tu Bla! = 89 B, (9) 


(C4, constante) (8) 
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où 3% a pour valeur -+ 1 ou — 1 suivant À et à. 
Nous pouvons prendre pour ['; des matrices à la 
fois unitaires et hermitiennes pour toute inversion 
HIMbAIRE AE ME EN? 

D'une seconde opération de [A] sur (8), nous 
obtenons 


[4]? d = C4 Ÿ. 


La condition de compatibilité (4) sera donc réa- 
lisée si 


(10) 


Ci =1, | 
cela indépendamment de À. 

Vis-à-vis de chaque inversion unitaire À, nous 
avons donc deux champs symétriquement self- 
couplés définis par C4 — + 1 et — 1 respective- 
ment. Dans le premier cas (€4 — 1), selon (7), les 
membres non couplés de masse propre bien définie 
sont 


La = (4 Di Ua de = (4— T'4 da), (12) 


qui satisfont aux conditions de symétrie suivantes: 


IN” doA = = Va; 3 UpA Sir Vy- 
Au second cas (€4 — —- 1), les membres non 
couplés sont 
dé = (b— Tata, = Et +Tava (19 
FOIE ST= Re A YA); MERS — A YA}; 
CU v2 | 
qui vérifient les conditions de symétrie 
Tadoa = — da Ta Ya = + dr. (15) 


Dans les deux cas, les champs a sont tous l’état. 


propre de la petite masse m, — m, et les champs b 
celui de la grande masse m, + m,. Mais, d’après 
(13) et (15), ils ont des propriétés de symétrie diffé- 
rentes suivant que {4 = + 1 ou — 1. Nous dirons 
que les champs d4 et 4, forment un doublet-A et 
d, et d, forment un doublet-4’. (C£. [1]). 


3. Inversions antiunitaires. — Soit B l’une quel- 


conque des huit inversions antinuitaires : B=€C, 
CP,T, PT, CM, CPM,TM, PTM. Dans cetcss 
par définition, [3] 


[8] V= 68 le V8, 

L'indice supérieur * désigne la transposition des 
matrices. 4 est le champ adjoint de à, à = dt 8, 
0 étant une matrice non singulière définie par 
B0 = — 68, B,6 — 08, 06000 


matrices non singulières l'} doivent vérifier les 
conditions [2], [3] 


(Ëg, constante). (16) 


De 6x ls = 6, (17) 


où également 39 — + 1 ou — 1 suivant B et À. 
Nous prendrons pour l'; des matrices unitaires 


N° 12 


(11) 


as 


_ No12 


- et hermitiennes en cas de toute inversion antiuni- 
Re Drm =r, 

Nous avons démontré [2] que, indépendamment 
de la représentation des matrices f;, les opérateurs 
. L'4 comportent les propriétés de symétrie intrin- 
sèques suivantes : 


T 
TB = ep l, 


Te 0715! = 280, 


= 1; (18) 


a — 1. (19) 


sg et «x prennent pour valeur + { ou — 1 sui- 
vant l’inversion B et le spin considérés. Des consi- 
dérations générales nous ont permis de déterminer 
les valeurs de € et «4 en fonction du spin en cas de 
chaque inversion antiunitaire B [2]. 

Cela étant, une seconde opération [B] sur d nous 
donnera, selon (16), 


[B1? ÿ = Esp os ep Ÿ. (20) 


Pour que les champs symétriquement self- 
couplés soient réalisables vis-à-vis de B, il faut donc 
que 


rte 
EméR= up er =)1. (21) 


Mais, des résultats d’une étude antérieure men- 
tionnée plus haut [2], nous tirons facilement 


(22) 


en cas de Cet PT M et indépendamment du spin 
considéré, et 
XB EB — (— 1)?5, 


DRE 


(S, spin) (23) 


en cas de chacune des six autres inversions anti- 
unitaires. On peut donc conclure que le champ 
symétriquement self-couplé peut être défini de 
façon compatible vis-à-vis de toutes les inversions 
antiunitaires en cas de spins entiers. Mais 
dans le cas des spins demi-entiers, un tel champ 
n’est réalisable que vis-à-vis de C ou de PTM. 
Dans le cas des inversions antiunitaires, quand 
les champs symétriquement self-couplés sont réa- 
lisables, les membres non couplés de masse propre 
bien définie sont 


x 1 : = 
d, = (4 +Esls V8), Ÿp ANT Ep LB VB), (24) 


qui remplissent, selon (21), les conditions de symé- 


trie suivantes : 


qe = LT Cnil c 
TB Vas = EB° Vo, Le Vos = — Ep V4. (25) 


4. Ondes planes monochromatiques. Significa- 
tion physique des conditions de symétrie. — Les 
champs symétriquement self-couplés vis-à-vis du 
groupe d’inversion G sont des champs linéaires. 
Chacun d’eux comporte deux masses propres, 
Ma = M — M2 et Mo = M, + M, Tout comme 
dans le cas du champ self-couplé vis-à-vis de MTL], 
on trouve facilement les solutions de ces champs en 
ondes planes monochromatiques. À une impulsion 
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p correspondent deux énergies Æ, et Æ, données par 
Eÿ = p? + mé. 


Représentons les deux sortes 
monochromatiques par 


Ei = pi + m, (26) 


d’ondes planes 


(s) 
a '(PEas rt, + m) = up, E, ; 


(s) ïi ps. 
bb, Ev;rt, + m) = uÿ”(p, E,:; + mjePb-"—Et), 


FE 


m) eYP-"—Eat), 


(27) 


où s désigne la valeur propre du spin longitudinal. 
Les spineurs u, et u, obéissent aux équations 


(iB.p— Ba Ea + ma) U/(p, Eu; + m) = 0, 


(iB.p— Ba Er + mu) uÿ(p, E5; + m) = 0. 


Ce qui caractérise la théorie des champs symé- 
triquement self-couplés, ce sont les conditions de 
symétrie (13), (15) et (25), que doivent remplir les 
champs non couplés d et d, qui sont les deux états 
de masse propre bien définie. Dans le cas des ondes 
planes, ces symétries introduisent des conditions 
sur les spineurs u, et u,. Nous résumons ces condi- 
tions de symétrie explicitement ci-dessous pour 
chacune des seize inversions du groupe G. 


(28) 


Inversions unitaires. 
(CA) 
L'est un cas trivial. 
TP ug (bp, Ea; + m) = + ur (—p, Es; 
(s 
b 


ro) 
Do up, Eb 3 + m) =— df—p, y; + m); | 
(29) 
l'or ua /(P, Ea; + m) = + ug (b,— Es; + m), 
l'or uÿ/(p, 5: + m) =— uÿ (p, — Ey 5 ÆLTAL: 


. 
L'opr u${p, E: + m) = + uÿ(—p,—E,; + m), 


mn) UD Erree 


Topr uÿ/(p, Ev; + m) =— uw (—p,—E;+m); 
(31) 
Dar 8) (b, Eos + m) = Æ ug'(p, Eas Fm), 
INT us) bp, EE; + m) = — uS(p, ES - TR) | 
(32) 
L'or 0 (P, Bas + m) = + ug (—P, Ea; Fm), 
pm Up, E5: + m) =— uÿ(— p,E;; F.mj;: | 
(33) 
l'orm u/(p, Ea; + m) = + us (pb, — Es; F m), 
l'orm uÿ(p, Ev; + m) =— uÿ (pb, — Er; Fm); | 
(34) 
l'oprm up, Eu: + m) = + ua (—b,— EE; rl 
Loprm uÿ (bp, Ev3 + M F m). 
(5) 
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Inversions antiunitaires. 
(s A Ep = ei?g) 


lou" (b, Eu; + m) 


Ds (36) 
ouf” (p, EE; + m) 


+ 
(ae) 
ea: 
+ 
Q 
SG 
| 
= 
= 
ue 
ÿ : 
== 
me — 


JE 
Topus (p, Ea; + m) 


mn (37) 
Top u)"(p, Es; + m) 


Ÿ 

li) 

a) 
SJ 
S 
La 
+ 


m), 


Tru” (p, E,; + m) 


HR Lo eur 


H 


de 38 
ai (38) 


Tr uÿ p, Es; E m) 


re 
=— € Ty (— PRES ILE TR) 


DAT 
T'pr uÿ) (bp: EG; +m) 


tro MU Dem) . 
Fes : 9 
Dpru (pb, 1 0) 


” : 
— e PT (p, Ep: + m); 


ee ? 


CE  “ 


ST 
Tomue”. (b, Ea; + m) 
ip (s) 2: 
= Le CMyE HE + mn), 40) 
+ 
T (s)T THE 
cmub (pb, Er; + m) 


= e Mu be Er +m); 


on P 
T'opmus (p, Eu: + m) 


An PNR Co 
(41) 


T 
T'opmuy (pb, E5: + m) 


iecrm (5) + 
—e Mug (p,— Er; + m); 
(ST 


Trmua (p, Es; + m) 


ip ne (5) à Æ 
SR <F e Lure (— p, E, US m), 


ee | 
Tomuÿ” (bp, E5; + m) 


ip (3) Ft 
—e Mur (— p, En: 


? 


= 

Z 

2 
TZ 


_(9T 
TrrM ua 


(b, Ea; + m) 
= + e *PTM 9 (p, DEEE 
Tru (b, Ep; + m) 
— — e PTM y (p, Ep, + m). 
Ces conditions de symétrie sur les spineurs w et 
u, en théorie des champs symétriquement self- 


couplés coïncident exactement avec les normali- 
sations uniformes des ondes planes monochroma- 
tiques vis-à-vis du groupe d’inversion G, normali- 
sations que nous avons étudiées en détail dans le 
cas du spin 1/2 dans un article récent [4]. Mais, 
à la différence des théories usuelles, les champs 4, 
dans la présente théorie, ne peuvent être unifor- 
mément normalisés que selon une seule possibilité 


et les champs b doivent l’être selon l’autre possi- 


bilité différente. 


Dans le cas des demi-entiers, on a, 


spins 


selon (23), ag eg = — 1 pour B Æ C ou PTMS 


Les normalisations uniformes des ondes planes ne 
sont done pas réalisables vis-à-vis de ces six inver- 
sions antiunitaires en cas de spins demi-entiers, 
comme nous l’avons d’ailleurs montré en dé- 
tail antérieurement dans le cas du spin 1/2 [4]. 
Par conséquent, les conditions de symétrie 


(37) — (42) ne sont pas remplies en cas de spins. 


demi-entiers. Dans le cas des spins entiers, au 
contraire, les normalisations uniformes des ondes 
planes sont réalisables vis-à-vis de toutes les 
seize opérations du groupe d’inversion G, et les 
conditions de symétrie (29) — (43) peuvent toutes 
être bien réalisées. 


Cela confirme, d’une façon explicite, les conclu- 


sions des deux paragraphes précédents. 

Les conditions de symétrie (29) — (43) sur les 
ondes planes des champs symétriquement self- 
couplés impliquent des limitations physiques sou- 
vent inadmissibles sur la structure même de ces 
champs. Considérons par exemple le champ 2P- 
symétriquement self-couplé. Les conditions de 
symétrie (29) doivent être réalisées sur les ondes 
planes de ce champ. Elles signifient (cf. [4]) que, 
dans le champ a comme dans le champ b, s’il ya 
une onde plane u®(p, Ei), il y aura aussi et néces- 
sairement l’onde plane correspondante u® (— p, 
E;) qui a d’ailleurs, à un signe près, une même am- 
plitude que uÿ (p, Æ;). Cette condition est, au 
point de vue physique, incompatible avec les 
conceptions classiques, car par exemple on peut 
bien avoir un faisceau d’électrons d’une quantité 
de mouvement p sans avoir dans le même faisceau 


un nombre égal d'électrons ayant une quantité de - 


mouvement — p. Le champP-symétriquement self- 
couplé ne représente donc pas la réalité physique 
générale. 

Des arguments d’ordre physique semblables 
excluent en fait la réalité possible de tous les 
champs symétriquement self-couplés vis-à-vis des 
huit inversions indépendantes de M. Mais de tels 
arguments ne s'appliquent pas à ceux définis vis- 
a-vis des huit inversions comportant HW. Dans ce 
cas, en effet, les conditions de symétrie n’intro- 
duisent des relations qu’entre les champs d(m) et 
d(— m). Ces conditions n’apportent donc pas de 
modification aux théories usuelles des particules 
qui ne sont basées que sur le champ V(m) sans 
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tenir compte de Ÿ(— m) dans une manière symé- 
trique. Les champs self-couplés définis vis-à-vis 
des inversions comportant M, qui tiennent compte 
de la symétrie de + m et — m, sont donc tous 
admissibles au point de vue physique. Le plus 
simple de ces champs est celui défini vis-à-vis de 
l’inversion M elle-même, qui est le champ fonda- 
mental d’un modèle des leptons que nous avons 
proposé et étudié récemment [1]. 

Cependant, si nous nous bornons à des états phy- 
siques particuliers et bien spécifiés, les chiamps 
symétriquement self-couplés vis-à-vis même des 
inversions ne comportant pas M pourront bien 
représenter chacun un système ou un problème 
physique particulier. Ainsi, au cas où l’on ne consi- 
dère que des états physiques qui comportent par 
exemple les états (p, Es, s) et (p, — Es, — s) d’une 

. même amplitude selon les conditions de symétrie 
(37), le champ symétriquement self-couplé vis-à- 
vis de CP pourra bien représenter ces états phy- 
siques particuliers. Nous en avons un exemple : les 
mésons À neutres dans l’interaction faible de 
désintégration. Dans cette interaction, en effet, les 
mésons X© directement observables sont les deux 
états propres de CP, K? et K%, définis par (voir 
(51) 


1 1 
D —— 4: CP]\) K9 — — {XK0 Ec ne K@ , 
, | + [CP]) FL + Ccr lcr 7e) | 
1 1 F0 
K9 —— (1 — [CP]) K° = (A0 — Ecp l'or Kcr). 
5 À [CP1) Sal cP 1cP 52) | 


(4) 


K2 et K2 obéissent dans ce cas aux conditions 
de symétrie (25) et (37). Nous dirons qu'ils for- 
ment un doublet-CP. Ils peuvent bien avoir des 
masses propres différentes dans notre formulation 
des champs symétriquement self-couplés. Un tel 
champ est d’ailleurs réalisable, car les X° sont des 
bosons. 

Si l’on définissait K? et KX2 comme les deux états 
propres de CPM (doublet-CPM) au lieu de CP, le 
problème des mésons Æ0 ne serait plus limité aux 
seuls états physiques particuliers et bien spécifiés 
que nous avons prescrits plus haut. 


5. Quelques cas spéciaux. — Nous examinons 

maintenant quelques cas spéciaux des champs 
symétriquement self-couplés en ce qui concerne les 
masses propres. 
40 m, — 0. Les deux membres du doublet-G 
auront la même masse propre m,. Les champs fon- 
damentaux d ne seront plus self-couplés et seront 
simplement les champs usuels. En fait, la structure 
de doublet de ces champs perd sa signification 
dans ce cas. : 

20 m, = 0. Les champs considérés sont toujours 
self-couplés dans ce cas, et on a toujours de vrais 
doublets, bien que leurs deux membres aient main- 
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tenant la même masse propre m,. Les équations de 
champ (1) et (3) s’écrivent maintenant 


Ba d1 Ÿ— m[G] Ÿ = 0, 
Ba DL GT Ÿ— m3 Ÿ = 0. 


Les deux membres du doublet-G obéissent alors 
aux équations 


Ba dx Va + Ma Va = Fi 
Ba dx Vo — Ma Vo = 0, 


da et dy étant définis toujours par (7). Il est à 
remarquer que le membre d, a pour terme de 
MASSE — Mo. 

Il serait intéressant de considérer A9 et K2 
comme un doublet-CPM dans le modèle (45), ce 
qui équivaut au doublet-CP dans la théorie usuelle. 
Nous avons déjà signalé au paragraphe précédent 
l’avantage de ce nouveau point de vue. 

30 Mi = Ma = 0. Les doublets perdent leur 
signification. Dans le cas G = M et du spin 1/2, 
ce sera le champ du neutrino à quatre composan- 
tes [1]. En cas de spin 1, on aura le champ du 
photon. 


(45) 


(46) 


6. Conclusion et remarques. — Les champs symé- 
triquement self-couplés peuvent être définis de 
façon compatible vis-à-vis de toutes les inversions 
unitaires et indépendamment du spin. De tels 
champs sont réalisables aussi vis-à-vis de toutes les 
inversions antiunitaires en cas de spins entiers, 
mais ne le sont que vis-à-vis de C'et de PT M dans 
le cas des spins demi-entiers. 

Des arguments d’ordre physique excluent la 
réalité physique générale des champs self-couplés 
vis-à-vis des huit inversions indépendantes de M. 
Mais ceux définis vis-à-vis des huit inversions 
comportant M sont bien admissibles au point de 
vue physique. De ceux-ci, le champ le plus simple 
est celui défini vis-à-vis de M elle-même, qui, dans 
le cas du spin 1/2, est le champ fondamental du 
modèle de doublets-WZ, un modèle des leptons que 
nous avons étudié récemment [1]. 

Cependant, si l’on se borne à des états physiques 
particuliers et bien spécifiés, les champs symétri- 
quement self-couplés définis vis-à-vis même des 
inversions indépendantes de M peuvent bien repré- 
senter chacun un système ou un problème phy- 
sique particulier. Dans ce sens, A? et A5 forment 
un doublet-CP. Mais on pourrait définir aussi X® 
et 5 comme un doublet-CP MW, de sorte que le pro- 
blème des mésons Æ neutres ne soit plus limité aux 
seuls états particuliers mentionnés plus haut. 

On pourrait penser à une généralisation possible 
de la théorie des champs symétriquement self- 
couplés en y comprenant aussi ceux définis vis- 
à-vis des rotations À ou des transformations de, 
Lorentz L. On aurait alors un type de champs qui 
peuvent comporter chacun un nombre quelconque 
de masses propres différentes au lieu de deux seu- 
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lement dans le cas des champs self-couplés vis-à- 
vis des inversions. Mais des considérations d’ordre 
physique comme celles évoquées au paragraphe 
précédent excluent en général la réalité physique 
de tout ce type de champs. On ne pourrait pour- 
tant pas exclure la possibilité qu’ils puissent repré- 


senter certains problèmes physiques particuliers et 
bien spécifiés. En outre, les champs self-couplés 
définis vis-à-vis de RM ou de LM seraient bien 
admissibles au point de vue physique. 


Manuscrit reçu le 23 décembre 1960. 
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REVUE DES LIVRES 


Puex (E. M.) et Puex (E. W.), Principes of electricity and 
magnetism. (1 vol. relié toile, 430 pages, 1960. Édi- 
teur Addison-Wexley, Reading  (Massachusset), 
8,75 dollars (en anglais)). 

D’allure essentiellement didactique, ce livre reprend les 
principes fondamentaux de l'électricité et du magnétisme. 

Après une introduction sur les systèmes d’unités et les 
notions d’analyse vectorielle et de calcul intégral néces- 
saires, les auteurs traitent des théorèmes fondamentaux de 
l’électrostatique à l’aide des opérateurs différentiels. La 
première relation de Maxwell est introduite dès le début de 
l’étude des diélectriques. 

Les auteurs étudient ensuite les lois classiques du courant 
continu, sous l’angle des vues modernes sur les conducteurs 
métalliques. 

Un important chapitre est consacré à l’induction électro- 
magnétique et aux deuxième et troisième relations de 
Maxwell. 

. La troisième partie décrit d’abord le champ magnétique 

d’un courant (quatrième loi de Maxwell), puis les corps 

magnétiques : quelques pages traitent notamment de l’im- 

portante question des électrets. On peut cependant regretter 

que le paramagnétisme et le diamagnétisme soient insufl- 
samment expliqués. 

Les deux derniers chapitres sont consacrés aux courants 
alternatifs et aux radiations électromagnétiques. 

À la fin de chaque chapitre de ce livre, écrit à l’usage des 
étudiants avancés, on trouvera un grand nombre de pro- 
blèmes d’application, certains avec leur réponse. 

En résumé, cet ouvrage permettra à ses lecteurs de 
manier avec aisance et sous leur forme mathématique 
moderne, les notions fondamentales de l’électromagnétisme. 

N. ZELLER. 


Kirrez (C.), Éléments de physique statistique, traduit de 
l’américain par A. Desprès. (1 vol. relié de 250 p., 
16 X 25 cm, Dunod, Paris, 1961, prix 32 NF.) 

Nous sommes heureux que cette traduction mette à la 
portée des étudiants français un des meilleurs livres de 
physique statistique publiés depuis une dizaine d’années. 
Il est clair et concis (Tolman est illisible en raison de sa 
prolixité), fait appel, bien entendu à la mécanique quan- 
tique, et à la notion de trace, mais il n’y a pas de surcharge 
de théorie. Des questions modernes telles que celle des tem- 


pératures négatives, sont exposées. Les divers chapitres 
sont suivis d'exercices ; ce qui rend cet ouvrage bien com- 
mode pour les professeurs. Voici les matières traitées : 
I. Principes fondamentaux de la mécanique statistique 
(entropie, équilibre, correspondance avec la thermody- 


namique, mécanique quantique, fonctions thermodyna- 
miques, distribution de Maxwell, potentiel chimique, molé- . 


cules diatomiques, aimantation, distribution de Fermi, de 
Bose, matrice densité, températures négatives). — IT. Fluc- 
tuations, bruits de fond, thermodynamique des phénomènes 
irréversibles. (Mouvement brownien, équations de Fokker- 
Planck, relation de Onsager. Principe de la production 


minimale d’entropie.) — III. Méthodes cinétiques et théorie 


du transfert. (Bilan détaillé, théorème H, noyaux, équation 
de transfert de Boltzmann, conductivité électrique et ther- 
mique d’un gaz d’électrons, magnéto-résistance, viscosité, 
etc.) Enfin quatre appendices, mathématiques et sur les 
calculateurs électroniques. 

J. WINTER. 


Jean Rossez (Professeur à la Faculté des Sciences de 
Neuchâtel, Directeur de l’Institut de Physique), Physique 
générale. Éditions du Griffon, Neuchâtel. Diffusion 
Dunod Paris. 4 vol. relié, 486 pages, 1960. 

Cet ouvrage s’adresse à un public bien informé de la 
physique élémentaire fondamentale et désireux de se fami- 
liariser avec les notions de base de la physique théorique 
moderne. 

L'auteur passe ainsi en revue les grandes divisions de la 
Physique, avec l’idée de concentrer un acquit classique 
équivalent au Certificat de Physique générale, et de laisser 
place à un exposé clair et concis des notions que le lecteur 
physicien acquerra plus en détail par la suite. 

C’est dans cet esprit que sont exposés en 510 pages la 
constitution de la matière, la mécanique, les phénomènes 
ondulatoires, la thermodynamique, et enfin l'électricité et 
l’atomistique. 

De nombreux problèmes et un mémento mathématique 
complètent très heureusement cet ouvrage, que sa densité 
d’information et sa clarté, jointe à une présentation impec- 
cable, nous font vivement recommander aux étudiants et 
aux ingénieurs de recherche. 


ZELLER. 
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CALCUL DE LA PROPORTION DE COLLISIONS INTERDITES DANS LES INTERACTIONS 
ENTRE II MESON ET PROTON LIÉ A L'ÉNERGIE DE 5,0 GeV 


Par CLAuDE BASTARD, JEAN GREA, Josepn LAFOUCRIÈRE et Georces PHILBERT 


Institut de Physique Nucléaire, Faculté des Sciences, Lyon. 


Résumé. — Nous avons calculé la proportion des collisions élastiques interdites dans les inter- 
actions t-protons liés à l’énergie de 5,0 GeV. Les hypothèses de base sont les suivantes : 

1. La section efficace de diffusion élastique x__-proton lié est égale à la section efficace x__-proton 
libre. Sont simplement exclues les collisions aboutissant à un niveau déjà occupé par un proton. 

2. La répartition énergétique des protons intranueléaires est celle d’un gaz de Fermi à la tempé- 


rature de 0 °K. 


Deux méthodes sont employées : la première utilise une variation continue des paramètres de 
choc ; l’autre est une méthode de Monte-Carlo ; elles fournissent des résultats en accord satis- 
faisant l’une avec l’autre et la proportion des collisions interdites est égale à 0,17. 


Abstract. — We have calculated the ratio of elastic collisions forbidden in interactions between 


r and bound proton at 5.0 GeV energy. 
The basic hypotheses are the following : 


1. The cross-section of elastic scattering x_-bound p is equal to cross-section tæ_-free proton. 
Collisions leading to a level already occupied by a proton are simply excluded. 

2. The energy distribution of intranuclear protons is that of Fermi gas at 0 °K temperature. 

Two methods are used. The first one utilizes a continuous variation of collision parameters ; 
the other is a Monte-Carlo method. Both give results in a satisfactory agreement and the pro- 


portion of forbidden collisions is equal to 0.17. 


Aux énergies élevées, les interactions entre les 
particules élémentaires (mésons, nucléons) et les 
noyaux sont considérées comme dues à une suite 
de collisions individuelles : la première entre la par- 
ticule incidente et un nucléon du noyau ; les sui- 
vantes entre des particules émises dans une des 
collisions et d’autres collisions. Un tel mécanisme 
peut être tenu pour valable si la longueur d’onde 
des particules en collision est faible par rapport 
aux « dimensions » du noyau. Les trajectoires 
sont décrites classiquement, les caractéristiques de 
chacune des collisions individuelles sont quan- 
tiques. Celles-ci sont dérivées de données expéri- 
mentales et, dans un certain nombre de cas, se 
retrouvent à partir de considérations théoriques. 

Cependant les collisions entre une particule et un 
nucléon lié ne sont pas rigoureusement semblables 
aux collisions entre la même particule et un nucléon 
libre. En effet : 

a) l’état du noyau intermédiaire (noyau après la 
cascade complète) doit être énergétiquement pos- 
sible ; 

b) l’état de chaque nucléon émergeant d’une 
collision doit être permis : autrement dit, cet état 
ne doit pas être occupé par un autre nucléon de 
même charge. 

L'étude de la condition a) exige un examen com- 
plet du mécanisme de cascade. Cette condition joue 
certainement un rôle négligeable dans les inter- 
actions avec un noyau moyen ou lourd ; en ce cas, 
en effet, la cascade est longue, le transfert d'énergie 
au noyau considérable et le continuum atteint 


presque toujours. La condition b) est donc la seule 
qui soit importante. 

Dans l’approximation semi-classique, le traite- 
ment habituel de celle-ci est le suivant : les proba- 
bilités absolues de collisions d’un type donné (défini 
par la nature des particules secondaires et leurs 
angles d'émission) sont tenues pour être les mêmes 
que si les nucléons étaient libres [1]. Sontsimple- 
ment exclues celles conduisant à des états interdits 
des nucléons [21]. Cette méthode est certainement 
insuffisante. Notre propos n’est pas de la discuter 
ici. Nous admettons qu’elle constitue une approche 
acceptable du problème particulier abordé dans 
cet article. 

Les nucléons dont l’ensemble forme le noyau 
sont supposés indépendants ; leur état est défini 
par leur moment linéaire. En première approxi- 
mation (modèle du gaz de Fermi à une température 
de 0 0K), tous les états compris entre l’énergie ciné- 
tique T = 0 et l’énergie maxima T — 22 MeV, 
sont occupés. Soient respectivement o» et or les 
sections efficaces totales d'interaction avec un nu- 
cléon lié déterminé et avec un nucléon libre ; 
dépend du nucléon lié en interaction ; en effet les 
paramètres de chocs sont fonction de l'impulsion de 
ce nucléon. Expérimentalement, c’est une valeur 
moyenne < 6» > qui est atteinte. Soient o$ et 0 
les sections efficaces de diffusion élastique. 

Nous nous proposons, dans le cadre de l’approche 
envisagée, de déterminer exactement le rapport 
< oÿof >. Nous utilisons deux méthodes diffé- 
rentes : 
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1. Une méthode par interpolation : Les o, sont 
calculés pour un certain nombre de paramètres de 
chocs couvrant tout leur intervalle de variation ; 
les valeurs intermédiaires sont obtenues par inter- 
polation. La précision sur < o6$/0f > est limitée 
par celle des calculs numériques et par la validité 
des interpolations. 


2. Une méthode de Monte Carlo : La seconde mé- 
thode est, du point de vue logique, aussi satisfai- 
sante que la première. Les calculs numériques, ce- 
pendant, y sont plus facilement précis. Son emploi 
permet d'autre part d'obtenir une indication sur la 
validité de la méthode de Monte Carlo appliquée à 
des phénomènes nucléaires. 

Les paramètres de chocs ont été tirés au sort ; 
les o$ sont calculés dans chaque cas individuel ; les 
valeurs s’en déduisent par une simple moyenne 
arithmétique, sans aucune pondération, tendant à 
corriger les distorsions des distributions tirées par 
rapport à des distributions homogènes. 


I. Cinématique d’une collision pion-nucléon lié. 
Paramètres définissant la collision. — Les para- 
mètres définissant analytiquement les caractéris- 
tiques d’une collision entre un pion et un nucléon 
lié se définissent le plus simplement ; pour la voie 
d’entrée dans un référentiel lié au noyau avant la 
collision (pratiquement système du laboratoire : 
LS) ; pour la voie de sortie dans un référentiel lié au 
centre de masse (CMS) du pion incident et du 
nucléon lié. Les paramètres que nous utilisons 
sont : 

— pour la voie d'entrée : les modules des impal- 
sions du pion et du nucléon en interaction et 
l'angle « de leurs directions ; 

— pour la voie de sortie : les angles € et 0 défi- 
nissant respectivement l’orientation du plan de Ia 
collision et la direction du pion diffusé. 

La distribution des € est isotrope, tandis que 0 
est donné par la courbe de section efficace diffé- 
renticlle de diffusion élastique 7—-proton libre. 

Une telle courbe (fig. 1) a été obtenue par 
Maenchen et col. [3] au moyen d’expériences por- 
tant sur un faisceau de pions ayant des énergies 
comprises entre 4,5 et 5 GeV. Des considérations 
semi-théoriques basées sur le modèle optique [4] 
permettent d’utiliser au mieux les données expéri- 
mentales. Deux observations importantes doivent 
être formulées quant à la courbe de la figure 1 et à 
l'emploi qui peut en être fait. 

a) Cette courbe ne correspond pas à la diffusion 
d’un pion d’énergie déterminée mais à celle d’un 
ensemble de pions d'énergies variant entre 4,5 et 
5 GeV (système du nucléon au repos). 

b) Dans le problème que nous traitons, les éner- 
gies dans le système propre du nucléon sont disper- 
sées entre (E + 1) GeV et (E — 1) GeV où E est 
l'énergie du pion incident (LS). À chacune. des 


énergies correspond une courbe de diffusion parti- 
culière. 

Dans une première approche nous avons négiigé 
les déplacements des courbes de diffusion autour de 
la courbe moyenne correspondant à l'énergie Æ. En 
fait les courbes varient continûment et en première 
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approximation l'élargissement dû aux variations 
de Æ vers la borne supérieure £ + 1 est compensé 


par la rétrécissement dû à ces mêmes variations vers 


la borne inférieure E — 1. 

Une approximation plus poussée tenant compte 
des variations de la courbe de diffusion avec l’éner- 
gie a donné d’ailleurs des résultats très voisins de 
ceux fournis par la méthode approchée. 

L’exposé de cette seconde méthode fera l’objet 
d’une publication ultérieure. La déformation de la 
courbe de diffusion avec l’énergie étant négligée il 
est licite d’utiliser les résultats de Maenchen sans 


- les corriger pour l’effet d’élargissement dû à la dis- 


persion des vitesses dans le faisceau. Nous avons 
donc admis que ces résultats donnaient correc- 
tement la section efficace de diffusion élastique 
pion-nucléon libre à toutes les énergiss comprises 
entre 4 et 6 GeV. 


IT. Première méthode de calcul du coefficient 
d'interdiction. — «) Considérons d’abord le cas 
particulièrement simple où l'impulsion du nucléon 
est nulle (LS). Les formules donnant l'énergie finale 
du nucléon après le choc dans ce système se ré- 
duisent à 

ER = MYy?— M{y? — 1) cos 8 


us re En 00e 
NE ER EE 


où M est la masse au repos du nucléon, Em et Pro 
l’énergie et l'impulsion du pion dans le système du 
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laboratoire avant la collision, EX l’énergie du nu- 
cléon dans le système du laboratoire après la colli- 
sion. 

La condition Ex > 961 MeV impose à l’angle de 
diffusion d’être tel que 


Myÿ®— M{y®— 1) cos 6 > 961 MeV. 


Cette inégalité définit un angle limite 6. Pour E 
(énergie du pion incident, LS) = 5 GeV, 6; — 0,139 
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radian (fig. 2). La section efficace totale de diffusion 
élastique pion-nucléon lié est alors donnée par 


Gh—È2T JE (do/dQ) sin 6 46 


Oi 
et le coefficient d’interdiction est égal à : 
‘A 81 (dod@Q) sin 6 40 
Ch 0 


AS D 
cf 1 (ds/dQ) sin 6 40 
0 


où dQ est l’élément d’angle solide dans le système 
du centre de masse. 

b) Si l'impulsion du nucléon n’est pas nulle, 
nous devons utiliser les formules générales (voir 
appendice). Elles font intervenir outre l’angie de 


diffusion 0 dans le système du centre de masse, 
trois autres paramètres à savoir l'énergie initiale du 
aucléon Ex, et les deux angles «et €. 

L'énergie du nucléon après le choc a été calculée 
pour un certain nombre de valeurs des paramètres, 
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à savoir Ex (énergie initiale du nucléon dans LS) : 
944 MeV, 949 MeV, 954 MeV, 959 MeV. 


œete: 09,309, 600, 900, 4200, 4500, 4800 


6: 40, 280 


Pour toute collision permise EX > 961 MeV. 
Cette condition se traduit pour un processus donné 
par un angle de diffusion 0 supérieur à un angle 
limite 6, fonction des paramètres e «et, Ex. Cet 
angle 6; est défini par l'intersection de la courbe 
EX = f(6) et deladroite Ex — 961 MeV. La figure 3 


80,420,1460, 200, 240, 


€" à MeV 


degres 


Ci = «h 
TGS 


correspond au cas particulier Ex, — 944 MeV, 
ei 909, 0e=300, 

Pour chaque couple de paramètres (œÆ£%,) nous 
introduisons un coefficient d'interdiction 


IF 1 (dofdQ) sin 6d0de 


2x [°° (ds/d@) sin 6 à 0 


K(xE x) 


Cette intégrale ne pouvant que s’évaluer numé- 
riquement, nous avons dû introduire un coefficient 
«partiel » X, correspondant à un angle € donné : 


15 (ds fdQ) sin 6 40 
K(e) — 


(do/d@) sin 6 40 
{ 


Nous avons ainsi construit le réseau des courbes 
K,(E) correspondant aux diverses valeurs des para- 
mètres « et Ex et Æ a été obtenu par intégration. 
La figure 4 correspond à « = 600, Ex, — 944 MeV. 

Le coefficient X correspond aimsi à une valeur 
déterminée de la direction et du module de l’impul- 
sion du nucléon par rapport au pion. Nous défi- 
nirons un coefficient global X4 en effectuant la 
sommation sur tous les nucléons, donc en faisant 
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intervenir la fonction de répartition des vecteurs- 
impulsions dans le noyau 


Ka = [ f4Exo) K(4Exu) de dExo 


et en désignant par f(«En,) cette fonction de répar- 


K4 Vo 


Pricer 
Pa VE 
20 
10 
rad 
% fai ni FL 2h sh 
6 3 2 3 6 4 
FrG.-5. 


tition. Nous avons supposé qu’elle était celle d’un 
gaz de Fermi au zéro absclu, donc : isotrope en « 
et homogène en Ex, pour Ex, < Ex, max. 

Les intégrations s'effectuent numériquement, 
d’abord sur «, puis sur Æx,. Les courbes obtenues 
dans les étapes de ce calcul sont données par les 


figures 5 (cas particulier Ex, — 954 MeV) et 6. 
Nous avons obtenu : Ka — 0,170. 


K3 À Te 
20 
10 
M eV 
à 
939 964 949 954 959 Ex 
Fic. 6. 
Précision obtenue. — L’angle 6; obtenu graphi- 


quement ne peut être défini qu'avec une certaine 


marge d'erreur. Nous avons constaté que les va-. 


leurs obtenues par répétition de la construction 
présentaient une distribution gaussienne et nous 
avons adopté comme écart maximum A6; — 20’ 
valable avec une probabilité supérieure à 99 %,. 

Les intégrations permettant de déduire le coef- 
ficient Ka ont été effectuées à l’intégrimètre. 
Compte tenu de la précision des mesures des aires 
et de la répercussion de l’écart A6; sur le résultat 
final, nous avons admis. 


RGO 1700 015; 


IIT. Emploi de la méthode de Monte Carlo. — 
1. PRINGIPE. — Un certain nombre de collisions, 
choisies au hasard, sont étudiées ; les paramètres les 
définissant sont tirés au sort, et les propriétés ciné- 
matiques des particules dans l’état final déter- 
minées ; si l’énergie du nucléon appartient au do- 
maine interdit par le principe d’exclusion, la colli- 
sion est considérée comme impossible. Du rapport 
entre les nombres des chocs exclus d’une part, des 
chocs étudiés d'autre part, se déduit immédia- 
tement le cofficient d'interdiction. 

Les distributions dans lesquelles se font les ti- 
rages de paramètres sont telles qu’il y a égalité entre 
la probabilité pour un type de collision d’êtré 
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choisie et sa probabilité physique. La précision de 
la méthode dépend essentiellement du nombre de 
cas étudiés. Dans notre travail celui-ci a été pris 
égal à 1 000. La séquence des opérations est la 
suivante : 

1. Tirage des données initiales, définissant le 
canal d’entrée ; ces données sont l’énergie Ex, du 
nucléon et l’angle «. Pratiquement, afin de per- 
mettre une généralisation ultérieure, la direction du 
nucléon par rapport au pion a été définie par les 
deux angles polaires Ëet Y (cf. appendice). 

2. Détermination de l’état cmématique des par- 
ticules dans le système de centre de masse. 


3. Introduction des paramètres du choc 0, € ; 
calcul de l’état cinématique des particules après le 
choc dans le système du centre de masse. 

4. Retour au système du laboratoire. 

La précision des calculs est telle que les conser- 
vations de l’énergie et de l’impulsion sont assurées 
avec une précision relative de 2 X 107. 


2. INTRODUCTION DES DISTRIBUTIONS. — Deux 
types de distribution apparaissent : les unes unifor- 
mes, les autres non uniformes. 

A la première catégorie appartiennent les distri- 
butions des énergies cinétiques des nucléons intra- 
nucléaires, celles des angles sphériques & et d et, 
enfin celle de l’angle € qui détermine la rotation du 
plan de choc dans le système du centre de masse. 
Dans la deuxième catégorie figure la distribution 


des angles de diffusion 0. 

Pour les distributions uniformes le tirage porte 
sur le rapport de la valeur du paramètre considéré 
à sa valeur maxima possible (soit wfu max). En 
l'absence d’une suite convenablement testée de 
nombres au hasard, nous avons utilisé un ensemble 
de 2 000 disques portant des valeurs comprises 
entre 0 et 1 et séparées par des intervalles égaux à 
0,0005. 

Quant à la distribution des angles de diffusion 
utilisée, elle est telle que la densité de probabilité 
pour un angle d’être choisie soit égale à la section 
efficace différentielle expérimentale. L’intervalle 
séparant deux classes a été pris égal à 1 degré ; la 
population étant d'environ 500. 


3. RÉSULTAT. — Le coefficient d’interdiction, 
obtenu ainsi qu’il a été dit, est égal à 0,207 + 0,03. 
Le terme 0,03 représente l’erreur probable, déter- 
minée ainsi qu'il est dit plus bas. 

Ce résultat est parfaitement compatible avec 
celui fourni par la méthode des variations conti- 
nues, à savoir 0,170 + 0,015. 


4, INTERVALLE DE SÉCURITÉ. — L’intervalle de 
sécurité exigerait l’utilisation d’une formule analy- 
tique donnant l’énergie dans le laboratoire après 
le choc 'en fonction des données initiales et des 
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paramètres de choc. A cette formule serait appliqué 
le théorème des probabilités composées. Le calcul, 
difficile pour un processus unique, serait impra- 
ticable dans le cas d’une cascade à processus succes- 
sifs. Considérant les résultats des tests auxquels 
nous avons soumis les distributions uniforme et 
non uniforme, nous avons supposé que la variance 
du coefficient d'interdiction était déterminée essen- 
tiellement par la variance partiellé relative à 
l’aléatoire 0. C’est donc de cette dernière que pro- 
vient l’incertitude sur le résultat final. L'erreur pro- 
bable qui en résulte est de 15 % sur le coefficient 
d'interdiction. 


IV. Discussion et conclusion. — La proportion 


des collisions élastiques interdites dans les inter- 


actions de r— avec les nucléons liés à l’énergie 
de 5 GeV est donc égale à 0,170 environ. Nous pou- 
vons comparer ce résultat à celui d’un calcul sim- 
pliste, dans lequel le nucléon cible est supposé 
immobile [5]. La proportion obtenue est alors 
de 0,23. La mouvement des nucléons dans le noyau 
a donc pour effet d’abaisser notablement le coefii- 
cient d'interdiction. 

Il est intéressant de chercher une indication sur 
l’évolution de ce coefficient lorsque l’énergie des 
particules incidentes s’accroît. Compte tenu de ce 
que nous venons de remarquer, la valeur calculée 
dans l’hypothèse du nucléon immobile doit être 
considérée comme une limite supérieure. Les don- 
nées expérimentales sur les sections efficaces de 
diffusion élastique manquent, elles peuvent être 
remplacées par les valeurs semi-théoriques tirées du 
modèle optique. En utilisant celles-ce1 nous avons 
trouvé un coefficient d'interdiction mdépendant de 
l'énergie des particules incidentes et égal à envi- 
ron 0,26. Aux énergies les plus élevées, il apparaît 
donc que le coefficient d'interdiction sur les colli- 
sions élastiques ne dépasse pas 0,26. 


APPENDICE 


Nous indiquerons la séquence des différents cal- 
culs. Le processus utilisé est celui qui a fourni la 
programmation la plus courte à l’ordinateur 
IBM 650 dont nous disposions. Nous avons posé 
C = 1. Le système initial du laboratoire (LS) étant 
désigné par Ox,y9Z0, nous noterons les impulsions 
respectives du nucléon et du r— par : PK, et Pr. 

Px, est repéré par ses angles polaires £ et 4 et 
Pr matérialise l’axe Ox,. On désigne par « l’angle 
(Por Pro) 

Les énergies et les masses au repos dans LS 
seront notées : 

Ex, M pour le nucléon ; 

Er & pour le pion. 

Le calcul d’un processus élastique s’effectue alors 
dans l’ordre suivant : 
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19 Passage au système Ox,7,2, ainsi défini : 
— 2, Oy, sera le plan du choc ; 
— Oz, est matérialisé par Pr. 


Fic. 7. — « angle (Px,, P) 


a 
Ÿ 


angles polaires 
repérant Px, 


Ce passage s’effectue par les rotations : 
O = (0z,, Oz,) autour de Or, ; 
D = (0x, Ox;) autour de Oz.. 


20 Passage au système Ox,y,z, lié au nucléon 
par la transformation de Lorentz sur Ox, ainsi 
définie 

Y1 = Exo]M B1 — PNolEN0. 
Les transformées sont données par : 


lPrrel = (v1/00S &) [|Prrol COS «à — 6; Errol 
Ert: = Y1lEnto — B1lPrro] COS a] 
ENo — "1" 


39 Passage au système du centre de masse : 


a) par une rotation «, autour de Oz, pour amener 


Pr, sur l'axe des x ; 


b) par la transformation de Lorentz définie par : 


82 = [Prrel/(Erre + M) 


m= (1 8% 


FiG. 8. — x,0y, plan de choc ; Ox, direction Px, ; impul- 
sion nucléon Px,; impulsion II P,,; @ rotation 
(020, 021) ; ® rotation (Oxo, Ox;). 


Les grandeurs dynamiques du système du centre 
de masse sont alors : 


— Px = Pr Ya Be M 
ER — Ye M 
EE YelErte — Bo Prrel. 


40 On se fixe les paramètres au choc : 0 angle de 
diffusion et e donnant l’orientation du plan des 
impulsions finales PK et Pr. 

On repasse alors au système du laboratoire par la 
séquence inverse. 

Manuscrit reçu le 9 juin 1961. 
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SUR QUELQUES PROPRIÉTÉS DE Cd,As, 


Par Nicore SEXER, 


Électrotechnique Générale, Faculté des Sciences, Fontenay aux Roses. 


Résumé. 


L'étude de la résistivité de l’effet Hall et du pouvoir thermoélectrique de Cd;As, 


en fonction de la concentration de porteurs, a permis de montrer que le gaz d’électrons était prati- 
quement dégénéré et qu’à la température ambiante la diffusion prépondérante était due aux 
phonons acoustiques. D’après ces résultats, la masse effective de densité d’état est : 


70/0) 08-0010) Err 


Abstract. — The study of resistivity, Hall effect and thermoelectric power of Cd,As, as a func- 
tion of the carrier concentration has shown that the electron gas is practically degenerate, and 
that electron scattering by acoustic phonons predominates:at room temperature. 

According to these results, the density-of-state effective mass is 


m* = (0.050 + 0.015) m. 


Introduetion. — L’arséniure de cadmium Cd;As 
a été signalé dans la littérature comme étant un 
semi-conducteur [1, 2] de largeur de bande inter- 
dite 0,14 eV [3, 4]. Malgré sa structure tétra- 
gonale [5]la mobilité électronique est extrêmement 
grande (de l’ordre de 10 000 cem?/Vs pour 
T = 300 K et 3- 1018 porteurs par cmÿ). 

Nous avons essayé d’obtenir des renseignements 
sur le comportement de ce composé à l’aide de 
mesures classiques de résistivité, d'effet Hall et de 
pouvoir thermoélectrique. 


Préparation du matériau (*). — La préparation 
de Cd;As, est rendue difficile par l’existence d’une 
transformation en phase solide [6] à 578 ©C, le 
point de fusion étant à 721 0C. La préparation se 
fait en plusieurs stades : 

1. Fusion sous vide des deux composés en pro- 
portion  stoechiométrique dans un tube en 
quartz [7]. Le corps obtenu est polyeristallin, très 
friable et très impur. 

2. Minéralisation dans un lent courant d’hydro- 
gène à une température inférieure à celle du point 
de transformation, il se forme des aiguilles ou la- 
melles monocristallines de très petites dimensions 
(de l’ordre du mm). Cette phase de la préparation 
constitue une purification du matériau. 

3. Fusion sous vide des petits monocristaux, 
avec dopage éventuel, pour obtenir des échantillons 
de taille convenable ; ou bien distillation sous vide 
et croissance du cristal à partir d’un germe. 

Les études ont été faites sur une série d’échan- 
tillons, en général polycristallins. Ces échantillons 


_- sont soit « purs », c’est-à-dire qu'aucune impureté 


n’y est intentionnellement ajoutée, soit dopés. 


(*) La préparation des cristaux a été faite par M. Kol- 
tirine, 


Études expérimentales. — Les mesures de résis- 
tivité et d’effet Hall ont été faites entre 78 °K et 


500 ©K sur des échantillons en forme de parallélé- 


pipède allongé, obtenu par rodage. Les contacts 
ont été soudés à l’étain. 


V2F777/} 
hermostat- 6  (i 
milhvoltmetre _. 


e ique 
a) lectroniqu 


F1G.1.— a) Mesure dela résistivité. b) Mesure de l'effet Hall 


_ millivoltméetre 


electronique 


Fil commulfateur 


a 


| ZN 19m 
D, 
Dir ; 


- [solant 
electrique 


FiG. 2. — Mesure du pouvoir thermoélectrique 
LEUT 
2— 3 — AT 
2— k AV. 


Les mesures sur un monocristal en forme de 
plaquette ont été obtenues par la méthode Philips 
[8]. Enfin le pouvoir thermoélectrique à la tem- 
pérature ambiante a été mesuré sur des échan- 
üllons contenant différentes concentrations de por- 
teurs. 
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400+- 
30 
eod- 
0 
x 
e 
= 
(et 
100;- 
IL L 1 js 1 sé 
(0) ACOM 200 300 400 500 
TK 
F1G. 3. — Résistivité en fonction de la température. 


Résultats expérimentaux. — Tous les échantil- 
lons étudiés sont de type n, les résultats exposés 
sont donc relatifs à des électrons. 


3 


L 
D Ur PUÈS n —1{Re 
Pt Be + 
2e EEE ainsi que la mobilité de Hall 
NÉ rs : ua = RÎe. 
OR EPAE TK Nous constatons, en particulier, que le dopage 
Fic. 4. — Effet Hall en fonction de la température. au germanium diminue la concentration en por- 
TABLEAU 1 
: e uQ.cm Q uV/50 R cm#/cbs er Ur cm?/Vs 
ÉCHANTILLONS 78 oK 300 oK $00 °K 300 °K n em 18 K 300 K OBSERVATIONS 
B 18 60 0,35 DOS 9 400 5 000 
C 99 110 33 0,9 A FIKULE 27 200 7 300 
M 36 120 1 5.101 27 800 8 000 petite lame 
monocristalline 
TO 162 440 68 2.1018 21 600 18 000 dopé Ge 
LV 85 255 92 1,2.1018 59 000 19 600 dopé Ge 
AT 50 60 0,07 11020 1 400 1 200 dopé In 
Ddv 84 336 ps) 2,6 2 RAALOES 31 000 7 800 « pur » distillé 
sous vide 
Bg100 92 260 65 &5 ADS 52 000 5 400 dopé Ge 
BDs 85 220 62 1,4 4,5.1018 53 000 6 400 dopé Ge 
Bg; 14 310 DE LOES 26 500 10 000 dopé Ge 


Les figures 3 et 4 représentent les courbes p(T) 


et R(T) que nous avons obtenues. 

La résistivité est une fonction linéaire de T au 
moins jusqu’à 300 ©K, ce qui est en bon accord 
avec les résultats des autres chercheurs [1, 2, 3]. 

Dans la gamme de température envisagée À est 
sensiblement constant, nous en déduisons la con- 
centration approximative de porteurs libres par la 
formule 
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teurs libres. Le germanium joue le rôle d’accepteur, 
mais il ne permet pas d’obtenir une conduction par 
trous car au-dessus d’une certaine concentration de 
germanium, il se forme un alliage. 

Les résultats obtenus sont consignés dans le 
pou { avec les caractéristiques des divers échan- 
tillons. 


Interprétation théorique. — La grande densité 
d’impuretés et la variation linéaire de la résistivité 
en fonction de la température, nous suggèrent que 
le gaz d'électrons est dégénéré. La densité de por- 
teurs, donnée par la formule générale : 


n 


& /2x m* kT\3/2 
( ] 172) 


se réduit dans ce cas à 


8 (= os on 
M = : ES 2 
SVT R? 2) 
avec m* masse effective de densité d'état, n niveau 
Ë 5 2 
de Fermi réduit = —; 
kT 


Far) = fo arfler-n + 1) de 


avec x — €/kT, énergie cinétique réduite. Par la 
suite nous poserons 


Ne = 2(2r m* AT h2)8/2. 


D'autre part, nous constatons qu’en passant d’un 
échantillon pur à un échantillon compensé, c’est-à- 
dire contenant moins de porteurs de charge mais 
beaucoup plus d’impuretés, la mobilité augmente. 
La diffusion sur les impuretés n’est donc pas pré- 
pondérante. : 

Or, nous avons : - 


Vr(£) L(E) proportionnel à {(E). 


D vitesse des électrons au voisinage du niveau 
de Fermi. 

l libre parcours moyen des électrons. 

(La mobilité est en effet proportionnelle au 
temps de relaxation + et Ür + — L.) 

Comme Vf est proportionnel à E/?et n à E??, 
n18 u est proportionnel à {(Ë). 

Dans le cas d’une diffusion par les phonons 
acoustiques, L est constant ; il en est de même du 
produit n1}$ x, ce qui se vérifie expérimentalement 
à la température ambiante (fig. 5). A la température 
de l’air liquide ce n’est plus vrai car d’autres modes 
de collisions apparaissent et deviennent prépon- 
dérants. 

En vertu de ces hypothèses — gaz d’électrons 
dégénéré, diffusion par les phonons acoustiques — 
nous pouvons déduire du pouvoir thermoélectrique 
à la température ambiante la position du niveau 
de Fermi pour différents échantillons. Nous avons : 
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Q = (ke) (x?/3n) et nous vérifions bien (fig. 6) 
que n est porportionnel à n°". 


gaz d électrons dégénérés 
diffusion par les phonons 


acoustiques 


418 19 20 


log n 


Fic. 5. — Mobilité de Hall (T — 300°K) 
en fonction de la concentration de porteurs. 
En ordonnées on porte log u. 


D’ailleurs, si nous nous rapportons aux courbes 
de la figure 7 donnant le pouvoir thermoélectrique 


log? 
-fh 
RUE 
Er 
1 ! Re 11 
18 19 log n 


Fic. 6. — Niveau de Fermi réduit 
en fonction de la concentration de porteurs. 


en fonction du niveau de Fermi réduit, nous voyons 
qu'il est justifié de considérer le gaz d’électrons 
comme dégénéré à la température ambiante dans 
l’étude du pouvoir thermoélectrique. 


Détermination de la masse effective et du libre 
parcours moyen. — L’étude du pouvoir thermo- 
électrique nous ayant fourni le niveau de Fermi 
réduit, nous pouvons atteindre la masse effective 
de densité d’état à l’aide de la relation (2), en utili- 
sant la statistique des gaz dégénérés. Ceci nous 
donne 

moin 010,00-E20;08: 
53 
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Si nous voulons un peu plus de précision, nous 
pouvons utiliser les courbes de la figure 8. Nous 


#. 


limite inférieure des”? 


limite supérieure des pouvoirs 
ther moélectriques 


LI ra-h cas dégénéré 


8 10 € 
ra 
F1G. 7. — Pouvoir thermoélectrique en fonction du niveau 


de Fermi réduit pour différentes valeurs de l’indice de 
diffusion. 


5 
Q(n = +1 É ï r) Fazer (0) 
Ë KE r) Fjjotr (M) 


k /3 Te 
ee de 


voyons icique l'erreur faite, en utilisant les formules 
approchées, est plus grande que dans l’étude du 
pouvoir thermoélectrique, nous obtenons alors la 
valeur corrigée : 


cas dégénéré. 


m*[m — 0,050 + 0,015. 


Cette valeur de la masse effective est en bon 
accord avec celle obtenue par Zdanowicz [9] 
(0,046 m à 290 °K) qui a étudié le pouvoir thermo- 
électrique en fonction de la température, mais elle 
est plus faible que celle obtenue par Turner, 
Fischler et Reese (0,1 + 0,05 m[3]). 

S1 nous supposons que la masse effective de den- 
sité d'état est égale à la masse effective réelle, ce 
qui revient à admettre des surfaces d’énergies 
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sphériques, nous obtenons pour l'expression du 
libre parcours moyen des porteurs à la température 
ambiante. ; 
1= (1 V2) EV2m*—V2ule, 
qui avec les données numérique de notre étude 
devient : 
1—1,5.10—5 cm. 


Des mesures du pouvoir thermoélectrique en 
fonction de la température sont actuellement en 


cas dégénéré Y 


ANS 
ne NX 7? 
——— ——— —1— Fit 7 
& 6 8 10 fi 
FiG. 8. — Concentration de porteurs! 


en fonction du niveau de’Fermi réduit. 


cours. Elles doivent permettre d'atteindre les méca- 
nismes de diffusion à basse température. 

C2 travail a été exécuté au Laboratoire Central 
des industries Electriques à Fontenay-aux-Roses, 
où il a été rendu possible grâce à l’aide du Mas- 
sachusets Institute of Technology. Nous tenons 
à remercier M° le Professeur Olmer qui a bien 
voulu nous accueillir dans ce Laboratoire, et 
M° le Professeur Aigrain qui a dirigé cette étude. 


Manuscrit reçu le 10 juillet 1961. 
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PROPAGATION DES ONDES DANS UN MILIEU PARTIELLEMENT IONISÉ 


Par ANNIE BAGLIN, 
Institut d’Astrophysique. 


Résumé. — On étudie la propagation des ondes de basse fréquence et leur amortissement dans 
un milieu partiellement ionisé, indéfini, uniforme, plongé dans un champ magnétique B,. On étudie 
uniquement le cas où le champ magnétique est longitudinal. 

On obtient deux modes longitudinaux et deux modes transversaux, découplés. L’un des modes 
longitudinaux correspond à l’onde acoustique, l’autre est rapidement amorti. 

Les deux modes transversaux ont des vitesses légèrement différentes encadrant la vitesse des 
ondes de Alfven. Les collisions entre les particules amortissent ces ondes et le coefficient d’amor- 
tissement est proportionnel au carré de la fréquence. 


Abstract. — The propagation and damping of low frequency waves is studied in an infinite 
uniform, partially ionised medium in a longitudinal magnetic field. 

Two longitudinal modes and two transversal uncoupled modes are obtained. One of the trans- 
versal modes corresponds to the acoustic wave, and the other is rapidly damped. 

The two transversal modes have speeds are very close to the speed of Alfven waves (obtained 


in the case of a fully ionised medium). 


These waves are damped by collision between particles and the coefficient is proportional to 


the square of the frequency. 


I. Introduction. Équation de dispersion. — Nous 
considérons un plasma partiellement ionisé plongé 
dans un champ magnétique uniforme B4. 

Nous désignons par : . la densité moyenne des 
électrons de charge ge. de masse me ; ri la densité 
mo yenne des ions positifs de charge gi de masse mi; 
ña la densité moyenne des atomes de masse ma. 

Pour étudier les ondes qui peuvent se propager 
dans le milieu nous écrivons à la manière usuelle[1]. 


1. LES ÉQUATIONS DE MAWXELL, soit en u. e. m 


dB , 
V.B= 0, (2) 
? AUDE 
VAS EnI (3) 
VRP =S LT C20. (4) 


On suppose le milieu non perturbé électri- 
quement neutre. On appelle ne + ne la densité 
perturbée des électrons, rm; + m celle des 1ons. Le 
courant est donné par : 


1 = (fe + ne) de Ve + (u + m) avi 
que l’on peut linéariser sous la forme : 
1 = le Ge Ve + Ti Qi Li 
puisque n. et nr; sont petits. La charge s’écrit : 
6 = Nee + Ni Gi. 
2. LES ÉQUATIONS DE CONTINUITÉ sous la forme 


simplifiée obtenue en négligeant les termes du 
second ordre : 


D + VA = 0 (5) 
DE + Ve De = 0, (6) 
DE + V Fa Un = 0. (7) 


3. LES ÉQUATIONS DU MOUVEMENT, linéarisées de 
la même façon 


dv 
le M — = — NY, 


dt 
+ Ne JeL(E Æ ve’ A Bo) + Peer Pie Pos: (8) 
ni Mi “ =— VY 
+ mn QE + vi À Bo) + Pie + Pii + Pin, 0 (9) 
na ra De = — VY, + Pie + Pai + Paa. (10) 


= [æ] =] 
où de, Vi, YA sont les tenseurs de pression cinétique 
de chaque type de particules, et les termes P,8 
représentent le transfert de la quantité de mouve- 
ment par les chocs des particules de type « sur les 
particules de type £. 

On fait sur ces grandeurs les hypothèses sui- 
vantes : 


— On suppose que l’on peut écrire V.#, sous 
la forme : 


VIES = Vhe + # grad div Va — Lo V° Va 


u est la viscosité interne du gaz d’espèce « et on 
suppose que les particules d’autre espèce n’ont pas 
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d'influence ; u, sera donc négligé dans tous les 
calculs HRERICINS Vp. s'écrit : 


Vha = V(nœ kTà) . 


Le passage de l'onde correspond à des variations 
adiabatiques de pression caractérisées par un coef- 
ficient de compressibilité y,,ce qui permet d’écrire : 


Va = Ya ATa Va = Ma Va Wa 
V,, étant la vitesse moyenne d’agitation thermique 
des particules d’espèce « soit : 
Vé = Ya KTa]Ma 
— les grandeurs P,., seront écrites sous la forme : 
(2) 
PB = — Na Ma VaB(Va — V8) 


ce qui définit alors les grandeurs v,, cemme fré- 
quence de collision des particules d’espèce « sur 
les particules d’espèce £. 

D’après la loi de conservation de la quantité de 
mouvement d’un système isolé on a : Pyg + Pex —0 
soit : 

Na Ma Vai — 78 MB VB 


et pour la même raison P,, = (. 


Nous appliquons ces équations au cas d’une 
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Un vecteur A a pour composantes sur les trois 
axes de référence de ce trièdre : Ax, Ay, Az. Le 
champ magnétique uniforme B, est décomposé 
en Br sur Oz et Br dans le plan perpendiculaire. 
Nous orientons le trièdre de façon que Br soit 
dirigé suivant Oy. Les grandeurs caractérisant 
l’onde sont écrites sous la forme 


he Aie NES 


w est la fréquence, k le nombre d'onde. 


On a alors ny — — Ty Vyy et VV, A pour CoMpo- 


0 
santes 
Li 

Ua À? Vax, Ua K? Va, Uo KR? Voz — FU V& Na Vas ; 

L’élimination du champ magnétique entre les | 
équations (1) et (3) donne : 

AE dj 

VV EE Rs 


= 0. (11) 

Les équations de conservation de la quantité de 
mouvement (8), (9), (10) et l’équation (11), pro- 
jetées sur les axes de coordonnées permettent de 
constituer un système de 12 équations linéaires 
dont les variables sont : 


onde plane se propageant suivant la direction de Dex Pix Paz; Per er Pin: 
l’axe Oz d’un trièdre trirectangle de référence. Dee Dis Vas PRE TE de. 
— j0 lle Me Vex + le Je Ex + le Ge(Vey BL — Vez BT) — Vea le Me(Vex — Vax) — Vei le Me(Vex— Vix) — We k? Lex = 0 (12) 
— j0 Te Me Vey + le Je Ey — le e BL Vex — Vea le Me(Vey — Vay) — Vei Re Me(Vey — Viy) — We k? Vey = 0 (13) 
— j0 De Me Vez + le Je Ez + lle e PT Vex — Vei lie Me(Vez — Viz) — Vea lle Me(Vez — Liz) — He K? Lez 
+ it Me Ve Ver — 0 TA) 

— jo Pi M Vix + Ti Qi Ex © Ni Qifviy Br — Vis Br) — Via Ni Mi(Vix — Vax) — Vie Ni Mi(Vix — Vex) — ui À? dix = 0 (15) 
— jo Pi mi Viy + Ni Qi Ey — hi qi Br Liz — Via Ni Mi(Viy — Lay) — Vie Ni M(Viy — Vey) — pu À? viy = 0 (16) 
— jo Ni Mi Viz + Ni Qi Er + ni Qi PT Vix — Vie Ni Mü(Viz — Vez) — Via Ni Mi(Viz — Var) — ui k? vi, 

—- je mi Vivi; = 0 (17) 
— j0 Pa Ma Vax + Ua À? Uax — Vai lila Ma(Vax — Vix) — Vae fa Ma(Vax — Vex) = 0 (18) 
Te Jo Na Ma Vay À Ua R2 Vay — Vai la Ma(Vay Te Viy) — Vae la Ma (Vay ei Vey) = 0 (19) 
— j0 Da Ma Vaz + Va K? Vas + _. V2 Pa Ma Vaz — Vai Ta Ma(Vaz — Liz) — Vae a Ma(Vaz — Vez) — 0 (20) 


OMTEHCE ’ à) 
_ 2 ( "2 1) Ex — Ne Je Vex — Mi qi Vix — 0, (21) 
Joke , sl d 
LT Al wo? ") 1) Ey — De e Vey — Ni qi Viy — 0, (22) 
10 à 
re: Ex — Te Qe Ver — Vi Qi Vis = 0. (23) 


Ce système est donc linéaire et homogène. Sa 
compatibilité s'exprime par l’annulation d’un dé- 


terminant du 12° ordre qui donne Péquation re- 
liant X et © c’est-à-dire l’équation de dispersion 
pour les ondes planes. 

Pour mettre en évidence dans ce déterminant les 
grandeurs physiques intéressantes nous avons di- 
visé les équations 


(12) 
(15), (16), (17) par ni mi 
(18), (19), (20) par na Ma 


) (13), (14) Par ne Me 
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ce qui fait apparaître les fréquences gyroma- 
gnétiques et les fréquences de plasma notées de 
façon évidente : 


de Br, qi Br, 
Or, = Q = — 
ke Me er mi 
RE Due LOT 
a __ AT C? lie G8 oz —4T mg 
Op re" à p NE Te NS ei à 
Me mi 


On remarque que, en l’absence de champ magné- 
tique transversal, l'équation obtenue se décompose 
et se met sous la forme d’un produit de deux fac- 
teurs qui sont deux déterminants : l’un constitué 
par les éléments communs aux 8 premières lignes 
et aux 8 premières colonnes, l’autre par les élé- 
ments communs aux 4 dernières lignes et colonnes. 

Ce dernier constitue l’équation de propagation 
des modes longitudinaux. On peut remarquer 
immédiatement que cette équation est du troisième 
degré en k?, donc qu’il existe 3 modes longitu- 
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. [RE Ve 
RUE 1) #S 
F (Vei + Vea) 
mn. M RON EE 
Vie 10 ( AE Ce) 1) 
— (Vie + Via) 
Vae Vai 
— De e — Ni Gi 


A. EN NÉGLIGEANT LES COLLISIONS (vs —= 0)on 
peut écrire cette équation sous la forme : 
— (A2 VE — w°) [(k? Va — w?) (K? V? — w?) 
+ of(k? V3 — w°) + Q(K? VE— w°)] = 0. 
On a donc un mode atomique : 4? Và — ©? — 0 
et deux modes correspondant aux oscillations des 


particules chargées. Quand «© est petit ces deux 
modes ont la même limite : 


OB(R2 V2 — «2) — QB(k? VE — w°?) — 0 
soit : - 
o? of V? + OQ Va FÉA e 
en supposant Ye = Yi, Te — Ti, et en posant 
s ; V2 = ykT Re +74 


TLe Me + ni Mi 


V est la vitesse de propagation du son dans le 
milieu moyen constitué par les particules chargées. 
L’onde est donc-une-onde acoustique, 
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dinaux. L'autre par contre, constitue l’équation de 
propagation des modes transversaux. Si on néglige 
les viscosités internes (ce que nous ferons toujours 
par la suite) cette équation est du deuxième degré 
en k?, ce qui signifie qu’il y a deux modes trans- 
VéTSaux. 

Une autre particularité se produit si le champ 
magnétique est nul. Alors l’équation des modes 
transversaux est le produit de deux facteurs égaux, 
du premier degré en Æ?. Les deux modes sont alors 
semblables, polarisés rectilignement l’un sui- 
vant Ox, l’autre suivant Oy. 

Dans ce cas du champ magnétique purement 
longitudinal il y a indépendance des modes longi- 
tudinaux et transversaux. 

Mais pour un champ magnétique quelconque 
l’équation de propagation ne se décompose plus. 
Elle est du 5° degré en k?. Les cinq modes sont 
couplés. 


II. Étude des modes longitudinaux. — L’équa- 
tion de dispersion des trois modes s’écrit : 


Vea Pc 
Me 
Via si 
mi 
10 
k2 V2 
jo | * ] 0 
A (ai SF Vae) 
Jo 
f &T c? 


B. LORSQUE LES FRÉQUENCES DE COLLISIONS NE 
SONT PAS NÉGLIGEABLES. — L’équation de dis- 
persion peut se développer : 


— (RE PE — oo?) [(K? V8 — o?) (K? V} — o°) 
+ @B(k? V?— &°?) + QB(K2 VE — w°)] 

— jo(kt Va — w!) [(wi + Va) (E? VE — !) 
+ (Vie + Via) (4? VE — ©?) + va QË + via of] 

= J0| Va Seul FE ARTE?) 
+ of(E? PV? — w?) + BA VE — w!)] 
+ ©(REVR c + V2 b, + VE as) — (a, + db + c2)] = 0 


ds V2, C2, Sont des polynômes du second degré par 
rapport aux fréquences de collision définis par : 


Ag = (Via + Vie) (Vai + Vae) — Via Vai 

= Vie(Vai +. Vae) + Via ae 
ba — (Vei + Vea) (Var + Vae) — Vea Vae 

= Vei(Vai + Vae) + Vea Vai 


Ca (vei + Vea) (Via + Vie) — Vie Vet- Eh 
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Nous étudions les solutions de cette équation 
pour les basses fréquences. Elles sont obtenues en 
annulant les deux binômes constitués par les 
termes en £? w et en k‘ d’une part, les termes en 
k? w et en w? d’autre part ; soit : 


D Lo Rae AE PA EL OA 
©? (ven O3 + via O8) VE + (vai + Vae (og V? + QZ Va) 
| (D 


k? te (2 Q2 + Via oÿ Vai + :) | (Il) 


RPC PET 


Compte tenu des relations existant entre les fré- 
quences de collision ve et vae d’une part, via et vai 
d’autre part la solution (1) peut s’écrire en sim- 
plifiant par (vai + vue) 4T c? ng°? et en supposant 


Re —= Ni soit de —= — Qi et Ye —= Yi — Ya, 


"Ra 7] 


[- 10 — (vai AE Vae) + J PS Vez + Vae Viz F Vai Vaz = 0 
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Fissnlie Le 
Se Ne + Ni + Na (1) 
2 KT neMe + Ni Mi + na Ma 
k2 . Vai + Vae le + Ni + Na / | 
CRE EEE 


On a donc ici encore une onde acoustique qui se 
propage avec la vitesse V? = yAT Ju où p est 
la masse molaire moyenne. La solution (Il) | 
s’amortit sur une longueur d’onde puisque 
RCIP) EE | 

On peut alors connaissant ces solutions calculer 
les rapports des vitesses à l’une d’entre elles. L’éli- 
mination du champ électrique permet d’écrire trois 
équations dont les inconnues sont Uez, Liz, Vaz : 


ven Dez + [fe — (vai + vs) + ji — o Dis + (va +2) Dax = 0 


o Jo 
Q2 de "LOVE Q 
Via Vez + [ve + à] Viz + [- Jo — (Vie +-Via) © — 7 ie] Vaz — 
RO) ao / 
donc : 
| © Vei — Vei © + jo 
De 2 © Via 10? + (Vie + Via) © — IR? Vi — j0ÿ 
Us PR ne — Vei © + Jo | 
— Vie © + 705 10? + (Vie + Via) © — jk? Vi— ja 


et l’on obtient de même v;/v,. Au premier ordre 
en w on a: pour (1) vw — v,, w — v, les trois 
sortes de particules vibrent en même temps avec 
la même vitesse ; le gaz se déplace en bloc ; 

pour (II) wfv, = vifvs = — nal(ne + m) les 
gaz d’atomes neutres et de particules chargées se 
déplacent en sens contraire et frottent l’un sur 
Pautre ce qui explique l’amortissement très rapide 
de cette onde. 

Pour l’onde acoustique (solution I) un dévelop- 
pement à l’ordre supérieur en w permet de calculer 
le terme d’amortissement. Il faut alors considérer 
la solution de l’équation : 


1@[(Vea O5 + Via 8) + (vai + Vae) (8 + Q3)] 
— je À { (ven 0 + va of) PA 
+ (vai + Vae) (oë VÉ + Q$ Vä)] 
— ©? [oi + Q$ + 02 + bi + a} + RVE(OS + Qÿ) 
+ of V? + O8 VI + Vice + V? ba + Ve ai] } 


Le Palo V? + 08 VE = 0 


; 2 
sous la forme . æ ms (L +) où » est mfiniment 


petit avec &. En désignant par n et x le nombre 
moyen de particules et la masse molaire moyenne 
on obtient : 


Re TOURS 
L LL Na (Vae + vai) 
ss 952 à WP 
[ra mi + Ma me —À (me + mi + ma) +#]-2 3 
ñ n u V 
avec 


WE = a{V2— V2) + ba(P?— V?) + cu(V— Pa). 


En négligeant m./mi devant 1 on aboutit à la 
forme simplifiée 


ON [at + na) À 
Vai Vae L(27i + na)? 


Vea Vae (ni + me] 

Q3 (27; + na)? 
__ 1 Vei (ni + na), 
Qp (2ni + na)? 
Le deuxième terme ne fait intervenir que la 
fréquence de collision des ions sur les électrons ; il 
s’annule quand le nombre d’ions s’annule. Il repré- 
sente l’amortissement dû au frottement du gaz de 
particules chargées sur lui-même. Le premier terme 

peut s’écrire de façon simplifiée : 


n = — 
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en remarquant que vue —+ vai est proportionnel 
à mr; et que vi est portionnelle à ra, À est homogène 
à un temps. 

ni 
Ni + la 
sation. Le coeflicient d'amortissement s'exprime en 
fonction de X par 


Désignons par X — 


2 1— X 
Roc dar 
Ce deuxième terme représente la friction du gaz 
- de particules neutres sur le gaz de particules 
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le degré d’ioni-- 
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chargées. Il s’annule lorsque X — 1 soit rm — 0 
et croît lorsque r, augmente. 3,5 


III. Étude des modes transversaux. — Nous 
étudions ici uniquement les ondes de basse fré- 
quence ; 1l est alors légitime de négliger les cou- 
rants de déplacement ce qui, dans les équations de 

: À : 1 dE 
départ, revient à négliger le terme en Pen. et 
permet alors d’écrire l’équation de dispersion des 
deux modes transversaux sous la forme légèrement 
simplifiée suivante, en annulant le déterminant : 


JO Vel — Ven Vei Vea _ Tr, OL 0 0 0 

Me 
, , gi Q 

Vie TSI A VEt Via RE 0 ° SL 0 0 
mi 

Vae Vai Cry" Jo— Vai —— Vae 0 0 0 0 0 

— krcne qé © —4nrc?nqto 0 jk?c? 0 0 0 0 

OT, 0 0 O — jo — Vei — Vea Vei rs Ze 

Me 

c qi 

0 — Qr, 0 0 Vie RO AVie EL Via Via = 

Mi 

0 0 0 0 Vae Vai — jo EL ee 7 0 

0 0 0 0 —4rc?n qèo —4nc? ni q?o 0 jkèc? 


On posera pour simplifier l’écriture 
RC CH. 


On obtient une équation du second degré en w 


et du sixième degré en w. On se propose de recher- 


: = k Me 
io i re 
] ei e de 
Vie — JO — Vie _. 
Ju 

7 Île Je O TT ugio AR: 
OT, 0 0 
0 — Qr 0 
0 0 0 


{Pour les très basses fréquences les solutions sont 
celles de l’équation obtenue en conservant unique- 
ment les termes en 2, uw? et «4 soit (26) 


+ of, Of, u? + 2u w°{of, Q8 + Q$ wÿ) 


— oUuÿ + Nÿ)? = 0 


mtO Era Vei 


TT le Je © 


cher la forme de la solution u(w) au voisinage de 
ou 0. 


A. MILIEU COMPLÉTEMENT IONISÉ. — L’équa- 
tion (24) se réduit à (25) 


Vie a OV I _— 


mn EC) 


qui compte de 
(of, Q$ + QL wp)?— (oi + Qp)°? of QŸ = 0 


s’écrit en divisant par — wf : 


2 


Peer 
&° NI GEMEE 2 
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On vérifie alors que les deux ondes identiques 
obtenues sont les ondes de Alfven. Leur vitesse de 
propagation .est 

Br: 
ne. 


o7, Qx 
(of, Q3 + of, op) /? 


LA (3) 


où p désigne la densité moyenne du gaz 
p = Te Me + li Mi. 


En fait ces ondes sont amorties par les collisions. 
Pour calculer le coefficient de cet amortissement 
il faut développer l’équation (25) à l’ordre 6 en ©, 
faisant intervenir les termes en u? wo, uw, et 
A io 


En divisant par — ww? QZ? on obtient (27) : 


ii ge) + 2 Eat Sue (ue _ à) 
DFA r) of, QË GNT 
__ 4 (Vei + Vie)? + où + Of 0 
w? op QÀ D 
Posons 
QË Vei + @É Vie Ya (Vei + Vie)? + Of Le Op Sp: 
Ë QÉ F of NÉ 


Cherchons la solution de l’équation (26) jusqu’au 
3° ordre en w, sous la forme u/&? = (1 +xw)/V&, 
x est alors solution de l’équation du second degré 
(28) 

x? + 2J4Ax— B = 0 


qui s’écrit : 


x = — jA 4 VB — A? 

mais 

(vei + vie)? + QÉ + OT, 

LES 
B — À AC 
__ (QŸ vei + ©f vie)? _ Q + of 
où OË horun 

est toujours positif, donc : 

x —=— JA + va? + oflor Qr. 


Puisque k? c?/w? — (1 + xw)/VX? on obtient 
deux ondes de vitesse légèrement différentes, amor- 
ties par les collisions. 
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La présence de particules neutres lève donc la 
générescence. Soit V la vitesse de ces ondes ;on a. 


2 OZ L n2\1/2 
Pa = (s : NA + of \ 
4 SD ME 
soit 
Va _. 1 VOi 2 ©, 
Ha 2 or OL . 
ou encore 
For 
Fe 207 


Le coefficient d'amortissement est 


No 
Mc Va 


LA 


donc 


,  (QÈ vei + of vie) (0$ of + œÿ Qf)7? 
PA — 
2c of, Qÿ 


qui compte tenu du fait que 


s’écrit 


SEA Met reVoie + Vei NET 
2c(Qz or) 7? &L.) 103 ® 
Domaine de validité. 


2 2 
QT Vei + Of, Vie 


<a | 
of, + Of 
soit 
Q? 
PR de Me 
Vie Vei 


et compte tenu de Me Vei — Mi Vi, le terme pré- 
pondérant est le 2° et la condition est © << &£ /ve. 
Pour nr; — 1023 il donne 


o € 108 Bf,. 


B. CAS GÉNÉRAL. — Le nombre d’atomes neutres 
n’est pas négligeable. 

On développe l’équation (24) en utilisant pi 
exemple la règle de Laplace [4] avec des mineurs 
d'ordre 4. Ce qui permet d’écrire : 


(24) = As D, — A Do + As O3 + Ag D — À D + A6 Ds = 0 


A, est le mineur-constitué par les éléments communs aux colonnes ASP 


ne 4 4 4 
"À; 4 n 4 
y 4 4 4 
AE 4 4 a 
à, 4 n 4 


4 


4 


4 


4 4 AE A 
: 4 4,5, 4:6 
a v : DS 
4 4 2,3,4,6 
° 
EN 4 


3, 4, 5, 6 
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et aux lignes 1, 2, 3, 4. 

5, est le mineur constitué par les éléments com- 
muns aux 4 lignes et aux 4 colonnes qui ne figurent 
Had dent Aer (le 2845 6) 

On remarque que 

A= 0 Ad Ad A5= 0 A=d As = 0 
et l’on obtient 
A — 21; 09 + A + AË + A = 0 


avec 
L] A? = &w°[— u? P? + 2jo uPQ + w°Q°? 
P = 0 © + Bio + y Q = a © + Ê2 
% = ] 
Bi = Si 
n=— ii +) 


Ma 
a = — (Q$ + of) = — wi — 
mi 


| 
HET te 7 El RE 
d=a+h=c 
S1 — Vai + Vae + Vei # Vea + Vie + Via 
e As 02 = u* io) + ju © q'(o) + w2 gs) 
qe) = o1 Qr[o? vei Vie — J@(Vie ba + Vei 42) — @2 bo] 
go) = oz O1 { o°(o8 vie + O8 vei) 
— jo[(vai + ae) (O8 Vie + Q$ vei) + © 42 + Q$ b:] 
— (vai + ae) (oë 4 + À 8) } 
go) = o1 O1 Q$ œi[jo + (vai + Vae)]? 
# AË = url) + Ju wo rslo) + ro) w?. 
On pose 


Gi = Vei + Vea + Vai + Vae 
bi = Vie + Via + Vai + Vae 
ral) = — Qf[owt — 2j0° a; — w°[ai + 2b;] 
+ 270 a1 ba + bo] 
olje® + w°%{a, + vai + Ve) 
— jo[(Vai + Vae) 41 + bol + (vai + Ve) 2] 
ro) = Q$ wf[jo® + (vai + vae)]? 


@ A3 = u? salo) + jo us3(w) + saw) w°. 


rlo) = 208 


On passe de À, à À, en changeant i en e et inver- 
sement 


© A3 = — u? Q? of[jJo + (vai + Vae)]?. 


La première approximation des solutions w(w) au 
voisinage de © — 0 est donnée par l’équation (29) 
obtenue en égalant à O0 la somme des termes 
en u?, uw? et w* soit (29) 

AU ©, Qf(vai + Vae)? + 2u ©°(vai + Vae) 


[of Q$ + QÉ wo) (vai + Vae) + QE O ea + @p Où via 


2 
are Gt [nes Q$: + Via op + rh ŒB( Var + el] -= 0 - (29) 


PROPAGATION DES ONDES 817 


et en divisant par — (vai + va)? on obtient 


u? of, QË — 2u ot[of, O8 + OZ of] [1 + we | 
1 


qui s'écrit : 


2 
Na 
u_!tù) Pa 
wo? La EE 


_ On obtient donc deux ondes identiques dont la 
vitesse de propagation est V définie par 


SO1b 
Pr, Br, 


14 = = . 
Na VE&rp 
anne me + mm) (1 +2) 


e désignant la densité du gaz. 

V est une fonction décroissante de "an et 
tend vers 0 quand ra /n; -> co. 

Pour obtenir le coefficient d'amortissement de 


F 


4 
A 


Ta/7i 
HTC 

ces ondes c’est-à-dire le développement de u(«) 
solution de (25) jusqu’à l’ordre 3 en 6, il faut tenir 
compte des termes en u? &w, uw?, &°, u? &?, uw“ 
et of. 

On obtient alors l’équation (31) en ayant divisé 
l’ensemble par — © (va + ve)? 4 


u 4 \2 ; ur < u 
Et 9 SN) DES 
(£ vi) + 274: RE JB: = 

+ 2jC1 © + As 5 + Biu + Ci? = 0 


après avoir posé : (31) 


on Qy, + a, (1 +7) 
1 


AR 
À OL Qr,( vai < Vae) 
je (2 7) d (1 me) 
B Ga mi er “ 2m ne 
s Qx, &r(Vai + ae) 
2 Mà na | 
Ci —= 1 


Of, OFTET + Yae 
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2 
d3 (1 +2) + QŸ ©? — 201 &L Vei Vie + QE(vei + Vea + Vai + Vae)? 


ni 
Az fe of, Q? (vai > Vae) ? 
se Of(Vie + Via + Vai + Vae)? + 2(Q$ D: + «wf, a) 
of, QZ (Vai SF Vae) ? 
— 27 og (1 +7) [vi + ve) 51 + de] — 2(of QE + QE of) 
B; == ) à ni 


&f, Q? (Vai ce Vae) ? c 


FR EE LRU PE 
\ @f, QE(vai + ae)? 


avec On recherche le développement, limité à l’ordre 3, 
APR RE 2 es de la Se te l’équation (32). Pour cela on 


PP et et on détermine x en portant 


mais A, — B, V? +C, V4 —=0 donc (31) peut 


encore s’ re sous la ce (32) cette valeur dans (32) et en annulant le terme de 


plus bas degré en o, c’est-à-dire le terme en &? ; 


u NE : u il u . : 
(y) +47) (ar €) à | 
se x2 + 2jA,(1 — KV?)x— A, B, V?— C, Vi — (33) 
= — — — 2 — « Q 
A d’où l’on tire la forme (34) : 
avec Ro ANNE) 
2 Ma A 2 APE 5 2\211/2 
= Ge Dm ETAPE CR A,(1 KV?) RE 
ARE mn) L'expression entre crochet dans (34), compte- 
oz QL + d (: à tenu de C’, — K? A peut se mettre sous la forme : 
Q$(vai + Vase + Vei + Vea)? + @p(Vai + Vae + Vie + Via)? — 20, Qr[(vai + Vae) 51 + (Vei + Vea) (Vie + Via)] 
of, (E(Vai + Vae)? 
soit en remarquant que ce qui correspond à deux ondes de vitesse V’ définie 
a 
(Vai + Vae + Vei + Vea) (Vai + Vae + Vie + Via) pet 
= (vat Vae) S1 + (Vei + Vea) (Vie + via) A 2 (1 F x, ©) 
[Q x (vai + Ve + Vei + Vea) HE @T.( Vai Vie rie Via) ]? 
of, Qf(vai + ae)? soit 
elle est donc constamment positive. x est donc un V'=Y (1 Fe “ <) 
nombre complexe de partie imaginaire : 2 
et l’on a 
OT, O1 + de (1 + 2)" 
CRC | 
x = — Al — KV?) — : ARE 
J X1 © 
&r, Qr (vai + Vae) (1 +2) wo? y? 2 1 ©) 
1 . 
Ce SOIL 
et de partie réelle + x, avec 
n rss f(iite ET) 
(vai + Vae) (Or — or) + _ (or, + Qr) (Vae — Vai) | cV - 
1 . 
ne or, QL(Vai + Vae) ou encore 
Les solutions de l’équation (32) s’écrivent donc en tn. ®. je o) 
1 CV?  92cV 
Ke 1 : Ÿ " . . . . a 
© PANNES k a donc une partie imaginaire ; l’onde.ést amortie 


2" 0e 0 1b lu — 0) 9u — 
(P%4 + tm) ASE | 
0 o Tea 984 
Er — (: Es ] ol 
#4 ll 
Et + — 
Tu “7 cu 
7% | es (r De ] ol 
8A 4! 
(so + 14) — 
eu 
% E2A F4 MES TN 
al RE 
eA 21 
0 0 0 0 (1 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 . 0 0 
0 0 LU — 0 
0 0 0 Lo 
ZT ZeQ Z1Q 200 


0 0 0 0 
0 0 0 0 
0 0 0 0 
}o! 0 0 0 0 
2 0 15 tu — 95 ®u — 
& 44) OÙ 
(9% a Te) AT 9 
0 yen — of — TEA LAN 
pu a (Pa + a) — _ 
1h 2) inbess of En 
2 LEP LCA (Pa + cg) Éuès 
°b 4 À — of — 
0 0 0 0 (r ie 
0 0 8 0 
0 0 Lie) 0 
0 0 0 FLO 
ÂT Ava Âla 49Q 


SHTAVIUVA 


z0 


39 z 


| 


0 0 0 (83) 
0 0 0 (03) 
0 LU + 0 (LT) 
0 0 IO — (31) 
0 0 0 (32) 
0 0 0 (61) 
0 Ty — 0 (91) 
0 0 To (81) 
to ou vu | td 
(84 + I) — 
at A — a, ASE Re (81) 
(QE a) pce 
VA FU re °1a (Gt) 
(84 + Fa) — 
a 19 
KEY 4 AE, — (TI) 
re Fa *°a SNOILV AD 
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et le coefficient d'amortissement est donné par : 


na na à 
a [er O1 + à (1 nl ] cr 
For ei 1/2 
#7 o(on On) (vai + we) (1 +2) 
ï à 


formule légèrement différente de celle donnée par 
Piddington (5). Elle redonne pour à — 0 l’expres- 
sion obtenue dans lecas du milieu complètement 
ionisé : 
Vel @Op 
2c| OL Q1)#2° 
Compte tenu de ce que 


Vae D Vai Vei > Vie Ma © Mi 
on écrit 
na 
d — Yae (ve — vi) 
Ni 


et l’on peut mettre le coefficient sous la forme: 


ne 1/2 s 
20 Vre (1 + 2) (or, Qr,)1/? 
il 


m Ta \ CIE 
(ru +2 vx) 0 (1 +2) 0! 
L Me Ni 


2c{or, Qr)7E 


Le terme prépondérant est en général le deu- 
.xième. Il s'exprime en fonction de la vitesse de 
propagation : 


y = Br 
Fr 1 
Vmin (1 né 2) 
Ni 
par 


Vr w°? 
ri q° V3 


(Mi Vei + Ma Via) 
À ; 27 

ou en fonction de la longueur d’onde À = V par 
(a) 


87° mi Vr CS 
PET (Vie + Via) 53 


= 
rs 


SMS. | 
= ma 
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Evaluation numérique de ce coefficient dans le cas. 
de la photosphère solaire. 


Na — 1017 
M —= TO 
TU ADE APR 


Vao — Te G Doi 2,4 40? 
ver = 25716 Te VS 9,4 — 2,540 


Le premier terme donne 10% B51 w?. 

Le deuxième 5.10—2.B5* w? en négligeant ve 
devant va. Il est en effet le terme prépondérant, 
tant que Br, < 3 105 Gauss. 

Le coefficient d'amortissement qui est donc 
essentiellement dû aux frictions du gaz de parti- 
cules chargées sur le gaz de particules neutres 
s'exprime par : 

5 10—2 Br? w°. 

Domaine de validité du développement. 


La condition est À, w 1 
d’où 


o1 Or + de (1 +7) 
16 CA 


OT Qr (vai + Yae) 


mais Mei >> Vie et Vae > Vai SOIT 
d> — Yae(Vei + Via) 


on exprime la condition par 


Moi en G A] 


Vae ox (x, 


Dans le cas de la photosphère vi € va 
Vae = 2,4 105 


s J Via —= 4011 
il vient : 

1 1 
0—? BÈ F9 10 


1 
=> + 
a) 1 
le premier terme est prépondérant et la condition 
est 
o € 10—? B}. 


Manuscrit reçu le 25 avril 1961. 
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UNE DÉRIVATION EXACTE DE L'ÉQUATION INDÉPENDANTE DU TEMPS DE SCHRODINGER 
DE LA THÉORIE DE LA STABILITÉ DES MOUVEMENTS DANS LA MÉCANIQUE CLASSIQUE 


Par Nixoza ST. KALITZIN, 
Institut de Physique de l’Académie Bulgare des Sciences. 


Résumé. — À partir de l’équation ondulatoire de Tchétaéff, donnée dans son article [1], qui 
détermine les mouvements stables des systèmes mécanique conservatifs, nous déduisons, dans le 
cas stationnaire, l’équation indépendante de Schrôdinger pour n points matériels. 


Astract. — From the Tchétaéff waves equation [1], which gives the stable motion of conser- 
vative mechanical systems, we deduce, in the stationary case, the Schrôdinger equation, inde- 


pendent of time, for z material points. 


Dans cet article nous déduirons rigoureusement 
l’équation indépendante du temps de Schrüdinger 
pour 7 points matériels de la théorie de la stabilité 
des systèmes mécaniques dans la mécanique clas- 
sique, laquelle théorie a été crée par Poincaré, 
Liapounoff et Tchétaéff. Nous nous servirons essen- 
tiellement de l’article de Tchétaéff [1] où il donne 
l'équation ondulatoire fondamentale du mouve- 
ment stable. Selon Hmélévski [2] nous avons 
l’image approximative suivante pour le mouvement 
ondulatoire autour des trajectoires stables 
Tchétaéff a montré qu’au voisinage du mouvement 
stable d’un système mécanique conservatif, les 
petits mouvements perturbés oscillent autour des 
trajectoires stables comme un pendule, qui, déplacé 
de sa position d’équilibre, oscille autour de sa posi- 
tion stable. 

Sur chaque corps, en plus des forces connues 
introduites dans les équations du mouvement, agis- 
sent aussi de petites forces perturbatrices qui nous 
sont inconnues. Suivant le principe de la stabilité, 
ces forces détruisent chaque mouvement instable, 
ne conservant que les mouvements stables. Tché- 
taéff a trouvé l’équation à laquelle doivent être 
soumis les mouvements stables des systèmes conser- 
vatifs. Cette équation est identique à l’équation 
ondulatoire. 

Dans la nature, le mouvement ne s’effeciue 
jamais sur une trajectoire stable. A cause des 
petites perturbations il existe toujours des dévia- 

* tions insignifiantes et le mouvement s’effectue dans 
une région étroite enveloppant la trajectoire stable. 
C’est pourquoi les trajectoires réelles doivent oscil- 
ler autour d’elle (on obtient ainsi une onde parti- 
culière.) 

L'idée directrice de cet article sera que les phéno- 

_ mènes quantiques ne sont que le résultat de la 
stabilité des phénomènes mécaniques dans des sys- 
-tèmes de dimensions suffisamment petites. Cette 

_ conceptiôn s'appuie aussi sur le principe suivant de 


Tchétaéff : tous les phénomènes qui se réalisent 
effectivement dans la nature sont stable. 

Considérons selon Tchétaéff [1]un système méca- 
nique, qui consiste en z points matériels, qui 
peuvent se mouvoir librement dans l’espace. Dési- 
gnons par $ — 3n le nombre des degrés de liberté 
dusystèmepard Jess 1 2 SC 0 EURE 
Lay Yns Zn) Ses coordonnées indépendantes. Soient 
(ti, Yi, Z:) les coordonnées cartésiennes du 1-ème 
point (pi, pi, pi) les composantes de sa quantité de 
mouvement, m, — sa masse. Nous introduisons 
aussi les composantes de la quantité de mouvement 
Pi = (p}, PE, ..., p?) correspondant aux degrés de 
liberté. Nous supposons que les forces agissant 
dans le système ont une fonction de force qui ne 
dépend pas du temps 


LE = AD Te) 


Par 


27 = À À [(pé}e + (pi) + pal 
1 7 
nous désignerons la force vive double du système 
matériel considéré. 
L’équation d’Hamilton aux dérivées partielles 
prend la forme : 


Bale) +6) +6)]=26-0 « 


où Æ désigne la constante des forces vives, V est 
une fonction de q;, 4», ... gs indépendante de é. 

L'intégrale complète de l'équation d’Hamilton (1) 
est donnée par la fonction 


2 OS) CONS. 


iBSurere 
qui satisfait à l’équation (1) et dépend des cons- 
tantes æ&;, ... as. Dans ce cas la constante des 
forces vives est une fonction quelconque des cons- 


tantes «; ; 
HER EC as) 
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Suivant le théorème d’Hamilton-Jacobi la solu- 
tion générale des équations du mouvement est 
donnée par les formules È 


où les £; sont des constantes quelconques. 

Les mouvements perturbés du système méca- 
nique se déterminent par les différentes valeurs des 
constantes «; et B;. 

Pour distinguer, entre les mouvements possibles 
du système mécanique, les mouvements qui sont 
stables par rapport aux variables en cas de pertur- 
bation seulement dans les données initiales, on con- 
sidère les équations différentielles aux variations de 
Poincaré 


dé, d’H d? H 
Pass È k = ét où 
MS G& dP; EN dPj dPr n): ie 
dr Ne Pol | 
LR De À 3 
dt z É Da op, dr ? 3) 


où &,, 7 sont les variations des coordonnées g, et 
des composantes de la quantité de mouvement p, 
NTM RUES 

Pour le mouvement stable non perturbé les équa- 
tions aux variations de Poincaré représentent un 
système d'équations différentielles linéaires qui se 
réduit par une transformation non singulière liné- 
aire des variables en un système d’équations diffé- 
rentielles linéaires avec des coefficients constants, 
où tous les nombres caractéristiques (nombres de 
Liapounoff) du système des solutions indépen- 
dantes doivent être égaux à zéro. (Théorème de 
Poincaré, Liapounoff, Tchétaéff.) 

Supposons que le mouvement perturbé n’est pro- 
duit que par la variation des constantes B,, ... By, 
en fixant les valeurs des constantes &,, ... as. La 
stabilité du mouvement non perturbé dans cette 
hypothèse sera aussi une condition nécessaire pour 
la stabilité dans la variation de toutes les cons- 
tantes «; et G,. 

Les variations des coordonnées et des compo- 
santes de la quantité de mouvement dans cette 
hypothèse auront des nombres caractéristiques zéro 
si le mouvement non perturbé est stable. 

Avec l’aide des équations différentielles de Poin- 
caré (3) et de la détermination des nombres caracté- 
ristiques, Tchétaéff [1] arrive à l’équation suivante 
pour le mouvement stable non perturbé : 


ARE Es 12 (ä) 


2 
dx? 


LUATRCIA RES 

L’équation (4) est une équation différentielle 
elliptique. Nous introduisons une fonction double- 
ment différentiable ®(— Et + V) dépendante de 
l'intégrale complète — Et + V de l’équation aux 
dérivées partielles d’Hamilton-Jacobi. 
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Pour le mouvement stable non perturbé, à partir 
de (4), (1) et (2) nous avons : 


é=1 Mu \ dti dyi dz; 
mA» PF Cor For 
— EE =: TE es 
2 mi ÜTe dy? dz | 


LA TON NE dV\? dV\? 

DANN EE N EE Que RÉ 
& À mi (2 + ee 24 | 
_%E—%4)2© 
E? dt? 


et par conséquence 


22 D 
dy} 


LT 5 


sh pv 


2(F—A)22®  %, 1 [DD 
V9 dt? pi) dx? 
L’équation ondulatoire (5) donnée par Tchétaéff 
dans son article [1] établit l’analogie entre la théorie 
mathématique de la lumière de Cauchy et les mou- 
vements stables des systèmes conservatifs. 
Puisque nous considérons des phénomènes sta- 
tionnaires, nous pouvons poser 


f 


Die a (6) 
« est une constante indépendante de &,, ... œs, 
B:,... Bs. Alors nous aurons 

# 1 /D2D 220 220 S 

Le Ce 5e +5) = 2 — ) 9,0) 


Ainsi « ne dépend pas de l’espèce des pertur- 
bations du mouvement stable, ni des paramètres du 
mouvement, et nous pouvons admettre que c’est 
une constante générale pour tous les mouvements 
et toutes les perturbations. 

Nous posons 


CR 10 (8) 


où À est la constante universelle de Planck. Alors 
l’équation (7) obtient l’aspect 


SRE 8T? 


< il = 
a : De <= (E — À) ® — 0, (9) 
où 
d2D DD 22 
AN: D — | 
dr? ou d2 


est l’opérateur de Laplace. L’équation (9) est 
l'équation indépendante du temps de Schrüdinger 
pour un système de » particules. 

À ma connaissance, ce qui précède donne la pre- 
mière dérivation exacte de l’équation indépendante 
du temps de Schrôdinger des principes de la méca- 
nique classique, sans l’usage de certains paramètres 
cachés, ou de champs cachés, ce qui la différencie 
essentiellement de la théorie de Bohm. 

Cette dérivation donne une nouvelle signifi- 


l 
| 
I 
| 
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cation de la mécanique quantique, tout à fait dans 
l'esprit des idées de Louis de Broglie [3]. 

L'auteur est persuadé que la mécanique quan- 
tique pourra être dérivée à partir de la considéra- 
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tion des mouvements stables dans le cadre de la 


mécanique classique. 


Manuscrit reçu le 2 juin 1961. 
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REVUE DES LIVRES 


Dicke (R. H.) et Wirrke (J. P.), Introduction à la méca- 
nique quantique. (Introduction to Quantum Mechanics.) 

4 vol. relié de 369 p., 15 X 23 cm, Addison-Wexley, 

Reading, Mass., U. $. A. et Londres, 1960, prix $ 8,75 

(en langue anglaise). 

Un bon cours de mécanique ondulatoire à une fonction 
d’onde avec spin (mais non relativiste). Comprend un 
chapitre sur les statistiques quantiques, avec définitions 
de la matrice densité. 

J. WINTER. 


Relations de dispersion et particules élémentaires. (École 
d’été de physique théorique des Houches.) 4 vol. relié 
de 672 p.,18 X 25 cm, Paris, Herman, 1960, en français 
et en anglais, prix 84 NF. 

Cet ensemble de cours concerne les bases de la théorie 
quantique des champs et celle des particules élémentaires. 
Le point de vue est le plus élevé : on cherche à fonder la 
théorie des champs sur des bases solides et pour cela on 
fait appel à la théorie des distributions, due à M. Schwartz, 
qui généralise la notion de fonction. Le lecteur devra avoir 
une excellente culture en analyse, et qui plus est une 
culture moderne. Il devra aussi, pour apprécier les chapitres 
de la fin, connaître les faits expérimentaux de la physique 
moderne concernant les particules élémentaires. Les auteurs 
de ces cours sont des spécialistes réputés. Voici un sommaire 
des chapitres. 

I. Gozp8erGEer (M. L.), Introduction et application des 
relations de dispersion (cas non-relativiste et relativiste, 
diffusion pion-nucléon, etc.) (anglais) — IT. WicaTMAN 
(A. $.), L’invariance dans la mécaniqué quantique relati- 
viste (français). — III. WicHrmAN (A. $.), Fonctions ana- 
lytique de plusieurs variables complexes (anglais). — 
IV. Omxës (R.), Démonstration des relations de dispersion 
(français). — V. KALLEN (A. O. G.), Propriétés des valeurs 
propres dans le vide des opérateurs de champ (anglais). — 
VI. CHew (G. F.), Relations de double dispersion et uni- 
tarité comme base d’une théorie dynamique des interactions 
fortes (représentation de Mandelstam, divers types de dis- 
persion, structure électromagnétique du nucléon) (anglais). 
— VII. TreimAn ($. B.), Les interactions faibles (neutrino, 
radiation bêta, capture du méson mu, etc...) (anglais). 
VIII Yamacucui (Y.), Interactions fortes des particules 
étranges (isospin, étrangeté, isosymétrie, parités de spin, 
collision et diffusion K — N, K — N, etc... (anglais). 

J. WINTER. 


Demers (Pierre) (Professeur à la Faculté des Sciences de 
l’Université de Montréal), lonographie. Les émulsions 
nucléaires. (1 vol., 834 pages relié toile bleue, Presses 
Universitaires de Montréal.) 

Bien connu comme pionnier de la photographie corpus- 
culaire, l’auteur présente une vue d’ensemble des problèmes 


concernant les émulsions destinées à la physique nucléaire. 

Les plaques photographiques spécialement préparées, ou 
émulsions ionographiques, sont l’un des plus puissants des 
moyers d'analyse des phénomènes corpusculaires. 

L'auteur étudie d’abord en détail dans les deux premières 
parties le processus photographique, la préparation des 
émulsions, et leurs propriétés. La troisième partie de l’ou- 
vrage est consacrée à un exposé complet des conditions 
d'emploi de ces émulsions et de la technique des mesures 
(dénombrement, mesures de parcours, diffusion angu- 
laire, etc). 

Dans les trois dernières parties, l’auteur traite des appli- 
cations de cette technique aux mesures de radioactivité, 
aux réactions de faible ou de grande énergie, et à l’étude 
des rayons cosmiques. 

Cet ouvrage peut être considéré comme un traité de pho- 
tographie corpusculaire. Il faut particulièrement louer 
l’auteur d’avoir joint aux aspects théoriques des questions 
traitées de nombreux détails expérimentaux et pratiques. 

La présentation agréable de l’ouvrage, abondamment 
illustré, en rend la lecture facile. 

ZELLER. 


Synce J. L.), Relativité. Théorie générale. (1 vol. relié, 
16 X 23 cm, de 505 p., North-Holland publishing Ce, 
Amsterdam, 1960. En langue anglaise), prix 55 florins. 
Ce livre est celui d’un géomètre qui se réclame princi- 

palement de l'inspiration de Minskowski. La notion de 

principe d'équivalence lui paraît obscure et mal fondés, et 
doit être écartée, même sielle a joué un rôle historique : L’exis- 
tence d’un champ de gravitation est une propriété absolue. 

Le premier chapitre est une introduction générale de calcûl 

tensoriel ; à côté du déplacement parallèle, on introduit le 

déplacement de Fermi, qui constitue la généralisation cor- 
recte de la notion-newtonienne de système d’axes non- 
tournant. Au chap. Il, on introduit la fonction uni- 
verselle Q, égale (à + 1/2 près) au carré de la distance géo- 
désique de deux points d’univers. Au chap. III, on traite la 
chronométrie dans l’espace de Riemann. Chap. IV, Le 
continuum matériel. Chap. V, Quelques propriétés des 
champs d’Einstein (notamment les ondes de choc gravi- 
fiques). Chap. VI, Lois de conservation intégrales, et équa- 
tions du mouvement. Chap. VIT, Champs de symétrie sphé- 
riques (champ de Schwarzschild, orbites et rayons du 

champ solaire, déplacement spectroscopique). Chap. VIII, 

Quelques univers spéciaux. Chap. IX, Ondes de gravitation. 

Chap. X, Électromagnétisme (notamment univers élec- 

trovac, avec champ électromagnétique et sans matière). 

Chap. XI, Optique géométrique. Tout l’ouvrage est clair, 

illustré de nombreuses figures, et fait appel à l’esprit géo- 

métrique. Il n’est pas (ce qui mérite d’être dit) alourdi de 
nombreux développements de théorie des ensembles. 
J. WINTER. 
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TAUX DE BRANCHEMENT ENTRE LES DIVERS MODES DE DÉSINTÉGRATION DU MÉSON MU 


Par Me Arix JACQUEMIN et M. Pauz KESSLER, 
Laboratoire de Physique Atomique du Collège de France. 


Résumé. — Nous avons calculé à l’aide de la méthode des processus quasi réels, les rapports de 
branchement entre les trois modes de désintégration observés du méson mu : 


WE NEC En 


et ANA nf AE RTE à 

u—e+v+v+er +er, 
ces rapports de branchement étant fonctions de l’énergie minimum du photon émis ou de la paire 
créée. Nos résultats concordent au point de vue ordre de grandeur avec les récentes mesures expéri- 


mentales de Crittenden, Walker et Ballam. 


Abstract. — Using the method of quasi-real processes, we have calculated the branching ratios 
between the three observed decay modes of the mu-meson : 


DAC RUE AV 
u—e + v+ 
DAC ENNE TN 


ES 
HE Mer. 


these branching ratios being functions of the threshold energy of the emitted photon or the pair 
created. Our results are in rough agreement with the recent experimental measurements of 


Crittenden, Walker and Ballam. 


Introduction. — Une expérience récente de 
Crittenden, Walker et Ballam [1] a permis de déter- 
miner dans certaines conditions les rapports de 
branchement entre les trois modes de désinté- 
gration observés du méson u*, à savoir : 


a) ut et Ÿ + y, 
b) ur et Ÿ ÿ + y 
c) ur er + v + uv 


b) 


TC: 


Ces trois processus sont essentiellement repré- 
sentés par les trois diagrammes de la figure 1. En 


toute rigueur, il faudrait ajouter pour chacun des 


effets (b) et (c) un diagramme supplémentaire où le. 


photon (réel ou virtuel) serait émis non par l’élec- 
tron, mais par le méson mu; toutefois, étant 
donnée l’importante différence de masse entre le 
méson mu et l’électron, la contribution de ces dia- 
grammes peut être considérée comme négligeable. 

La méthode d’approximation dite des processus 
quasi réels [2] se prête de façon particulièrement 
simple au calcul des rapports de branchement 
entre les trois diagrammes considérés. 


Taux de branchement entre les processus 
u—e+v+v + yet um —e + v + v. — Pour 
comparer les diagrammes b) et a) de la figure 1, 
nous avons à utiliser la formule (17) de la réfé- 
rence [2], formule donnant le spectre correspon- 
dant à un diagramme élémentaire du type de la 
figure 2 (émission d’un photon par un électron rela- 
tiviste) : 


ax € 
; Pmax 2 


9 nl 
PE Rare (in É ;) 
T n k Heure 


ne k a. 
où E et m sont l'énergie initiale et la masse de 
l’électron (nous avons pris c — 1), k et o l’énergie 
et l’angle d'émission du photon. 

Nous utiliserons d’autre part la formule de 


Michel [3], donnant le spectre d’énergie de l’élec- 
tron créé dans le processus a) (en adoptant un 
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système de référence où le méson y est au repos) : 
LE? L 
ÎLE) dE = ya LO(W — E) + 26(4E — 8W)] dE, 


où W est l’énergie maximum de l’électron (environ 
la moitié de l’énergie de masse du méson pu); pa 
été défini à l’origine comme un paramètre compris 


£* À 
0 


E- 
FIG, 2 


entre 0 et 1, dépendant de la nature du couplage. 
Nous adopterons la valeur 6 = 0,75 couramment 
admise à l’heure actuelle. La formule devient 
ainsi : 
2 
HE) dE = 2 (8W— 2E) dE. 


Ce spectre étant normalisé à l’unité, le rapport de 
branchement Ôpa — 6/0, sera donné simplement 
par : 


PE JE T {(E) dE J, PE pi: à) dE. 


Il convient de définir les limites d’intégration : 
k, est évidemment le seuil expérimental de détec- 
tion du rayonnement de freinage interne ; dans un 
calcul précédent [2], nous avons supposé que ce 


._ seuil était de l’ordre de quelques masses électro- 


niques seulement, ce qui nous avait permis de sim- 
plifier beaucoup les calculs ; nous renoncerons ici à 
cette hypothèse. Nous pouvons prendre d’autre 
part Æmax — Æ, et par ailleurs E4 = k4. 

Nous devons également discuter du choix de quax, 
qui comporte un certain arbitraire. Nous prenons 
habituellement œuax & 1 dans les applications de 
la méthode des processus quasi réels. Ce choix 
suppose toutefois que le seuil d'émission du photon 
est relativement bas (autrement dit que k est 
essentiellement petit par rapport à Æ). Si ce n’est 
pas le cas, un choix plus correct serait 

max Le (E AE k) LE 

(tenant compte du fait que les angles 0 et œ de 
la figure 2 doivent tous deux rester petits, pour 
que notre formule soit valable). Cependant, nous 
introduirions ainsi des complications considérables 
dans le calcul. Nous préférons donc prendre 
max — 1 tout en admettant que ce calcul cesse 
d’être valable pour 4, > 0,5 W. 
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Dans ces conditions, l’intégration sur 4 donne : 
5 VND OL BRUT ENTER RL CUS 
in USE k\ dk HE (ny—;) (np: nes co) 

Posons : 
De FT, To —= kolW. 


Compte tenu de : In Wm = 4,6, il nous reste 
pour l’expression de na : 


ka fl 


Spa = e (3x? — 2x5) (In x — 4,1) 


Lo 


(in z— In mn. +2) dæ. 


L'intégration, directe mais assez longue, donne 
(en posant À — In (1/x,)) : 
Spa = LA(1,84 — 0,3328 + 0,084) 


— 9,08 + 2,890 — 1,422 + 0,9078 — 0,294]. 


La courbe de la figure 3 donne les valeurs ainsi 
calculées de x pour k, compris entre 1 MeV et 
25 MeV. 


Ôba ne 
5 


10: 


2 (MeV) 
0 5 10 15 20 25 


TC: 


Taux de branchement entre les processus 
u—e+v+v+e-+et et pm —e + v + v. — 
Pour comparer les processus c) et a), nous avons à 
considérer à la fois un diagramme élémentaire du 
type de la figure 2, et un autre diagramme élémen- 
taire du type de la figure 4 qui représente la créa- 
tion d’une paire électron-positron par un photon 
de haute énergie (4 > m). À ce dernier diagramme 
est associé le spectre (voir 2, formule 35) : 


P(k, E*) dE* 


ü (RUES 4\ LE*X2 + (k— E*)?] dE* 
ne Pmax j) E 1] Q 


kè 


54 
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Nous avons intégré cette formule dans un travail 
précédent [4], obtenant ainsi le coefficient de con- 
version interne d’un photon de haute énergie : 


k 2 & k 19 
* | Mc FR nes 
p P(K, E*) dE* — (in S). 


ss 


Fi. 4. 


Nous avons donc : 
W 
Po fo + f(E) dE 


ke 
[ PEL DIE 4) dk [ P(K, E*) dE* 
(2 m >. 


2& W kmax k 19 
= f. HE) dE [PE à (in = — 53) 


En intégrant sur k (avec les mêmes conventions 
pour les limites que plus haut), et-en posant à 
nouveau: EN; = RIW; À Mn(l/r), 
nous obtenons : 


2a 
da = — 
ca 3x 


ER (3x2 — 2x5) (In x + 4,1) 


TH Lo 
(inte—ÿine A +5 +4)] + 3,02 a L 
2 A z a? 


où l’on a défini : 
A=—5% 2: EURE CES 1) À 
L'intégration finale donne : 
Sea — (8042/3772) [A2(— 0,92 + 0,3428 — 0,09x4) 
+ (5,55 — 2,89x0 + 1,422$ — 1,83a$ + 0,544) 
— 3,89 + 5,79z0— 3,542$ + 2,552 — 0,874]. 
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La courbe de la figure 5 donne les valeurs ainsi 
calculées de Sa pour k, compris entre 3 MeV et 
25 MeV. | 


Comparaison avec l’expérience. — Le tableau 
ci-dessous permet de comparer nos valeurs pour les 
taux de branchement avec les valeurs observées, 
dans les conditions de l’expérience citée [1], 


Obs. Théor. 
1) Sba(ko — 10 MeV) (1,4 + 0,4).10—2 1,31.10—2 M 
2) Spa(ko = 20 MeV) (3,3 + 1,3).10—$ 5,3 .10—% 
3) Sca(ko = 10 MeV) (2,2 + 1,5) 40702 :23,9100 


On constate que nos résultats concordent, du 
point de vue ordre de grandeur, avec les valeurs 
expérimentales. Dans le cas 1) l’accord est excel- 
lent. Quant au cas 2), où le seuil k, est déjà assez 
élevé, le choix que nous avons fait de omax (Voir la 
discussion plus haut) peut expliquer le léger désac- - 
cord. Enfin, pour d.., la valeur théorique se trouve 
également légèrement au-dessus de la limite expéri- 
mentale ; ceci n’est pas surprenant, étant donné 
que notre calcul comporte ici deux approximations 
successives. 

: Nous tenons à signaler que le présent travail a 
bénéficié de l’aide du Commissariat à l’Énergie 
Atomique, ainsi que du Centre National de la 
Recherche Scientifique. 


Manuscrit reçu le 7 juillet 1961. 
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COLLOQUE « LASERS ET PHÉNOMÈNES DE COHÉRENCE » 


organisé par la 
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Ce Colloque, qui s’est tenu dans l’Amphithéâtre de Physique de la Sorbonne, a rassemblé près de 300 


personnes. 


La séance du matin a porté sur les problèmes de cohérence spatialé et temporelle, celle de l’après-midi 


sur les Lasers. 


Quatre rapports généraux, ainsi que des communications originales plus brèves, ont été présentés. 


Rapports. — Problèmes de cohérence temporelle, par P. AIGRAIN. — La cohérence partielle en Optique 
Instrumentale, par À. MarÉcHAL. — Masers hertziens, par P. GRiver. — Les masers optiques, par A. Kas- 


TLER. 


Un résumé de ces rapports, ainsi qu’un Compte-rendu plus détaillé du Colloque, paraîtra dans le Bulletin 


de la Société Française de Physique. 


Communications. — On en trouvera le texte ci-dessous. 


SUR UNE DÉFINITION MATHÉMATIQUE 
DE LA COHÉRENCE 


par Pierre VERNOTTE (Paris) 


Pour autant qu’on puisse la définir d’un mot, la 
cohérence, en optique, se rattacherait, de plus ou 
moins près, pour diverses ondes, à une communauté 
d’origine. La cohérence mathématique dont je veux 
parler signifierait l’appartenance de diverses données 
à une même loi, ce qui est tout différent mais marque 
cependant quelque cousinage entre les deux notions. 

Une même loi, cela suppose qu’il y a une infinité 
de données à cohérer, car on pourrait toujours écrire, 
fût-ce par un polynôme, une relation entre des 7; en 
nombre fini. Et ceci m’amène à dire qu'il ne peut y 
avoir incohérence ou désordre qu’à une échelle infinie, 
et l’on rejoint, ainsi, le point de vue aléatoire. 

J’ai donc dit que la cohérence mathématique entre 
une infinité de données signifierait que toutes ces 
données rentraient dans une seule et même loi. Mais 
on peut concevoir une infinité de lois pour relier une 
infinité quelconque de données, de sorte que, si l’on 
veut réserver une notion de cohérence, le problème 
sera de choisir un genre de loi caractéristique d’une 
suite. Cette loi apparaîtrait sous la forme d’une ex- 
pression mathématique ayant un sens quelle que soit 
la variable x, et non pas seulement lorsque x reçoit 
les valeurs entières correspondant aux numéros d'ordre 
successifs des données y. Autrement dit, le problème 
est un problème d’interpolation. On doit donc chercher 
linterpolation la plus naturelle, mais le problème re- 
bondit : que faut-il entendre par là ? Si je considère 
une expression valant zéro pour x = 0,1,2,3,...n,.…., 
faut-il interpoler par y — 0, ou par y = sin xx ? Car 
.es deux lois seraient raisonnables. 


Il est donc nécessaire d'introduire une condition par 
ticulière, et la condition à laquelle je me suis arrêté 
est celle-ci : si, ayant fait subir aux donnnées y; d’abs- 
cisse entière & une certaine transformation o(:;), on 
applique à ces données transformées la même méthode 
d'interpolation qui avait conduit à la fonction F(x), 
on obtient la transformée o[F{(x)], c’est qu’on a réalisé 
l’interpolation que j'appelle idéale [1], et c’est l’exis- 
tence d’une telle interpolée idéale F(x) qui sera, par 
définition, la preuve de l'existence d’une loi reliant les 
diverses données y;, et donc la preuve de la cohérence 
de la suite y. | 

Ainsi, soit la suite des factorielles y — à ! On réalise, 
par les produits infinis indiqués par Gauss, une inter- 
polation F(x) que l’on sait être égale à l(x + 1). Si 
l’on forme les carrés des factorielles en posant 
o(yi) — y}, le même calcul conduit à l’interpolée 
[F(x)}?, donc la fonction Gamma réalise l’interpolation 
idéale entre les factorielles, qui constituent donc bien 
une suite cohérente. On peut montrer que, dans le cas 
de y; = 0, c’est y = 0 et non y = sin rx qui réalise 
l’interpolation idéale. 

Cela étant, supprimons une des données, soit yx ; 
construisons l’interpolation idéale entre les autres, 
soit F(x) et soit donc F(k) la valeur interpolée rem- 
plaçant la donnéefsupprimée. Si F(k) n’est pas égal 
à Yr, C’est que yx était un terme aberrant, sans cohé- 
rence avec les autres y;. Plus généralement, suppri- 
mons diverses données Yx, Yr, Ym Yn» +, Et CONStruISONS 
l’interpolée idéale entre les ordonnées conservées. Si 
l'on à y Æ F(k), y À F(1), Yn Æ FM), …, Ye, Yu 
Yms -, SOnt étrangères à la suite cohérente définie 
par les autres, et ont le caractère de quantités arbi- 
traires. 

On parviendrait doné, ainsi, à donner un sens à 
l’idée vague de cohérence, et, par opposition, à l’idée 
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d’arbitraire ; et cela pourrait aller très loin, si loin 
même que je ne veux pas insister : on entrevoit l’appli- 
cation au problème du,déterminisme, de la liberté, … 
Je dois avertir, toutefois, que si j’ai bien donné [1] une 
formule pouvant servir de base au calcul de l’interpo- 
lée idéale, la solution effective de ce problème pré- 
sente actuellement de grandes difficultés. 


L'idée de cohérence ainsi introduite a toute la. 


richesse des conceptions mathématiques, l’imagina- 
tion du mathématicien n’étant pas freinée par toutes 
les contingences qui entravent le physicien soucieux 
de garder le sens du concret. En vertu de cette richesse, 
il pourra se présenter des manières diverses de conce- 
voir la cohérence. Car l’interpolée idéale (même après 
que l’on aurait éliminé de la suite les termes visible- 
ment aberrants, peut ne pas exister au. sens strict, la 
fonction qui tendrait à la réaliser n’étant pas alors 
univalente — ce serait le cas de la fonction (— 1) et 
lPinterpolation ne pourrait alors être effectuée avec le 
caractère d’idéalité, c’est-à-dire la conservation par 
transformation, que si l’on fait appel simultanément 
à plusieurs branches de la fonction que l’on se permet- 
trait de prendre l’une pour l’autre [2]. Cette circons- 
tance se présentera, par exempie, avecles suites ondulées. 
L'application pratique essentielle se rapporte à 
l’étude de l’allure à l’infini des développements en 
série, convergents et divergents. En effet, moyennant 
une généralisation éventuelle de la notion de valeur 
à l'infini, qui, d'ordinaire, n’a pas à jouer dans les 
applications, on peut énoncer qu’il existe dans le plan 
complexe un argument tel que la série entière 


A9 — A3 3 + A9 2? — A3 2 + ... (1) 
soit nulle pour [2] = co, si la suite À; est cohérente. 


Il en résulte que si la suite est telle qu’il ne puisse s’agir 
que de valeurs réelles, l'expression 
A, 3 — À3 2° + A3 28 — ... 

(on peut énoncer un résultat plus général quand cette 
expression est multipliée par 2?) a pour limite, pour 
z — co, la quantité 4, obtenue en extrapolant idéa- 
lement à z — 0 la‘suite cohérente À;, ce qui est sou- 
vent facile d’une façon au moins approchée. A lopposé, 
pour exprimer que l’on n’a introduit, dans la suite À;, 
aucune valeur aberrante, on n’aura qu’à écrire que la 
série (1) est nulle à l'infini. 

Ce problème d’étude à l’infini se pose, en particulier, 
quand on essaie de ramener l'intégration d’une équa- 
tion aux dérivées partielles de la physique, en x et £, à 
celle d’une équation simplement différentielle en z, 
l’argument z unique groupant x et ? par la formule 
z = x"]t. Cette stylisation peut être permise quand les 
conditions aux limites s'expriment : pour x — 0 quel 

ue soit {(z — 0) ; pour { — 0 quel que soit x(z — co). 

n en tire la particularisation cherchée de la solution, 
obtenue [3] sous la forme d’une série [1] dont le pre- 
mier terme a une valeur B, aberrante par rapport 
aux coefficients suivants A;-9 qui forment, eux, une 
suite cohérente : à cette fin, on écrit que B, est tel que 
la somme B, + A, ait la valeur imposée, à l'infini, àla 
solution cherchée. 
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ÉTUDE EXPÉRIMENTALE DE LA LUMIÈRE 
ÉMISE PAR UN LASER A RUBIS 


Par M. PauTHIER, R. GAUTIER, S. DEINESS, 
et G. AMAT 


1. Dispositif expérimental. — Avant d'exposer en 
détail les observations qui ont pu être faites sur la 
lumière émise par un laser à rubis, il est souhaitable 
de donner quelques indications sur le dispositif expé- 
rimental qui est d’ailleurs analogue à celui décrit par 
Collins et alii (Phys. Rev. Letters, 17 oct. 1960, p. 303)... 
Le rubis est un cristal de corindon (Al, O;) où quelques 
atomes de chrome se sont substitués à des atomes d’alu- 
minium. Le système de cristallisation est rhomboé- 
drique. Du point de vue optique, c’est un cristal 
uniaxe. L'indice ordinaire est 1,769, l’indice extraor- 
dinaire 1,760. 

Les rubis qui ont été utilisés sont roses, leur concen- 
tration est d'environ un atome de chrome pour 2 000 
atomes d'aluminium. Les atomes de chrome sont en 
moyenne distants les uns des autres d'environ 30 ang- 
strôms, ce qui est petit par rapport à la longueur d’onde 
de la lumière émise (7 000 angstrôms). 

Les cristaux ont été taillés en cylindres d’environ 
40 à 50 mm de long et 5 à 8 mm de diamètre. Les deux 
faces extrêmes sont polies, planes à moins d’une frange 
près et parallèles entre elles à une fraction de minute 
près. L’axe du cylindre est, soit parallèle à l’axe op- 
tique du cristal, soit perpendiculaire. 

Des mesures de diffraction aux rayons X ont montré 
que dans les meilleurs échantillons, il peut se produire 
des variations fortuites d’orientation de l’axe de 
l’ordre de 40. 

Les deux faces du cylindre sont argentées par évapo- 
ration sous vide, l’une de facon à être complètement 
opaque, l’autre de façon à avoir une transparence 
d'environ 10 %,.. Le cylindre est ensuite placé au centre 


IMPULS1Iori 
D'AMORCAGE 


= 
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d’un tube hélicoïdal en quartz rempli de xénon. On 
décharge dans ce tube, en l’ionisant au moyen d’une 
électrode auxiliaire, une batterie de condensateurs sui- 
vant le schéma de la figure 1. On observe la lumière 
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émise à travers la face semitransparente du rubis sur 
un écran placé perpendiculairement à l’axe du cylindre 
et à environ 2 mètres. 


II. Premiers résultats. — Pour de faibles énergies 
dissipées dans la charge, on n’observe qu’un éclaire- 
ment uniforme de l’écran dû à la fluorescence du rubis. 
Si, à capacité constante, on fait croître la tension de 
décharge, il existe pour chaque rubis et chaque tube à 
éclair une valeur bien déterminée de cette tension pour 
laquelle on voit apparaître sur l’écran une tache rouge 
extrêmement intense d’un diamètre de l’ordre de 10 
à 20 mm suivant les rubis. Si on place un détecteur 
photoélectrique à l’endroit où apparaît la tache, on 
peut relever les oscillogrammes suivants placés dans 
l’ordre des tensions de décharge croissantes (fig. 2). 

1. Fluorescence normale du rubis avec décroissance 
exponentielle qui correspond à une durée de vie 
moyenne à l’état excité de 4 millisecondes. 

2. Au maximum de la fluorescence, une légère pointe 
apparaît annonçant ainsi le seuil de l’effet laser. 

3. Pour une tension supérieure de 1 %, l'effet laser 
se traduit par une très brutale discontinuité de la 
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courbe. L’ordonnée maximum est des milliers de fois 
plus grande que le maximum de la fluorescence, ce 
qui correspond à la concentration de l’énergie lumi- 
neuse du rayonnement dans un angle solide très réduit. 
La fluorescence est encore perceptible après la dispari- 
tion de l'effet laser, quoique plus faiblement que s’il 
n'avait pas eu lieu. 

Un balayage plus rapide de l’oscilloscope montre que, 
pendant que se produit l'effet laser, l’énergie est émise 
sous forme d’impulsions brèves, bien individualisées 
d’une durée de l’ordre de la microseconde. 

4. Pour une énergie de pompage supérieure de 20 % 
à l’énergie du seuil, la durée et l’intensité de l’effet 
laser sont fortement accrues. On n’observe alors prati- 
quement plus de fluorescence après la disparition de 
l'effet laser. 

Un balayage plus rapide de l’oscilloscope montre une 
forêt d’impulsions se recouvrant les unes les autres. 


III. Étude de la face émissive du cristal. — Au 
moyen d’un objectif, on peut former l’image agrandie 
de la face émissive du cristal sur un papier sensible. 
Lorsque l’énergie de pompage dépasse légèrement le 
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seuil, on constate un ensemble de taches distinctes 
situées de préférence au milieu du rubis (fig. 3) et d’un 
diamètre moyen de quelques dizièmes de millimètres. 
Lorsque l’énergie de pompage est bien supérieure au 
seuil, on constate alors qu'il n’y a qu’une grosse tache 
centrale (fig. 4), mais dont la structure granuleuse 
laisse penser qu’elle est constituée par la superposition 
dans le temps et dans l’espace d’une multitude de 
petites taches élémentaires. Cette supposition rend 
compte de la directivité relativement faible (de l’ordre 
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du 1/100 de radian) du faisceau observé. Il est donc 
possible de considérer le phénomène obtenu comme la 
superposition au centre du cylindre d’un grand nombre 
de fibres lasers parallèles à son axe. Si on diaphragme 
la face du rubis de façon à ne laisser rayonner qu’une 
petite ouverture de 0,02 mm de diamètre, l’oscillo- 
gramme enregistré montre la présence d’un certain 
nombre d’impulsions séparées, ce qui correspondrait 


à des oscillations de relaxation dans l’émission d’une 
fibre. 


INSTa 


IV. Étude du rendement lumineux. — Si on intègre 
le flux lumineux émis au cours de l’effet laser, on cons- 


tate que l’énergie lumineuse concentrée dans le fais- 


ceau est environ 1/2 joule (fig. 2). La puissance dissi- 
pée dans le tube à éclair peut être très variable suivant 
la disposition du tube et du rubis. Admettons qu’elle 
soit de l’ordre de 500 joules. On voit que le rendement 
est de l’ordre de 1/1000 à peine. 

Toutefois, un examen plus serré fait apparaître que 
ce piètre rendement est principalement dû à l’insuf- 
fisance de l'efficacité du dispositif d’illumination, mais 
que le rubis lui-même a une efficacité tout-à-fait con- 
venable. On peut diviser le rendement total en 4 ren- 
dements partiel 

1. Rendement lumineux de la lampe à éclair, c’est- 
à-dire rapport de l’énergie lumineuse rayonnée dans 
la bande spectrale utile à l’énergie électrique totale 
dissipée, soit R.. 

2. Rendement géométrique de la lampe, c’est-à-dire 
rapport du flux lumineux tombant sur le rubis au flux 
lumineux total rayonné, soit A. 

3. Rendement d’absorption du cristal, c’est-à-dire 
rapport du nombre d’atomes excités au nombre de 
photons verts incidents, soit À. 

4. Rendements d'émission du cristal, c’est-à-dire 
rapport de l’énergie en lumière rouge émise en laser 
à l’énergie absorbée en lumière verte, soit R4. 

R;, R:, R; sont déterminés par la technologie et 
susceptibles de grandes améliorations. Leur produit 


est de l’ordre de 2 %. R, est plus spécialement le ren- 
dement laser. Comparons ses valeurs expérimentales 


et théoriques. 

Dans un cylindre de rubis, il y a environ 2.101 
atomes de chrome. La fréquence de pompage est de 
6.10" et l’énergie nécessaire au pompage est donc de 
Pordre de 8 joules. La fréquence d'émission est 4,3.1014 
et l’énergie maximum qu’on peut penser recueillir est 
environ » joules. On voit donc que le rendement expé- 
rimental est de l’ordre de 10 % du maximum théo- 
rique, ce qui peut être considéré comme bon pour un 
dispositif qui est encore en étude dans les laboratoires. 

La cause la plus évidente de perte est que seule une 
fraction inférieure à 50 % de la face du rubis rayonne. 
Les autres causes sont : l’absorption et la diffusion 
parasites dues aux défauts du cristal ; les pertes aux 
réflexions sur les faces latérales et particulièrement la 
perte à travers la face semi argentée ; la fluorescence 
qui concourt avec l'effet laser pour désexciter les 
atomes et qui se produit seule tant que le pompage n’a 
pas atteint la valeur critique. 


V. Étude de la polarisation de la lumière émise. — 
Deux cas sont à considérer : barreau taillé perpendi- 
culairement à l’axe optique du cristal et barreau taillé 
parallèlement. Dans le premier cas (fig. 5), si l’on place 
devant ce barreau un prisme de Wollaston, on observe 
pendant l'éclair deux taches lumineuses sur un écran. En 
tournant le Wollaston, on obtient l’extinction complète 
de l’une des deux taches. Pour tous les barreaux qui ont 
été observés jusqu'ici, la tache qui reste lumineuse 
correspond à la vibration polarisée rectilignement dont 
le vecteur électrique est perpendiculaire à l’axe op- 
tique. 

ve d’autres termes, l'effet laser se produirait préfé- 
rentiellement, sinon exclusivement, suivant la vibra- 
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tion ordinaire. Notons que si on émet de la fluorescence 
pure en déchargeant dans le tube une énergie inférieure 
au seuil, on s’apervoit que la lumière de fluorescence 
est partiellement polarisée de la même façon. On peut 
proposer une explication simple de ce résultat en re- 
marquant que, dans le rubis, l’absorption de la lu- 
mière à la longueur d’onde de la fluorescence est 
moindre pour la vibration ordinaire que pour la vibra- 
tion extraordinaire. 


AXE OPTIQUE 
TACHE LUMINEUSE 


| BI-RÉFRINGENT 


RUBIS 


DIAPHRAGME TACHE ÉTEINTE. 
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Dans la compétition qui s’établit entre les différents 
modes au départ de l’oscillation maser, ceux qui sont 
établis sur la vibration extraordinaire sont nettement 
désavantagés et sont ainsi étouffés par les autres. 

Dans le cas maintenant d’un barreau parallèle à 
l'axe optique du cristal, le montage précédent donne 
deux taches d'intensité égale, quelle que soit l’orien- 
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tation du Wollaston. L’interposition d’une lame 1/4 
d’onde pour la lumière de fluorescence ne change rien 
au résultat. La lumière n’est donc ni polarisée rectili- 
gnement dans son ensemble, ni polarisée circulaire- 
ment dans son ensemble. On peut penser que chaque 
impulsion élémentaire qui compose la lumière émise 
possède une polarisation déterminée. Il faut donc ana- 
lyser le phénomène dans le temps (fig. 6). On dia- 
phragme donc la face émettrice du rubis, de façon à 
n’avoir que quelques impulsions séparées. On place 
deux cellules photoélectriques à l’emplacement de 
chacune des taches du montage précédent et on photo- 
graphie sur un oscilloscope à deux faisceaux l’oscillo- 
gramme fourni par chaque cellule. De façon à éviter 
toute cause d’erreur avant chaque mesure, on fait un 
essai préalable : on place en avant du Wollaston et 
à 459 de ses directions un polaroïd. On doit alors obte- 
nir la même courbe sur les deux faisceaux. On retire 
ensuite le polaroïd et on refait un enregistrement. 

Les oscillogrammes obtenus sans polaroïd présentent 
alors des différences notables dans l’amplitude des 
impulsions homologues, certaines impulsions ayant 
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une amplitude quasi nulle sur l’un des faisceaux (fig. 7 
et fig. 8). Ces résultats semblent bien montrer que la 
lummièré émise reste polarisée d’une certaine façon 
pendant la durée d’une impulsion. La nature de cette 
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polarisation est plus difficile à préciser. Une polarisa- 
tion rectiligne dont la direction de polarisation serait 
variable avec chaque impulsion élémentaire ne serait 
pas incompatible avec ces résultats. M. Kastler nous 
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a fait remarquer qu’une polarisation initialement cir- 
culaire à l’intérieur du cristal serait vite déformée si le 
cristal possédait même une très légère biréfringence 
dans la direction de propagation. En effet, pour ne pro- 
voquer sur un trajet de b cm qu’un retard inférieur 
au 1/10 de longueur d’onde, il faudrait supposer que la 
biréfringence est inférieure à 10. Or ceci est bien 
inférieur à ce que l’on peut pratiquement obtenir à 
cause par exemple des variations d'orientation de l’axe 
optique ou des contraintes internes. 


VI. Conclusion. — Cette étude expérimentale n’est 
qu’un début d'investigation. De très nombreux points 
restent à éclaircir. Le problème expérimental le plus 
important consiste à démêler les impulsions élémen- 
taires qui composent la lumière émise et à étudier sépa- 
rément leurs propriétés. Il est probable que ces pro- 
priétés se révèleront beaucoup plus près de la cohé- 
rence idéale que les propriétés globales de la lumière 
émise n’ont permis de le montrer jusqu’à présent. 
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INFLUENCE DE DIVERS PROCÉDÉS 
D’EXCITATION OPTIQUE 
SUR'LES POPULATIONS 
DES NIVEAUX D'ÉNERGIE 
DES IONS Crô+ DU RUBIS 


Par François GiREs et Guy MAYER 
Département de Physique Appliquée, C. $. F., Corbeville. 


Nous utilisons deux types de dispositifs pour mesurer 
optiquement le degré d'occupation des divers niveaux 
électroniques des cristaux fluorescents en présence 
d’excitations brèves ou périodiques. 


A. Excitations périodiques. — La figure 1 montre 
le schéma de principe du dispositif conçu pour mesurer 
des effets relativement faibles dus à une excitation 
périodique. 

Après passage dans le monochromateur M, la lu- 
mière issue de la source s continue et stable est modulée 
en intensité par l’assemblage des deux polariseurs P 
et À, de la lame quart d’onde 1 et d’un cristal transpa- 
rent C- rendu périodiquement biréfringent par une vi- 
bration mécanique. La fréquence de modulation Q 
peut être choisie dans une gamme de 50 à 500 kilo- 
cycles parmi les fréquences propres de résonance méca- 
nique du cristal. 

Nous avons utilisé un cristal de PO4 H? NH, trans- 
parent entre 2 500 À et 13 000 À, dont la piézoélec- 
tricité permet d’exciter commodément les vibrations, 
mais n'importe quel cristal transparent collé sur n’im- 
porte quel cristal piézoélectrique excitateur est à peu 
près aussi efficace, 
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Le cristal est muni de plusieurs paires d’électrodes ; 
l’une d’entre elles sert à l'excitation ; les autres servent 
à ajuster la fréquence de modulation et à en contrôler 
Pamplitude dont la stabilité fixera en définitive la 
qualité des mesures. 


Fig. 1. — Dispositif de mesure de Cépopulations faibles et 
périodiques. s : lampe à filament de tungstène ; M: 
monochromateur ; L: lentilles ; P et A: polariseurs 
rectilignes ; Z : lame quart d’onde ; Cr : cristal de phos- 
phate d’ammonium; R': parallélépipède de rubis ; 
S : lampe au xénon 150 watts ; M: miroir sphérique ; 
C : cellule photoélectrique ; A, : amplificateur sélectif ; 
D, : détecteur synchrone ; V : voltmètre. 


La lumière passe ensuite dans l’échantillon étudié R, 
éclairé périodiquement par une source S, à la fré- 
quence &, choisie nettement inférieure à la fréquence Q 
de la modulation due au cristal biréfringent. Si par 

exemple S est une lampe à haute pression de xénon 

alimentée sur le secteur, la lumière produite contient 
une forte composante à 100 périodes. Concentrée sur 
un rubis R par un miroir sphérique M, elle y créera 
une dépopulation de l’état fondamental périodique 
de même fréquence. 

Si la lumière modulée à la fréquence Q a une lon- 
gueur d’onde correspondant à une transition relative 
au niveau fondamental, son amplitude se trouvera 
modulée à la fréquence 100 (fig. 1). 

La profondeur de cette modulation et son déphasage 
par rapport à la lumière excitatrice peuvent être mesu- 
rés à l’aide de l’ensemble cellule photoélectrique C, 
amplificateur sélectif A, détecteur synchrome D, et 
voltmètre V. 

Le déphasage donne des indications relatives aux 
durées de vie + quand wr = 1. La profondeur de modu- 
lation indique le degré de dépopulation du niveau fon- 
damental. 

En effet, considérons le cas simple où un seul groupe 
de niveaux a une durée de vie notable ; cela semble être 


le cas du rubis où les niveaux Æ(?2Æ) et 2A(2ÆE) sont 
caractérisés par la même durée de vie [1] 3 1055 à 
la température plus élevée que l’ambiante où nous 
opérons. 

Si ©, est la puissance moyenne de la lampe, le flux 
émis sera (x/2)p0lcos wt/2| où w/2 est la pulsation 
du secteur et sa composante périodique la plus impor- 
tante sera (2/3) ©, cos wt. Si N, est le nombre d’ions 
Cr$+ par cm$ de rubis, N le nombre de Cr*+ dans l’état 
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fondamental, »r le nombre de Cr%* dans les états EÊ(E) 
et 2A(2E), on aura N = N,—net 


(dnfdt) + (nfr) = (213) x@0 cos ot + A6, (1) 


où À est une constante dépendant de la géométrie de 
l'expérience, de la nature de la source S et des propriétés 
du rubis. 

La dépopulation du fondamental est périodique ; 
elle oscille de la valeur moyenne A9, + avec une ampli- 

2 ÀPoT 
tude de à TS — 

VI + ox? 
à la lumière excitatrice de arc tg or. 

Si la lumière transmise (à la fréquence Q) est sus- 
ceptible d’exciter les ions Cr°* dans leur niveau fon- 
damental, on aura en appelant 6 la section efficace 
de cette excitation 


I = I,e—No et AIJI = — (ANJN) Nox. (2) 


et elle est déphasée par rapport 


La quantité Nox est fournie par des mesures clas- 
siques d’absorption. La figure 2 représente les gran- 
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Fi1G. 2. — Absorption optique d’un rubis à 0,03 %, de CroO3. 
En trait plein le domaine spectral de nos expériences, 


deurs V, pour diverses longueurs d’onde dans les 
rubis à 0,03 % de Cr? 03 (N, = 10! cm3), le vecteur 
électrique étant perpendiculaire à l’axe optique. Le 
trait plein correspond au domaine où nous avons déjà 
expérimenté ces dispositifs. 

A l’aide d’une lampe à xénon Osram de 150 watts, 
nous avons pu mesurer une dépopulation périodique 
relative d’amplitude 2.10—4 dans un rubis de 
5 X 5 X 12 mm contenant environ 3.1018 ions, ce 
qui correspond à une dépopulation moyenne de 
D Le 

Le rubis ne recevait que 1/12 du flux total de la 
lampe. De ce flux reçu, il n’absorbait que 1/3, compte 
tenu des divers polarisations de la lumière incidente ; 
le facteur t/\/1 + w?r? vaut 1,5.10—88. 

Si on admet [2] que chaque photon absorbé amène 
un ion dans les états excités E(?2Æ) ou 2A(?E), on cons- 
tate que 8 % de l'énergie électrique fournie à une 
lampe au xénon à haute pression peut servir à exciter 
la fluorescence du rubis. 


83% 


B. Mesures d’effets forts en présence d’excitations 
brèves. — Si quelques centaines de joules sont trans- 
formées en lumière par un tube «flash » au xénon en 
un temps petit devant +, on gagne un facteur de l’ordre 
de 10% dans le degré de dépopulation du niveau fonda- 
mental et il peut atteindre quelques dixièmes. 

En fait, avec un flash en spirale entourant des rubis 
de 40 mm de longueur et de diamètre 3 ou 5 mm, les 
conditions géométriques d’illumination sont moins 
bonnes que dans l’expérience précédente et nous 
n’avons obtenu le fonctionnement en «laser » [3] que 
pour des seuils de 4,5 kilojoules (1). 

L'utilisation d’un réflecteur en papier d'aluminium 
nous à permis d’abaisser ce seuil de 17 %. 

Nous nous sommes proposés de mesurer à l’aide 
du dispositif de la figure 3 la population du niveau 
fondamental au’seuil critique de fonctionnement en 
« laser ». 


As 
L se) vers oscilloscope 
|  |---—- - 
Fic. 3. — Dispositif de mesure de la population du niveau 


fondamental en présence d’une excitation optique brèpe. 

S : lampe au xénon de 150 watts ; P et A: polari- 
seurs rectilignes ; l : lame quart d’onde ; Cr : cristal de 
PO* H? NH; L: lentille ; R: rubis ; M : monochro- 
mateur ; C: photomultiplicateur ; As : amplificateur 
sélectif. 


Le parasite électrique créé dans le laboratoire par la 
décharge d’une batterie de condensateurs chargés à 
4,5 kilojoules est une difficulté de l’expérience. Le 
parasite optique dû à la lumière du flash diffusée ou 
réfléchie sur les parois de l’enceinte d’expérience en est 
une autre. La bandefde l’amplificateur sélectif doit 
être assez étroite pour assurer une bonne protection 
contre ces perturbations, mais assez large pour per- 
mettre l’analyse de phénomènes rapides. Nous avons 
choisi 8 kHz. Les lectures se font finalement à l’oscil- 
loscope. 

Beaucoup d'expériences nous sont encore nécessaires 
pour chiffrer l’imprécision de nos mesures, dont voici 
le résultat préliminaire : 

Dans les conditions exactes d’excitation optique 
qui nous permettent d'obtenir l'effet laser de concen- 
tration de la lumière de fluorescence dans un cône d’axe 
perpendiculaire aux faces planes du cylindre de rubis, 
et d'ouverture angulaire inférieure au degré, avec le 
même rubis dépourvu de ses argentures pour permettre 
la mesure de transmission, nous mesurons pendant les 


(1) Les expériences étaient faites avec une capacité de 
1 600 &F. Nous avons depuis obtenu le fonctionnement en 
« laser » avec 750 joules en utilisant une capacité de 80 UF, 
toutes choses égales par ailleurs. 
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flashes une population minima du niveau fondamental 
égale à 79 % pour une excitation inférieure au seuil 
(4,3 kJ), et une population minima de 75 % pour 
pour une excitation supérieure au seuil (4,7 kJ). 

Si le fonctionnement en laser se produisait entre 
deux niveaux simples, il est peu douteux que le seuil 
correspondrait à un appauvrissement du fondamental 


de l’ordre de 50 %,. Si le niveau Æ(2E) qui est double 
existait seul, on pourrait attendre le fonctionnement en 
laser pour une population du fondamental, dégénéré 
à l’ordre 4, s’élevant encore à 66 % ; mais il faut aussi 
approvisionner 2A (2£) qui ne joue pas de rôle dans” 
l’amorçage des oscillations, bien qu’il contribue peut- 
être à les entretenir si l’équilibre de Boltzmann s’éta- 
blit assez rapidement entre les deux niveaux. 
Ceci nous ramène donc à attendre un chiffre de l’or- 
dre de 50 0/4, là où nos premières mesures donnent 75 0/. 
Nous poursuivons nos efforts d'amélioration de ces 
dispositifs pour mieux comprendre les phénomènes du 
rubis et étudier les possibilités d’autres cristaux (?). 


(2?) De nouvelles expériences ont confirmé les valeurs 
données ici pour la population du fondamental au seuil de 
fonctionnement en « laser ». 
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LES SEMI-CONDUCTEURS 
ET LEUR UTILISATION POSSIBLE 
DANS LES « LASERS » 


Par C. BENOIT À LA GUILLAUME et Mme TRrc 
École Normale Supérieure. 


Parmi les diverses possibilités qu’offrent les semi- 
conducteurs pour la réalisation de « Lasers », on peut 
distinguer deux catégories : 

a) Les « Lasers » de type classique, où il faut réaliser 
une inversion de population entre deux niveaux. 

b) Les « Lasers » à deux bosons qui utilisent une 
transition au second ordre entre deux niveaux avec 
émission simultanée d’un photon et d’un phonon. On 
montrera que dans ce cas, il n’est pas nécessaire de 
réaliser une inversion de population entre les deux 
niveaux. 

D'une façon générale, l'intérêt des dispositifs à semi- 
conducteurs résiderait dans la facilité d’excitation, le 
pompage s’effectuerait généralement par injection de 
porteurs en excès à l’aide de jonctions np polarisées 
dans le sens direct. On aurait aussi des possibilités 
immédiates de modulation par contrôle du courant 
dans la Jonction. En contrepartie, le milieu actif con- 
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tient de nombreux porteurs libres, et il apparaît que 
dans bien des cas, absorption due aux porteurs libres 
sera une des causes principales de pertes. 


I. Laser du type classique. -— Nous citerons deux 
possibilités : 


a) TRANSITION ENTRE NIVEAU DE BASE ET NIVEAU 
EXCITÉ D’UNE IMPURETÉ. 

Prenons par exemple le cas de l’impureté As dans 
Ge (fig. 1). L’état de base 1 s de l’électron lié à As+ est 


situé à 14 meV en-dessous du bas de la bande de con- 
duction Æ4 et l’état excité 2p + 1 à 1,6 meV en-des- 
sous de Æ.. La transition optique entre ces 2 niveaux 
est permise ; elle a été observée en absorption vers 
99 u [1]. Le calcul du rapport des populations a2/a1 
des états 2 et 1 en fonction de la densité de trous injec- 
tés p donne : 


d> __p Cns Cp1 + À exp (E1 — E,JRT) 
Le es pCn Cp: + À 


Cn, et Cn, représentent les probabilités pour 
qu’un As+ capte un électron libre dans l’état 1 ou 2, 
Cp;et Cp, les probabilités pour qu’un trou libre se 
recombine avec un électron lié à As* dans l’état 1 ou 2. 

On voit que a, peut devenir supérieur à 1 si 
Cno Cp1lCn: Cpa 1. D’après les calculs de M. Lax [2] 
on peut affirmer que Cn, > Cn,. D'autre part, nous 
avons étudié en détail les transitions entre les trous 
libres et les électrons liés aux donneurs. Dans le cas 
de As, on trouve Cp, - Cp2, à cause de l’existence 
d’une assez forte transition directe entre E, et la bande 
de valence. Finalement Cn,Cp,]Cn; Cp, est probable- 
ment supérieur à 100, mais cependant la valeur de p 
nécessaire pour observer l’inversion de population peut 
être assez considérable. Il est à craindre que l’absorp- 
tion due aux porteurs libres (qui est proportionnelle 
à A2 et à 1/T) rende ce projet irréalisable. 


b) LASER UTILISANT UN SEMI-CONDUCTEUR DÉGÉ- 
NÉRÉ, — Considérons le cas de InSb dégénéré de 
type n (fig. 2). L'étude de la recombinaison radiative 
dans InSb confirme que la transition entre les états 
1 et 2 de même k dans la bande de valence et de con- 
duction est permise [3]. Il est facile de voir que si 
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l’état 2 est au voisinage du bas de la bande de conduc- 
tion, on peut réaliser une inversion de population entre 
les états 2 et 1. Il suffit pour cela que la différence entre 
les positions des quasi-niveaux de Fermi pour les élec- 
trons et pour les trous soit supérieure à la largeur de 
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bande interdite. Par exemple si le niveau de Fermi est 
situé à 5 meV au-dessus du bas de la bande de conduc- 
tion à 4 CK, ce qui correspond à une densité d’électron 
no — 105, on‘trouve que l’inversion de population 
en bas de la bande apparaît pour une densité de trous 
injectés, supérieure à 10! par em. 

On peut espérer que dans ce cas, vu la longueur 
d'onde de fonctionnement 5,3 4, l’absorption par les 
porteurs libres ne sera pas un obstacle insurmontable, 


IT. Laser à deux bosons. — Dans Ge et Si, la re- 
combinaison radiative intrinsèque s'effectue à basse 
température sous forme d’annihilation d’excitons indi- 
rects. Il s’agit d’une transition au second ordre où 
2 particules sont émises ou absorbées : 1 phonon et 
1 photon. Si on désigne par M et N les nombres de 
photons et de phonons par mode, la probabilité d’émis- 
sion induite d’un photon dans ce mode, avec émission 
induite ou spontanée d’un phonon, est proportionnelle 
à M(N +1) P, et la probabilité d'absorption de 
même photon, avec absorption d’un phonon est pro- 
portionnelle à MNP,, P, et P, désignant les nombres 
de places occupées et libres dans la bande d’exciton. 
On voit que l’émission induite sera supérieure à l’ab- 
sorption si P,/P, > NJ(N +1). Or les phonons 
dans Ge ont une durée de vie très courte (mode for- 
tement amorti) et par suite, la population des phonons 
demeure celle de l’équilibre à la température 7, d’où 
la condition P, P, > exp — Aw/kT où io est l’éner- 
gie du phonon. 

Si on représente la population de la bande d’exciton 
par une température effective 7%, on trouve la condi- 
tion 7, > Ex T[ho. Donc on n’a pas besoin de réaliser 
une T7, négative. Par exemple, pour Ge à 40K, on 
trouve 7, > 120 0K ce qui correspond à une densité 
de porteurs injectés extrêmement faible. 

En fait, ce sont les autres causes de pertes qui vont 
fixer le produit np minimum pour l’observation d’émis- 
sion cohérente. Prenons comme exemple une struc- 
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ture possible (fig. 3) où le problème de la sélection des 
modes a été pris en considération. On obtient les fac- 
teurs de qualité partiels suivants : pour les pertes par 
réflexion : 


4/g: = À —r)[27L, 


24 


Surfaces 
argentees 


F6: 


— pour les pertes dues à la diffraction et au placement 
latéral des modes non purement longitudinaux 


4 [do = À? (nx + ny + 1)/arD?, 
— pour l’absorption dans le matériel 
1 [qe = œX[2r (x coefficient d'absorption). 


On choisit Z de façon à ce que l’écartement entre 
deux modes longitudinaux consécutifs (correspondants 
àn.et n, + 1) soit du même ordre que la largeur natu- 
relle de la raie, d’où ZL = 0,4 mm. 

1 /q2 est alors fixé par r. 


Wade Tr sine "1008; 


On choisit 1/g, (pour nx7 — ny — 0) voisin de 1/g, 
de façon à amortir notablement les modes pour les- 
quels nr; ou », sont différents de zéro, d’où D = 0,1 mm 
1 /g, est dù essentiellement à l’absorption due aux por- 
teurs libres (électrons et trous). Cependant à 4 °K, 


Ag, 10/7 


tant que n — p < 1018 par cm. 

La condition d’oscillation obtenue en écrivant que 
l’on à au moins {1 photon induit par mode conduit à 
n—p> 107 par cem*. Il semble donc que le fonctionne- 
ment de ce « Laser » est possible. Les densités de 
porteurs considérées paraissent réalisables, au moins 
en impulsions longues (— 1 ms) dans la structure 
pin envisagée. 


RÉFÉRENCES 
[1] Fan (H. Y.) et Fisner (P.), J. Physique Chem. Solids, 
1959287270: 


[2] Lax (M.), J. Physique Chem., 1959, 8, 66. 
[3] Bexorr (C.) et LavazzarD (P.) à paraître. 
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POSSIBILITÉS DE LASERS 
A SEMI-CONDUCTEURS 


par Maurice BERNARD (*) et Georges DURAFFOURG (*). | 


Dans la plupart des publications consacrées depuis | 


quelques années aux «lasers », le cas habituellement 
considéré, à notre connaissance, est celui de l’émission 


stimulée entre des niveaux localisés. Ainsi, il est bien | 
connu que la condition nécessaire pour pouvoir ampli 


fier des ondes électromagnétiques par émission induite 
est qu’il existe entre les deux niveaux localisés une 
inversion de population c’est-à-dire que le niveau 


d'énergie supérieure soit plus peuplé que le niveau | 
d'énergie inférieure. Cependant, en 1958, P. Aigrain | 


[1] a montré que dans le cas de l’émission de deux 
bosons, l’un stimulé, l’autre thermalisé, cette condi- 
tion n’est plus nécessaire. 


En fait, dans l’approximation à un électron, un état | 
électronique d’un cristal est représenté par une onde | 
de Bloch étendue à tout le cristal. Considérons dans | 


un semi-conducteur deux tels états, l’un de vecteur 


d’onde X; d'énergie Æ, (ki), l’autre de vecteur d’ondek;, | 
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Î 


d'énergie Ee (k;) ; supposons que le premier appar: : 


tienne à un groupe d'états (v) et le second à un : 
groupe (c) ; nous admettons en outre que l’on puisse : 


passer d’un état ÆE, (k;) à un état Æc (k;) par une tran- 


sition mettant en jeu d’une part un photon d'énergie Av, : 
d’autre part part un ou plusieurs bosons (photons où 
phonons) que nous supposerons en équilibre thermody: 


namique avec un thermostat. 

Si les électrons occupant les états de l’un des groupes, 
par exemple le groupe (c), échangent de l’énergie avec 
le thermostat dans un temps court comparé avec le 
temps nécessaire pour que la population de ce 


groupe [2, 3] revienne à l’équilibre par des transitions 


vers des états de l’autre groupe, alors l’occupation 
des états (c) peut être décrite par un pseudo-niveau 
de Fermi F4 associé à ce groupe ; de même l’occupa- 
tion des états (v) est décrite par un pseudo-niveau de 
Fermi F, ; à l'équilibre Fe = Fy = F, où F, est le 
vrai niveau de Fermi. Un calcul assez simple montre 
alors que la condition pour que les transitions induites 
Pemportent sur les absorptions est : 


Fe ES Fy = Rv. (1) 


Cette relation est générale : elle est l'expression du 
deuxième principe de la thermodynamique ; en efet, 
en appliquant un raisonnement déjà utilisé dans le cas 
du maser à trois niveaux localisés [4] on montre la 
relation suivante : 


Fe — Fy = h\y [1 — (T[Tb)], (2) 


où hv, est l'énergie des photons fournis par la transi- | 


tion « pompe » considérée comme source chaude, 7, 


la température de cette source et T la température 


du thermostat qui sert de source froide. La relation (1) 
peut donc s’écrire : 


TE à. ON 


où 1 est le rendement du laser et 1 le rendement de 


(*) Gentre National d'Études des Télécommunications, 
Département Physique-Chimie-Métallurgie, 


LL dun 


Le, LE. 
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Carnot ; la relation (3) montre qu’une machine ther- 


. mique fonctionnant entre deux sources de tempéra- 
à tures 7 et T, et fournissant de la lumière cohérente 


a un rendement inférieur au rendement de Carnot. 

Le résultat précédent peut être appliqué en identi- 
fiant les groupes d’états (v) et (c), soit à la bande de 
valence et la bande de conduction d’un semi-conduc- 
teur, soit à l’une de ces bandes et à un niveau d’impu- 
reté. 

L'intérêt pratique de la condition (1) apparaît dans 
les remarques suivantes : 

19 Transition directe bande à bande : la condition (1) 
suggère dans ce cas l’emploi à basse température d’un 
semi-conducteur à faible masse effective (In Sb ou In 


Vecteur d'onde 


} As) ; la figure 1 montre la transition pompe et la tran- 


sition signal. 
29 La condition (1) sera plus facile à satisfaire si, 
pour un semi-conducteur donné, on peut diminuer Av ; 


» c’est Justement ce qui se produit pour les transitions 
indirectes ; on retrouve là un résultat établi par Ar- 


GRAIN [1] qui, le premier, a proposé en 1958, d'utiliser 
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la lumière de recombinaïson dans le germanium pour 
obtenir une émission de lumière cohérente. 

30 La condition (1) sera encore plus facile à satis- 
faire si on utilise une transition vers un niveau d’impu- 
reté ; la figure 2 montre les transitions pompe et signal 
dans le cas du second niveau accepteur du zinc dans 
le germanium. 


Energie 


Fic.°2. 


Un article plus détaillé sera publié prochainement ; 
des expériences sur In Sb, In As et sur Zn-Ge sont en 
Cours. 
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RÉSONANCE MAGNÉTIQUE NUCLÉAIRE 
A DOUBLE QUANTA DU DEUTÉRIUM 
DE L'EAU ZÉOLITIQUE LOURDE 


Par Avant P Duchos XX PARE, 


Laboratoires d’Électrostatique et de Physique du Métal 
et de Physique Générale, Grenoble. 


Résumé 


Le travail présenté a été effectué sur un monocristal 
naturel d'Edingtonite-alumino-silicate hydraté de la 


famille des zéolites. Nous avons partiellement substitué 
à l’eau du cristal, de l’eau lourde et étudié la résonance 
magnétique nucléaire des noyaux de deutérium. La 
figure de résonance se compose d’un triplet de raies 
d’égale intensité et d’égale largeur (0,2 Oe) et d’un 
doublet de raies identiques (0,2 Oe de large). L’en- 
semble est symétrique par rapport à la raie centrale du 
triplet qui correspond au champ de résonance du 
deuton libre à la fréquence utilisée (6 MHz) [1]. 

Une interprétation quantitative complète du spectre 
est donnée. Le triplet est attribué à la résonance des 
deutons des molécules D,0 et le doublet de HDO, Les 
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raies correspondant aux transitions habituelles (absorp- 
tion d’un quanta d’énergie radiofréquence), sont inob- 
servables par suite du mouvement moléculaire et de 
l'existence du moment quadrupolaire du deuton. Les 
raies observées correspondent toutes à des transitions 
à double quanta (AM = + 2). 
Section Sud-Est, Grenoble, 
Séance du 19 avril 1961. 
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ÉTUDE DE LA MOBILITÉ 
DES IONS ÉCHANGEABLES DANS UNE ZÉOLITE 
: PAR DES MESURES 
DE RELAXATION DIÉLECTRIQUE 


Par P. Ducros, 


Laboratoires d’Électrostatique et de Physique du Métal 
et de Physique Générale, Grenoble. 


Résumé 


Il est bien connu que les zéolites sont des échangeurs 
d'ions. La mobilité des ions échangeables est à l’origine 
de l’important domaine de relaxation diélectrique de 
la chabasie (Ca Al, Si O,>, 6H,0). À température 
ordinaire le maximum de la partie imaginaire de la 
constante diélectrique correspond à une fréquence de 
mesure voisine de 100 cycles par seconde pour la 
chabasie calcique. En faisant varier la température on 

eut mesurer l’énergie d’activation du phénomène : 
0,70 + 0,04eV [1]. 

Par l’échange d'ions Ca++ — 2Na+ le domaine de 
relaxation diélectrique est déplacé vers les hautes fré- 
quences. Ainsi à température ordinaire, le temps de 
corrélation du mouvement des ions (assimilé en pre- 
mière approximation au temps de relaxation diélec- 
trique) passe de 10—2 s pour Cat+ à 10—-7s envi- 
ron pour Na+. Cet effet est interprété par les liaisons 
électrostatiques dépendant principalement de la valence 
des ions. 

La modification de la teneur en eau des cristaux par 
pompage fait varier considérablement le temps de 
relaxation diélectrique. Une explication qualitative en 
est donnée. Elle fait intervenir la mobilité dans le 
cristal des molécules d’eau polaires. 

Section du Sud-Est, Grenoble. 
Séance du 19 avril 1961. 
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SUR LA CHALEUR SPÉCIFIQUE HYPERFINE 
DES MÉTAUX DE TERRES RARES 
DANS LE DOMAINE DE He LIQUIDE 


Par B. Dreyrus, 
Laboratoire des Basses Températures, Grenoble. 


La propriété de l’isotope 3 de l’hélium d’avoir une 
courbe de tension de vapeur saturante supérieure à 
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celle de l’hélium ordinaire a été utilisée pour réaliser au 
Laboratoire des Basses Températures [1], [2] un 
cryostat fonctionnant entre 0,3 ©K et environ 1,2 °K. 


La possibilité de se maintenir, en réglant le débit du | 


fluide, à toute température dans cette gamme, fait de 
cet appareil un instrument particulièrement adapté à ia 
mesure des chaleurs spécifiques. Nous avons ainsi 


entrepris une série de mesures sur les métaux de terres | 


rares. Ces métaux sont en général magnétiquement 


ordonnés avec des comportements d’ailleurs très | 


variés, allant du ferromagnétisme usuel à des struc- 
tures en hélice antiferro ou ferromagnétiques. Le cou: 
plage hyperfin entre le moment nucléaire et la valeur 
moyenne du moment électronique donne lieu à une 
décomposition du niveau fondamental L’anomalie du 
type Schottky, pour la chaleur spécifique, se traduit 


dans le domaine de température étudié par une 


« queue » en A/T?. Les valeurs expérimentalement 
trouvées pour le samarium et l’erbium sont : 


ASS) AIRE 
A (Er) = 25 + 


2 mj mole—1 deg. 


3 mj mole—1 deg. 


Il est intéressant de comparer ce couplage hyperfin 
dans le métal avec celui déterminé pour les sels des 
métaux de terres rares par la résonance électronique. 
En supposant que la valeur du paramètre < 1/r > 


dans la constante de couplage hyperfin est la même. 


dans le sel et dans le métal, on peut rechercher les 
fonctions d’onde de l’état fondamental électronique et 
calculer a priori le rapport des constantes Aser et 
Améta. Lion métallique est soumis d’une part à un 
champ cristallin, probablement faible et d’ailleurs peu 
connu, sauf quant à son signe qui est l’opposé de celui 


du sel, d'autre part à un fort « champ moléculaire ». 
Sous l'influence de ce champ le fondamental est pro= 


bablement un état de M, maximum. 
Nous donnons ci-dessous le rapport théorique ainsi 
calculé et le rapport expérimental des constantes À. . 


TABLEAU 
Tb Dy Sm Er 
A ce RAR 
sel ] 0:87 1.55. 1,89 99 
Amét. exp 
ser | 1 1,86 1,32 3 
Amét./ th 


On peut considérer que l’accord est satisfaisant, vu 
la simplicité du modèle théorique. Il semble, en tout 
cas, qu’il ne soit pas nécessaire de faire intervenir la 
polarisation des électrons de conduction. 


Section du Sud-Est, Lyon. 
Séance du 19 mai 1961. 
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NOUVEL APPAREILLAGE DE MESURE 
DES AIMANTATIONS ; 
APPLICATION AU CRISTAL DE COBALT 


Par R. PAUTHENET, 


Laboratoire d’Électrostatique et de Physique du Métal, 
Grenoble. 


Variation thermique des constantes d’anisotropie. — 
Nous rappelons que dans un cristal de cobalt la direc- 
tion de facile aimantation est suivant l’axe c à basse 
température et dans le plan de base à des températures 


_ supérieures à 600 ©K. Nous représentons l’énergie 


magnétocristalline par une expression de la forme 
K,.sin? 0 + X,.sin4 0 dans laquelle 0 désigne l’angle 
de l’aimantation spontanée J, avec l’axe c. La déter- 
mination des constantes À, et X, a été faite entre 2 
et 700 ©K en utilisant la méthode de W. Sucksmith et 
J. E. Thompson [1]. La constante X, est égale à 
6,6.105 ergs/cm® aux très basses températures, décroît 
lorsque la température croît, s’annule et change de 
signe à 518 0K ; X, décroît régulièrement entre 1,55 et 
0,3.105 ergs/em$ ; l’aimantation à saturation absolue 
est égale à 1,72 magnétons de Bohr par atome de 


_ cobalt. 


Étude du mécanisme de la variation isotherme de 
l’aimentation en fonction du champ. — Désignons 
par o l’angle du champ extérieur H, avec l’axe c ; 
nous avons mesuré pour différents angles ©, dans des 
champs extérieurs jusqu’à 30 000 Oe, les variations 
des composantes Jyj, suivant //., et J1, suivant une 
direction perpendiculaire à },, de l'intensité d’aiman- 
tation résultante J de la substance. Pour ceci, nous 
avons utilisé un électro-aimant de mesure des aiman- 
tations par extraction axiale dans l’entrefer duquel 
est placé un cireuit de mesure constitué par trois 
bobines induites dont les axes sont à 909 les uns des 
autres ; l’un de ces axes est suivant celui de l’électro- 
aimant. Les résultats sont interprétés quantitati- 


- vement par la théorie des phases de domaines élémen- 


taires [2]. En l’absence de champ, à l’intérieur du 
cristal, les aimantations spontanées sont le long de 
l'axe c mais dans des directions opposées ; on divise 
ainsi les domaines élémentaires en deux phases dis- 
tinctes. Faisons croître À, à partir de zéro, en main- 
tenant l’angle © constant ; par déplacement de parois, 
le volume de l’une des phases augmente aux dépens du 
volume de l’autre phase ; les aimantations spontanées 
tendent à s’aligner dans le sens de 4 ; Paimantation 
résultante J sera telle que le champ intérieur 


H; —= fe — nJ 


(n : coefficient de champ démagnétisant) soit perpen- 
diculaire à l’axe c pour que les deux phases puissent 
subsister en équilibre. Ce mode d’aimantation subsiste 
jusqu’à ce que le champ atteigne une valeur pour 
laquelle il n’y a plus qu’une seule phase de domaines 
élémentaires ; pour des champs supérieurs, la variation 
d’aimantation est uniquement due aux rotations de 


|  l’aimantation spontanée de la phase restante. 


Étude du processus du changement d’axe de facile 
aimentation. — Dans un cristal uniaxe, parfait, en 
l’absence de champ extérieur, la position d'équilibre 
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stable de l’aimantation spontanée est définie par le 
minimum de l’expression de l’énergie magnétocristal- 
line K,.sin? 0 + X,.sin4 0. En général, 0 est égal à 0 
ou 7/2 ; cependant, si X, et K, satisfont à l’inégalité 
O<— K,/2K, < 1, J, fait l'angle 
0, = arc.sin /— K,/2K, 

avec l’axe c. Les aimantations spontanées sont alors 
dirigées suivant les génératrices de deux cônes préfé- 
rentiels, opposés par le sommet et d’angle au som- 
met 0,. Cette condition peut être réalisée pour le 
cobalt entre 518 et 598 °K. Nous avons mis en évidence 
expérimentalement cet angle 0,. L’échantillon mono- 
cristallin est placé au centre du système de trois 
bobines induites à axes perpendiculaires ; l’axe c fait 
l'angle © avec 71, ; on suit l’évolution de la compo- 
sante J,. Lorsque la température croît, l'amplitude 
de J1 décroît, s’annule pour 0, — o puis change de 
signe. Connaissant +, en notant la température pour 
laquelle 0, — +, on déduit (0,, T) ; la courbe expéri- 
mentale peut être représentée de manière satisfaisante 
par la variation (are-sin |/— K,/2K,, T). Dans cette 
zone de températures à cône préférentiel, le mécanisme 
d’aimantation s’effectue suivant trois modes successifs. 
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SYMÉTRIE CRISTALLINE 
ET FAIBLE FERROMAGNÉTISME 


Par E. F. BErTaur, 


Laboratoire d’Éectrostatique et de Physique du Métal, 
Grenoble. 


L’auteur retrace d’abord la théorie macroscopique de 
configurations de spins, développée par Dzialo-. 
shinski [1] et Turov [2] à l’aide de l’exemple d’un 
orthoferrite tel que LaFeO;, pérovskite dont le groupe 
d’espace est Pbhnm. Cette théorie consiste à construire 
un hamiltonien //, invariant dans les opérations de 
symétrie du groupe cristallographique. Au lieu d’expri- 
mer {7 dans les variables spins, il est plus commode 
d'envisager certaines combinaisons linéaires privi- 
légiées de spins à savoir celles qui forment la base de 
représentations irréductibles : dans le cas concret des 
pérovskites TFeO;, (T = terre rare ou Y), cette base 
est formée par les quatre vecteurs, F, G, C, À (1) corres- 
pondant 


F=s, Fs, ss Es, 
G=Ssi—5S; +S3—S4 
C=s1 +se —s; —s, 


A = S1 So SRE: S4 (1) 


aux spins s; (J = 1, ..., 4) des atomes de Fe aux 


: : ARE ENS. 1 Tage 
points respectifs 00 555 () 5 ; 0 55 0 5 0. 

Le vecteur G caractérise la configuration colinéaire 
(+—+—) pour laquelle il est maximum. F (1) est 
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maximum pour la configuration (++++) etc. 
Lorsque les interactions d'échange sont isotropes, les 
modes F, G, C et À apparaissent isolément. 

Des interactions anisotropes mélangent ces confi- 
gurations. Toutefois, dans l’approximation d’un hamil- 
tonien d’ordre 2 on peut combiner deux modes seule- 
ment lorsqu'ils appartiennent à la même représen- 
tation, c’est-à-dire lorsqu'ils possèdent les mêmes pro- 
priétés (ou « caractères ») de transformation dans 
toutes les opérations de symétrie du groupe. 

Le tableau 1 résume les résultats ainsi obtenus par 
Turov [2]. Les modes se trouvant sur une même ligne 
peuvent être couplés deux à deux. Cette théorie est 
très puissante pour indiquer les couplages possibles, 
mais elle ignore leur mécanisme. 

Récemment l’auteur a développé une théorie matri- 
cielle, théorie microscopique qui permet de descendre 
aux interactions élémentaires [3]. II a montré notam- 
ment que les vecteurs propres T(K) d’une certaine 
matrice C(k) forment précisément la base de repré- 
sentations irréductibles. On peut scinder la matrice C(k) 
en une partie isotrope Gi(k) et une partie aniso- 
trope Can(k). La partie isotrope contient les intégrales 
d’échange ordinaires. 

Dans la partie anisotrope Ca(k) on peut distinguer 
celle due au 

a) champ cristallin ; 

b) couplage dipolaire ; 

c) couplage pseudo-dipolaire (linéaire dans linté- 
grale d’échange et quadratique dans le couplage 
spin-orbite) ; 

d) couplage antisymétrique de Dzialoshinski [1]- 
Moriya [4] (bilinéaire dans l’intégrale d'échange et le 
couplage spin-orbite). 

La condition générale de l’association de 2 modes 
Ti(k) et Trr(k) est que 


Ti(K) . Ean/k) . Tr(k) 0. (2) 


Les tenseurs d’anisotropie correspondant aux causes 
a) b) c) sont symétriques et donnent lieu à des inva- 
riants symétriques de la forme 


Tia Tug + Tig Tire = (3) 


où 
a ÆP; 


œ, Ê = à, y, z. 

Si Tri représente un mode antiferromagnétique et T1 
un mode f. ferromagnétique (f. f. = faiblement ferro- 
magnétique) nous définirons un f. f. symétrique lorsqu'il 
est dû aux causes a) b) c). C’est le cas de NiF,[5]oùle 
f. f. est dû au champ cristallin. 

Le tenseur correspondant à d) est antisymétrique et 
donne lieu à des invariants de la forme 


Trio Tuig — Tag Tia = r (4) 


c’est le cas de &-Fe,0,[1],[3]. On définira un f, f. anti- 
symétrique lorsqu'il est dû au couplage D-M [1], [4]. 


Tandis qu’à g (3) est toujours associé une aniso-' 


tropie dans le plan («f), r (4) par contre ne contribue 
pas à l’anisotropie basale. On peut donc s’attendre 
a priori à une anisotropie dans le plan de base (1) plus 


(1) C'est-à-dire une direction privilégiée dans le plan Oxy 
perpendiculaire à l’axe c du cristal. 
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forte dans NiF, (due à g (3)) que dans «-Fe,0, où elle 
ne peut être due qu’à des termes d’ordre quatre au 
moins dans À. 

Un autre moyen de distinguer entre f. f. symétrique 
et antisymétrique est la dépendance du f. f. de la 
température (cf. [4]). 

Le f. f. antisymétrique, linéaire dans l'intégrale 
d'échange, « durera » jusqu’au point de Curie tandis 


‘ que le f. f. symétrique dû au champ cristallin (?) ne 


se fera sentir qu'aux basses (azote liquide) et très 


basses (hélium liquide) températures où d’ailleurs les 


deux effets peuvent être compétitifs. 


L'auteur ap ne les considérations précédentes aux 
3 


pérovskites TFeO, où les associations observées par 
Bozorth et al. [6] et attribuables au réseau de Fe sont 
G, Fa aux températures élevées et G;F, aux basses 
températures. 

Il montre que la matrice Caipolaire(K) est diagonalisé 


par la matrice des vecteurs propres (1) de sorte que ” 


l'énergie dipolaire ne donne ici lieu à aucun couplage 
entre modes. Le tenseur pseudo-dipolaire, en se limi- 
tant aux interactions entre proches voisins (1-2 et 1-4) 


pourrait donner lieu aux associations Cy Fy + Cy Fo; 


Ga Ay + As Gy interactions 1-2) et 4, F, + À, Fi 


Gy Cr + G: Cy (interactions 1-4) non observées. Il ne 
peut donc être responsable du f. f. observé. 
Par contre, aussi bien les tenseurs du champ cris- 


tallin (‘) que du couplage D-M (4) [1], [4] sont: 
susceptibles de donner lieu aux associations observées … 


de sorte que l’explication donnée par Bozorth [7] en 
termes du seul champ cristallin est correcte en ce qui 
concerne la symétrie, mais physiquement incomplète. 

L’auteur penche pour l'interprétation suivante. Le 
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“ 


f. f. observé est dû au couplage D-M aux températures : 


élevées. Le renversement du spin aux basses tempé- 
ratures peut être dû au champ cristallin. 


TABLE 1 


REPRÉSENTATIONS IRRÉDUCTIBLES ET VECTEURS DE BASE 


T; Az Gy Ca 
Te F3 Cy G 
T, F, #E 
F Ge Ay F, 


L’ambiguité — f. f, symétrique ou antisymétrique — : 


n'existe que dans les groupes comportant seulement des 


opérations de symétrie binaires, c’est-à-dire dans les 


groupes orthorhombiques et monocliniques. 

Dans les groupes ternaires (axes 3, 3, 6, 6) il n’y a 
aucune ambiguité : le f. f, ne peut être qu’antisymé- 
trique. Enfin dans les groupes quadratiques il est 
facile de séparer les cas possibles. Ces conclusions sont 
en plein accord avec une note récente de Turov [8]. 
Un mémoire plus détaillé est sous presse [9]. 


(*} On peut généralement négliger les causes b) et c). 
(*) Par l’intermédiaire d’un dyadique A, À, + A, AÀ.. 
. (9) Par lintermédiaire d’un vecteur A, A est propor- 
tionnel à la matrice moyenne du moment orbital pris entre 
l’état fondamental et tous les états excités. 
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ÉTUDE D'UN MASER A RUBIS 


Par M. SouTir, 


Laboratoires d’Électrostatique et Physique du Métal, 
Grenoble. 


Une équipe mixte comprenant un assistant de la 
Faculté, J. Galland et trois ingénieurs de la S. A. C. M. 
R. Buisson, R. Chicault et D. Descamps, a systéma- 
tiquement étudié dans le cadre du laboratoire les condi- 
tions du fonctionnement des masers à rubis travaillant 
dans la bande « $ » de fréquence signal, soit autour 
de 3 000 MHz. 

Le principe de fonctionnement de tels amplificateurs 
est bien connu et l’usage du rubis a été développé 
grâce aux travaux des « Bell Telephone Laboratories » 
et de l’Université de Stanford [1]. Le niveau fonda- 
mental de l’ion chrome dans une matrice d’alumine 
conduit, en présence d’un champ magnétique con- 
tinu À,, à quatre niveaux dont l’espacement dépend 
à la fois de 77, et de l’angle 0 entre Æ, et l’axe trigonal 
du champ cristallin. Grâce à une irradiation sufti- 
samment intense à la fréquence de résonance de la 
transition 2 — 1, il est possible (à très basse tempé- 
rature) de rendre égales les populations des deux 


. niveaux correspondants. Ce « pompage » peut conduire 


à une inversion des populations de 1 et 0 et par suite 
à une amplification à la fréquence correspondante, 
avec un facteur de bruit extrêmement faible qui confère 
un grand intérêt à ce dispositif. 

Le problème était d’amplifier un signal très faible à 
une fréquence de 3 000 MHz. La fréquence de pompage 
choisie a été 9 000 MHz pour des raisons de commodité 
d’appareillage. 

Une disposition convenable des niveaux (fig. 1) est 
obtenue lorsque l’axe Ô du cristal et le champ magné- 
tique appliqué Æ, font un angle petit, compris entre 0 
et 120 suivant la valeur exacte de la fréquence signal 
et pour des valeurs de 1, voisines de 1 000 Oe. 


1 
é 0 
13 
2 
3 
Fig: 
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L’échantillon est placé dans une cavité résonnante 
qui doit présenter le minimum de perte à son accord 
sur la fréquence signal et dont le couplage à cette 
même fréquence doit être ajustable pour permettre de 
régler le gain de l’amplificateur. En raison de son 
faible encombrement, une cavité semi coaxiale a été 
choisie, le couplage avec le coaxial d’alimentation 
étant obtenu par une boucle magnétique de surface 
variable (fig. 2). 


Coaxe 10 cm 


FAGN2A 


La surtension propre de la cavité est de 8 000 à 
l’ambiante et de 12 000 à la température de fonction- 
nement dans l’hélium liquide. Le résonateur est égale- 
ment. accordé sur la fréquence de pompage, ce qui 
permet d’assurer une saturation complète de la tran- 
sition 1-2 avec des puissances de pompage de l’ordre 
de 20 mW ; un couplage variable avec le guide d’onde 
est obtenu grâce à un iris et une antenne dont l’enfon- 
cement est commandé par un excentrique. 

Plusieurs auteurs [2] ont observé que le fonction- 
nement Maser ne pouvait pas être obtenu pour un 
assez grand nombre de configuration de niveaux théori- 
quement acceptables dans la bande de fréquence signal 
de 2 800 à 3 200 MHz. Ils ont attribué ce défaut à des 
phénomènes de cross-relaxation entre transitions dont 
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les fréquences seraient entre elles dans des rapports 
plus ou moins élevés. Pour l’étude systématique de ce 
phénomène, nous avons pu montrer que si le Maser ne 
fonctionnait que dans certaines configurations de 
niveau, c’est au contraire parce qu'il intervient alors 
des phénomènes de cross-relaxation beaucoup plus 
simples qui sont cette fois bénéfiques pour l’ampli- 
fication [3, 4]. Ainsi, par exemple, lorsque la tran- 
sition 3 — 2 de la figure 1 égale approximativement 
celle de 2 — 1 (pompage) et, bien que le pompage 
direct 3 — 2 soit interdit par les règles de sélection, la 
cross-relaxation provoque des retours de 1 vers 2 
accompagnés de montées de 3 vers 2 qui équivalent 
à un véritable pompage de 3 vers 2 avec égalisation 
des trois populations de 1,2 et 3. On peut alors montrer 
que l’excès de population disponible pour l’amplifi- 
cation est quatre fois plus élevé que dans le cas clas- 
sique à 3 niveaux [3]. 

L'utilisation de ce phénomène impose une condition 
supplémentaire qui peut être satisfaite à condition de 
jouer à la fois sur À; 0 et la fréquence de pompage et 
définit les conditions optima de fonctionnement. Cepen- 
dant, il existe une assez large plage de réglage autour 
de ces valeurs, grâce au fait que la condition de cross- 
relaxation Vos — V,, n’est qu'approximative, une 
certaine différence entre les deux énergies pouvant 
être absorbée par des réarrangements des spins entre 
eux avec échange d’énergie vers le réseau. 

Pour mieux définir les conditions de cross-relaxation, 
nous avons étudié les valeurs du temps de relaxation T; 
spin-réseau de la transition 2 — 1 (à 9 000 MHz) en 
fonction de la disposition des niveaux. Lorsque la posi- 
tion du niveau 3 est convenable, T, chute brusquement 
à cause de l’apparition de ce deuxième processus de 


LEres T É 
100ms 


ZE 


100ms aa 


0 mon x 


1 L 
ARE 


relaxation. Ce phénomène est naturellement très sen- 
sible à la concentration en ions Cr : il joue très peu et 
reste négligeable devant T; si les spins sont trop 
éloignés les uns des autres, maïs il se complique au 
contraire de cross-relaxations secondaires défavorables 
«lorsque la concentration est trop grande. Nous avons 
observé un optimum pour un dopage de 0,07 % en 
ions. La figure 3 montre le temps de retour à la normale 


d’un signal saturé à 4,2 0K en absence de cross-rela- 
xation (7, — 150 ms) ou en présence (7, — 20 ms). . 
Les performances propres du Maser sont dans ces 
conditions intéressantes. Elles se caractérisent par une 
valeur du produit du gain en tension par la bande 
passant de 12,5 MHz à 1,5 °K et une température 
équivalente de bruit de l’ordre de 65 °K. 


Section du Sud-Est, Lyon. 
Séance du 19 mai 1961. 
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POTENTIEL ET FONCTIONS 
D'ONDE DE MOLÉCULES DIATOMIQUES 
EN DEHORS DES CAS CLASSIQUES 


Par R. GRANDMONTAGNE, 


Institut de Physique Générale de l’Université de Lyon. M 


Les travaux du Pr J. Janin [1] ont permis de photo- « 
graphier les bandes du système B? EX — X?X5 de Nf «# 
jusqu’au nombre quantique de vibration v = 24. Le M 


potentiel de l’état inféricur suit bien la loi de Morse 


V = D[1 — exp (— ax)?] dont la fonction d’onde « 


est connue 


Fi -Ne 420): avec se ea 
Oe Le 
v ST 
DR PeRT Lo rE L(z) = Y ( Fe d le 
oirs v 1-0 DES) TA 


L'état supérieur s’en écarte nettement [2] dès que v 
est supérieur à 9. Des études ont donc été entreprises 


pour tenter d’obtenir des fonctions d’onde correctes en 


dehors des cas classiques (harmonique ou Morse). 

La première partie du travail consiste à rechercher 
la courbe de potentiel à partir des termes spectraux 
expérimentaux. L'application du principe de Franek- 
Condon n’ayant pas donné satisfaction ; on a utilisé 
la quantification de l’action préconisée par Rydberg [3] 
On dessine une courbe approchée (de Morse) et on 
évalue graphiquement l’aire : 


$ mv ds = Ê V2m(E — V) ar 


qu’on sait être égale à (0 + 1/2)4. La mesure graphique * 


n'étant pas en accord avec la théorie il faut corriger. 
Les tâtonnements nécessaires sont considérablement 
réduits par une méthode suggérée par R. Eido [4]. On 
assimile les ares de courbe y = V/E —— V à des ares de 
paraboles et on calcule aisément la correction de 
lPabscisse d’un point à partir de la différence entre 
l'aire théorique et l’aire mesurée. Une deuxième appli- 


cation est rarement nécessaire. La valeur absolue de la 
’ a 
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correction étant déterminée, on précise encore sur 
quelle branche on doit lappliquer en recalculant la 
constante de rotation B [5]. 

L'utilisation de la courbe de potentiel pour en tirer 
les fonctions d'onde a donné lieu à plusieurs essais. La 
distorsion de Gaydon et Pillow [6] s’est avéré inuti- 
lisable même après perfectionnements [9]. Les « r cen- 
troids » de Nicholls [7] ne s'appliquent qu’au cas de 
Morse et aux bas niveaux. Après un recours à la 
méthode B. K. W. qui est longue ; on a monté une 
- récolution de l’équation de Schrüdinger sur une ma- 
chine à calculer analogique [9]. Le principe en serait 
relativement simple : un « lecteur de courbe » donne 
une tension proportionnelle à £ — V ; un ; multi- 
plieur » fournit alors le produit de E£ — V par la fonc- 
tion Ÿ calculée elle-même par une double intégration 
de la machine à partir de Ÿ ; on annule enfin la somme 


9 / 
L' ne V)E = 0. 


Mais l'équation a une infinité de solutions non con- 
vergentes alors que les solutions cherchées sont les 
fonctions propres qui convergent vers zéro quand r 
tend vers zéro ou l'infini. La machine étant dans l’im- 
possibilité de faire elle même la distinction il faut lui 
fournir des « valeurs initiales » soit ŸW”’W#” pour un r 
donné. On a pu calculer à la main la pente relative du 
point d’inflexion gauche par 


: z — b L° 
PF af 17) 


où Z est le polynôme de Laguerre de la fonction 
d’onde et L’ sa dérivée par rapport à z. Le calcul a été 
possible jusqu’au degré 21, grâce à la translation pro- 
posée précédemment [9] z = y + b + v. Moyennant 
la fourniture d’un point et de sa tangente la machine 
permet l’enregistrement graphique du reste de la 
courbe. 

La méthode, mise au point pour le système négatif 
de l’azote est applicable à n'importe quelle transition 
de mojécule diatomique potentiel et fonctions d’onde… 


Section du Sud-Est, Grenoble. 
Séance du 30 juin 1961. 


BIBLIOGRAPHIE 


4] Janin (J.) et Eyraun (I.), J. Physique Rad., 1954, 15, 
88 $. 

2] Douccas (A. E.), Canadian J. Physics, 1952, 30, 302. 

3] RypserG (R.), Z. Physik, 1931,178, 381 et 1933, 80, 514. 

4] Eipo (R.), Thèse de 3° cycle, Lyon, 1959. 

5] GRANDMONTAGNE (R.) et Eino (R.), C. R. Acad. Sc., 

1959, 249, 366. 

[6] Przcow (M. E:), Proc. Phys. Soc., 1951, À 64, 772, 

7} Nicnozzs (R. W.) et JARMAIN (R. W.), Proc. Phys. 

Soc., 1956, À 69, 253. 

[8] Eipo (R.), Cahiers ce Physique, 1961, 15, 00. 

9] GRANDMONTAGNE (R.), J. Physique Rad., 1958, 19, 151. 
GRANDMONTAGNE (R.), Cahiers de Physique, 1961, 
15, sous presse. 


COMMUNICATIONS A LA SOCIÉTÉ FRANÇAISE DE PHYSIQUE 


843 


ÉTUDE EXPÉRIMENTALE 

DES FONCTIONS D’ONDE 
ET DES NIVEAUX D'ÉNERGIE DES ÉLECTRONS 
DANS DES MODÈLES SIMPLIFIÉS DE SOLIDES 


Par Guy MESNARD, 
Institut de Physique Générale de l’Université de Lyon. 


L'étude a été faite à l’aide d’un calculateur analo- 
gique du type à courant continu sur des modèles à une 
dimension.-On suppose le solide défini par le potentiel 
périodique V(x) dans lequel se déplacent les électrons ; 
ce potentiel, de période a, maille du cristal à une 
dimension, est pris identique pour tous les électrons 
considérés, quelle que soit leur énergie Æ. Les fonc- 
tions d’onde sont alors données par l’équation de 
Schrôdinger 


Len 
ee 2 


[W{x) — ET 9 = 0, 


le cristal étant supposé isolé et l’énergie d’un électron 
restant invariable, toute transition entre niveaux élec- 
troniques étant exclue. Cette équation se prête bien à 
une résolution analogique, le temps simulant l’abs- 
cisse æ. 

La fonction V(x) peut, dans certains cas, être obtenue 
à l’aide d’un dispositif électronique ou électromé- 
canique, mais on peut aussi la définir d’une façon très 
générale par une courbe; on trace une période de 
cette courbe sur le papier d’un cylindre enregistreur, 
de façon à en faire exactement le tour ; la tension élec- 
trique représentant V(x) est alors obtenue à l’aide d’un 
montage « suiveur de courbe » et l’on retrouve à 
chaque tour du cylindre les mêmes valeurs. On peut 
introduire des modifications à la périodicité, simulant 
des défauts de structure ou des impuretés, mais nous 
n’envisageons ici que le cas du cristal idéal indéfini.…. 

Les bandes d’énergie permises et interdites sont 
obtenues aisément et avec une bonne précision. Il 
suffit d'observer, pour les diverses valeurs de Æ, si la 
fonction d’onde reste finie ou prend finalement une 
amplitude de plus en plus grande. On peut toujours 
écrire 

d = AF, exp ik,x + BF, exp 1k,x 


avec k, + k, — 2nr/a, F, et F, étant des fonctions 
spaciales périodiques de période a et r un entier. 

Dans les bandes permises, on prend #, et k, réels et 
opposés et l’on a 


F;=f +18 


f et g étant des fonctions réelles de période a. La solu- 
tion la plus générale fournie par la machine correspond 
alors à des valeurs de À et B imaginaires conjuguées 
et peut se mettre sous la forme 


d = (af — Bg) cos kx — (ag + Pf) sin kx, 


Fa = f — ig, 


« et étant des constantes qui dépendent des condi- 
tions initiales choisies. En prenant la valeur initiale 
zéro pour 4, on a, pour toutes les valeurs de x mul- 
tiples de a, 


d = À sin x, 
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À étant une constante ; on en déduit aisément la valeur 
. correspondante de k. Ainsi on pourra trouver la rela- 
tion entre Æ et k dans toutes les bandes permises et en 
déduire, aux divers niveaux, les vitesses associées aux 
électrons et les masses effectives correspondantes. Dans 
les bandes interdites, les quantités X, et k, peuvent 
s’écrire 


k, = ik + Ina k, = — ik" + mrla 


k' étant réel, / et m étant des entiers de même parité. 
On peut en déduire que la solution réelle la plus géné- 
rale de l’équation de Schrüdinger est de la forme 


d = Afiexp (— Ka) + Bf, exp (Kx), 


où les fonctions réelles f, et f, sont telles que 


f(@ + a) = (— 1) f(x) ; 
elles ont donc une période égale à a ou à 24. Avec la 
condition initiale 4 — 0 on trouve que 4 = + Ash 4x 


pour toutes les valeurs de x multiples de a et il est 
facile d’en déduire les valeurs de X’ pour toutes les 
énergies des bandes interdites. 

On peut par ailleurs déterminer les fonctions f et g 


des bandes permises. Si on substitue Ÿ = F exp 1kx. 


dans l’équation de Schrôdinger, on trouve pour déter- 
miner F l’équation différentielle 


2m 


F'+9kIF REF = 


Nu) 


En prenant F = f + ig, on trouve que f et g sont 
donnés par le système différentiel 


d?f dg 5. a : 

D RE — [x A (F E)| j 
de df. [ème : 
LS + 2 = [x + (W —#)] : 


Ce système est facilement résolu avec la machine 
analogique. On choisit pour chaque valeur de ÆE la 
valeur de À déterminée préalablement. 

On obtiendra les fonctions /, et f, des bandes inter- 
dites d’une façon analogue. En partant de l’expression 
d = Fexp (— k'x) on trouve que la fonction F pst 
donnée par l’équation ! 


(FEI | 


RO M TE 
facile à résoudre à la machine. En choisissant, pour 
chaque valeur de Æ, les valeurs correspondantes, oppo- 
sées, de k”, on obtient, pour cette équation, deux solu- 
tions particulières périodiques, que l’on peut enregistrer ; 
elles fournissent les fonctions f, et f,. Si on ne connaît 
pas à l’avance la valeur de k’, il est possible de l'obtenir 
en examinant les solutions de l’équation précédente 
pour les diverses valeurs de #’. On obtient en général 
une solution dont l’amplitude tend vers zéro ou finit 
par croître indéfiniment ; entre ces deux cas, pour la 
valeur négative de k£’ que nous cherchons, on obtient, 
après un régime transitoire, une solution périodique. 

Ainsi sont aisément obtenues à l’aide d’un caleu- 
lateur analogique toutes les caractéristiques des fonc- 
tions d’onde des électrons, dans un cristal idéal indéfini 
à une dimension, 
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SUR LA PRÉSENCE DE LA RAIE H, 
DANS LE SPECTRE DU CIEL NOCTURNE 


Par Maurice DurAY, 
Institut de Physique Générale de l’Université de Lyon. 


Au cours des années 1959-1960-1961, une cinquan- 
taine de spectres du ciel nocturne couvrant la région 
5 500-7 000 À ont été obtenus à l'Observatoire de 
Haute-Provence en collaboration avec M. Nguyen H. 
Doan et sous la direction de M. J. Dufay. Le spectro- 
graphe utilisé comporte une chambre photographique 
à miroir ouvert à F#/1 (O0 — 14 cm) dont les aberrations 
de sphéricité sont éliminées par deux lentilles correc- 
trices. Le système dispersif est un réseau à échelons 
Bausch et Lomb (blaze à 1,6 u). Les spectres sont pris 
sur film et la dispersion est voisine de 50 A/mm 
(spectre du second ordre). 

La raie H,, dont la présence dans le spectre du ciel 
nocturne en l’absence de toute activité aurorale a déjà 
été signalée par de nombreux observateurs [1, 2, 3], 
est visible sur de nombreux clichés. On la distingue 
aisément au milieu des composantes de la branche P 
de la bande 6-1 de OH qui est entièrement résolue. La 


‘raie observée ne possède pas le profil élargi de la 


raie H, du spectre des aurores. C’est une raie fine dont 
la largeur est définie par le contour instrumental 
(1 à 2 À), et dont l’origine ne peut donc être, ccmm 
dans les aurores, un flux de protons animés de grandes 
vitesses. En outre, l'intensité de H, est très variable. 
Parfois totalement absente, la raie peut être nota- 
blement plus forte que la composante P, (3) de 
OH 6-1 et avoir une intensité comparable à celle 
de [OI] 6 363 À. Aucune corrélation n’a été trouvée 
entre son intensité et l'agitation magnétique ; il 
n’existe pas non plus de covariance entre H, et d’autres 
émissions du eiel nocturne (en particulier 6 300- 
6 364 À de [OI)). 

Par contre une relation manifeste avec le passage 
de la Voie Lactée dans le champ de l’instrument nous 
a conduit à reconsidérer l'hypothèse déjà proposée [4] 
d’une origine galactique. En effet, sur de longues poses 
d'hiver (20 à 30 heures) au cours desquelles la Voie 
Lactée n’est restée dans le champ de l’instrument que 
pendant une fraction minime du temps de pose (envi- 
ron © %), la raie H, est très faible. Au contraire sur 
des poses plus courtes (environ 8 heures) H, est très 
intense lorsque le spectrographe a vu défiler les régions 
H 1 de la galaxie pendant la quasi-totalité de la pose. 
Une vérification semi-quantitative de l’origine galac- 
tique de la raie H, a pu être effectuée sur une série 
continue de clichés (novembre 1960 à juin 1961) pris 
dans des conditions identiques (visée vers le Nord à 150 
au-dessus de lhorizon, pose 6 à 8 heures). Sur ces 
chichés, l'intensité de H, présente une variation assez 
régulière et à longue période. Admettant en première 
approximation que les régions H Z7 sont unifor- 
mément réparties le long du plan galactique entre les 
latitudes + 109 et — 109 et connaissant le champ du 
spectrographe, il est facile de tracer le graphique 
donnant, en fonction de la date, la durée de passage 
des régions H 77 dans le champ de l'appareil. Le gra- 
phique obtenu est bien parallèle à celui donnant aux 
mêmes dates l'intensité de H,, avec un maximum 
prononcé en avril 1961, Une vérification du même 
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ordre a été faite sur huit spectres pris par Kvifte en 


Norvège [2], l’auteur ayant publié en détail ses condi- 
tions d'observation, et il est apparu encore que les 
régions H Z7 étaient à l’origine de la raie H, quand 
elle était observée. 

Les résultats obtenus vont à l’encontre des idées 
souvent émises concernant l’origine de H,. On relie en 
effet généralement [5, 6] l’émission de H, dans le 
spectre du ciel nocturne à l'observation de H, dans les 
régions de la haute atmosphère non éclairées par le 
soleil [7] et on est ainsi conduit à admettre que la 
raie H, est émise par diffusion fluorescente dans l’espace 
interplanétaire ou dans une géocouronne entourant 
l'atmosphère terrestre. 

Il nous paraît au contraire manifeste que la raie H, 
observée comporte une composante galactique impor-- 
tante, parfois prépondérante. Des mesures précises de 
longueur d’onde sur le spectre le plus intense nous ont 
permis d'évaluer à — 0,4 + 0,1 le décalage en 
longueur d’onde de la raie observée et ce décalage est 
en rapport avec les vitesses radiales des régions H 77 
qu’on peut estimer comprises entre — 10 et — 50 km/s 
dans les régions visées. En outre, ainsi que Kvifte 
[2] l'avait déjà signalé, nous avons observé sur 
les spectres où H, est très intense la composante la 
plus forte 6583 À du doublet 3P —17) de [N 17]. 
Ce fait est également très favorable à l’origine galac- 
tique de H,. 


Section du Sud-Est, Grenoble. 
Séance du 30 juin 1961. 
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UN PONT POUR L’ÉTUDE 
EN BASSE FRÉQUENCE 
DE LA POLARISATION ÉLECTROLYTIQUE 


Par Je FreurY, 


Laboratoire de Physique $. P. C. N,., 
Faculté des Sciences de Bordeaux. 


L'appareil présenté a pour but de permettre étude, 
en courant alternatif de basse fréquence, du compor- 
tement d’une cuve électrolytique contenant, entre 
électrodes de divers métaux, des solutions aqueuses 
plus ou moins diluées. ve 2: 

On sait qu’il apparaît alors, superposée à la résis- 
tance classique de la cellule, qui dépend de la concen- 
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tration, une capacité de polarisation atteignant pour 
certaines conditions de concentration, de fréquence et 
de nature d’électrodes, des valeurs élevées, de l’ordre 
e dr dizaines de microfarads par em? d’élec- 
rode, 

Le montage utilisé est un pont de Nerst (type 
Sauty-parallèle) qu'on a perfectionné, en le compli- 
quant peu à peu, pour permettre la compensation des 
influences perturbatrices et l’utilisation d’une cellule de 


B 


ee 


ne ME 


mesure comportant une électrode avec anneau de 
garde. 


PRiNcipe. — Le circuit principal du pont est AEBF 
et comprend les impédances z, et z, formant la tête de 
pont (dont le rapport sera pris égal à 1) z, inconnue à 
mesurer et z;, impédance réglable pour l’équilibrage. 
Lorsque ce circuit principal est équilibré, on doit avoir 
entre les valeurs complexes la relation 2, z3 = 2% z 
d’où la valeur de z4. 

En réalité ces impédances ne sont pas seules pré- 


- sentes et, en plus des courants 1, to, t3, là qui les par- 


courent, il faut tenir compte d’admittances de fuite par 
rapport à la masse, indiquées sur la figure ci-dessus 
par 4;, a, 4, et a, et parcourues respectivement par £;, 
1, 1 et &. Contrairement à ce qui se passe habituel- 
lement, les admittances de fuite ne peuvent pas être 
faibles, à cause de la faible impédance de la solution 
et de la nécessité de disposer un anneau de garde pour 
uniformiser le champ à l’intérieur de la zone utile de la 
cellule. Les admittances a, et a; seront ainsi fortes par 
nature et a; et a, seront du même ordre que les précé- 
dentes, pour des raisons d'équilibre. 

L'équilibre du pont principal exigeait seulement 
l'égalité des potentiels instantanés Va — bp; mais 
l’électrode centrale doit être au même potentiel que 
l’anneau de garde, 4’où va — 0 — (0. Cette condi- 
tion va — 0 est réalisée en équilibrant le pont de 
Wagner &, 31 2 4, ; ceci se fait en agissant sur a;, le 
détecteur étant entre A et la terre, et les impédances 
z, et z, étant laissées intactes depuis l’équilibrage préli- 
minaire du pont principal. La condition d'équilibre du 
pont de Wagner s'écrit : 4, 3, = @ Z: ; Cet équilibre 
Va = Up — 0 supprime l'influence perturbatrice des 
admittances a; et a,, puisque ces impédances ont leurs 
deux extrémités au potentiel du sol. 

L'inconvénient de ce montage classique est que les 
deux réglages précédents du pont principal et du pont 
de Wagner réagissent l’un sur l’autre ; on ne peut 
procéder que par approximations successives, ce qui 
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est acceptable dans le cas d’admittances parasites 
faibles mais devient absolument impraticable lorsque, 
comme dans notre cas, les termes perturbateurs sont du 
même ordre que ceux du pont principal. 

Pour diminuer les tâtonnements, nous avons utilisé 
une remarque de Ogawa, reprise par R. Hérou (Revue 
Générale de l'Electricité, décembre 1944, p. 295). On 
calcule facilement que l’équilibre du pont principal, 
c’est -à-dire va — v8 est obtenu lorsque : 


VA 


——— (ai Ag — G9 Go). 
RE =NUX 


GTS nr Ua 

Pour annuler le premier membre, c’est-à-dire revenir 

à la condition d'équilibre sans perturbations, qu’on 
emploiera pour calculer z,, il suffit bien de réaliser 
va = 0, c’est-à-dire d’équilibrer le pont de Wagner, 
mais on aura davantage à annuler aussi la parenthèse ; 
on voit que cette nouvelle condition 4; 43 = 4, &9 
correspond à l’équilibre du pont Terre-AEB-Terre, 
qu’on obtiendra en agissant sur l’admittance variable 


‘En équilibrant successivement ces trois ponts, le 
second membre, produit de deux facteurs qu’on annule 
l’un après l’autre, restera petit, et la réaction mutuelle 
des ponts sera faible. En pratique, il suffira souvent de 
régler une deuxième fois le pont de Wagner et le pont 
principal pour avoir l’équilibre correct, qu’on vérifiera 
en passant successivement sur les trois ponts. Naturel- 
lement, chaque équilibre correspond au réglage de 
deux termes en quadrature l’un par rapport à l’autre. 


RÉALISATION. — Un contacteur à trois positions et 
5 circuits, noté 1-2-3 établit les connexions permettant 
le branchement de chacun des trois ponts. Un autre, 
xB permet de monter dans la branche inconnue soit (x) 
la cellule électrolytique, soit (B) une impédance fixe de 
comparaison. Les impédances variables sont cons- 
tituées, pour la branche étalon, par un condensateur 
variable de précision type G 701 Sullivan 0-700 PF, des 
condensateurs fixes à faibles pertes associés en décades 
sur contacteurs à galettes, boîtes de résistances du 
commerce. Les impédances destinées à l’équilibrage des 
ponts auxiliaires, dont la valeur doit être stable mais 
n’a pas besoin d’être connue exactement, sont réalisées 
par des condensateurs et résistances type radio asso- 
ciés en décades sur contacteurs classiques. 

On a ajouté en parallèle avec la branche inconnue z, 
une impédance réglable z, qui est destinée à compenser 
les capacités parasites et résistances de fuite intro- 
duites par le montage même dans la partie étalon 3, ; 
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en particulier, les commutateurs de décades, compre- 
nant chacun trois galettes à 12 positions présentent 
une assez forte capacité. Pour éliminer cette pertur- 
bation de la branche étalon, il suffit de placer dans la 
branche inconnue une impédance connue, par exemple 
une résistance pure (position $), donner précisément 
à z, cette valeur, et réaliser les équilibres en réglant 
le pont principal par z. On a ainsi ajusté dans la 


CPP Ampli. Oscillo. 


nee 


branche inconnue une admittance parasite qui contre- 


balance celle de la branche étalon. Quant à la légère 
variation de capacité des commutateurs dépendant 
de la position du contact, elle n’a pas d'importance 
puisque incorporée d’une manière fixe à des capacités 
qui sont étalonnées sur place. 


La source est un générateur basse fréquence sinu- 


soiïdal type GMWI1 L E A 20 Hz-200 KHz, en réalité 
l'utilisation de la sortie symétrique par transformateur 
incorporé limite à environ 50 KHz les fréquences utili- 
sables. Comme détecteur, on utilise un oscilloscope 
cathodique, relié encore par transformateur à écrans 
électrostatiques. 

Les résultats obtenus sur diverses cellules électro- 
lytiques seront publiés ultérieurement. 


Section du Sud-Ouest, Bordeaux. 
Le 2 juin 1961. 
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SPECTRE D'ÉNERGIE DES NEUTRONS 
ÉMIS LORS DE LA RÉACTION d(n,2n)p A 14 MeV 


Par H. JÉRÉMIE, M. HEymaAN et J. KAHANE, 


Laboratoire de Physique Atomique et Moléculaire 
du Collège de France. 


Le spectre d'énergie des neutrons émis lors de la 
réaction d(n, 2n)p à 200 dans le système du Laboratoire 
a été déterminé par temps de vol. La section efficace 
différentielle intégrée entre 2 et 8 MeV est de 


. 80 + 12 mb/sterad. Ce résultat est en accord avec la 


mesure du spectre d’énergie des protons émis dans la 
même réaction. 


A) Appareïllage. — Un faisceau de neutrons de 
14 MeV est défini à l’aide des particules alpha associées 
dans la réaction d(t, n)«. Ces particules alpha déter- 
minent aussi l’instant zéro dans la mesure du temps de 
vol des neutrons. 

Ces neutrons sont détectés par un cristal plastique 
o X 5 em monté sur un photomultiplicateur 56 AVP 
et les particules &« par une feuille de plastique mince 
montée sur un 5819 suivi d’un amplificateur distribué. 
Le retard entre les impulsions provenant de ces photo- 
multiplicateurs est mesuré par un convertisseur temps- 
amplitude du type à 6GBN6. Dans le faisceau est placé 
un récipient contenant le diffrseur de D,0. Pour pou- 
voir soustraire l’efflet de l'oxygène, chaque mesure est 
répétée avec H,0, la différence donnant les neutrons 
diffusés par le deutérium. Ceci entraine que la mesure 
n’a pu être faite qu’à un angle faible où l’énergie des 
neutrons diffusés élastiquement par l’hydrogène est 
encore voisine de l’énergie des neutrons diffusés élasti- 
quement par le deutérium. La base de temps de vol 
est de 60 cm ; le seuil de détection des protons de 
recul, de 1,2 MeV. 

La section efficace absolue est calculée en comparant 
les résultats obtenus avec ceux d’une diffusion (n, p) 


-sur le diffuseur H,0 à 459 dans les mêmes conditions, 


étant donné que la section efficace (n, p) est connue. 
Cette méthode évite la détermination de l'efficacité 
absolue du compteur. En transformant le spectre de 
temps de vol d\/dt en un spectre dNJ]d£E on a pas 
tenu compte du pouvoir de résolution de l’appa- 
reillage (largeur à mi-hauteur du pic y de la réaction 
(n, n’y)-est égal à 2.10% sec.). 


Résultats. — La figure 1 donne d0/(dQ dÆ) pour des 
neutrons. La figure 2 permet de comparer avec le 
spectre de protons mesuré dans la même réaction par 
K. Ilakovac et ali [1]. 

A. Migdal [2] a pour la première fois suggéré que 
dans les réactions avec le deutérium on pouvait 
s'attendre à l’émission de groupes de particules ayant 
des énergies peu différentes à cause de leur force 


attractive dans l’état final. C’est ainsi que Ilakovac 
trouve dans le spectre des protons un pic très prononcé 
à l’énergie maximum correspondant à deux neutrons 
de peu d'énergie chacun, et un pic plus large aux envi- 
rons de 5 MeV correspondant à un neutron et 1 proton 
d’énergies peu différentes, les deux particules ensemble 
emportant l'énergie maximum permise par cette réac- 
tion. Dans le spectre obtenu pour les neutrons on 
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F1G. 2. — Spectre des protons par K. Ilakovac et alii. 
observe le même groupe aux environs de 4,5 MeV, 
mais la section efficace est plus forte. Malheureusement 
on ne peut comparer directement les deux mesures 
parce que la mesure pour les protons a été faite à 100 
tandis que pour les neutrons on n’a pas pu descendre 
au-dessous de 200. La section efficace intégrée entre 2 
et 8 MeV pour les protons à 100 est de 37 mb/sterad 
et pour les neutrons à 200 de 80 + 12 mb/sterad. 
Comme Seagrave [3] a constaté une diminution lente 
de la section efficace des protons avec l’angle, il est 
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probable qu’elle diminue à 200 de telle sorte que la 
section efficace pour les neutrons est environs trois fois 
plus grande que la section efficace pour les protons. 
Ceci amène la conclusion suivante : 

Il y a la même probabilité d'émission pour un groupe 
proton-neutron que pour un groupe de deux neutrons 
d’énergie peu différente, de sorte que pour chaque 
proton, trois neutrons sont comptés. Le spectre d’éner- 
gie des neutrons a la même forme que celui des protons 
ce qui confirme l’hypothèse de l’égalité de la force 
proton-neutron dans l’état singulet. 

(Ce travail a été effectué grâce à une bourse du 
Commissariat à l'Énergie Atomique.) 


Lettre reçue le 26 juillet 1961. 
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SUR UNE EXPÉRIENCE 
DE CONVERSION ÉLECTRONIQUE 
UTILISANT UNE LAMPE A VAPEUR DE CÉSIUM 


Par MM. M. Perrot, J. MiLLET et G. JARRY, 
Institut de l'Énergie Solaire, Alger. 


Dans le cadre des recherches entreprises pour conver- 
tir directement l’énergie solaire en énergie électrique 
nous avons été menés à réaliser et à étudier en Labo- 
ratoire des diodes à vapeur de césium, le césium étant 
choisi à cause de son faible potentiel d’ionisation. 

Hatsopoulos [1] a montré que la conversion directe 
de la chaleur en électricité est possible à partir de 
l'effet thermo-électronique, si l’on met en regard deux 
métaux de travail d'extraction très différents. Le phé- 
nomène est limité par la charge d’espace au voisinage 
de l’émetteur et il faut un espace inter-électrode extré- 
mement faible pour obtenir un rendement satisfaisant. 
Wilson, Nottingham [2], ont montré à la suite des tra- 
vaux de Taylor et Langmuir, qu'il est possible de 
réduire cette charge d’espace et donc de prendre des 
espacements plus grands, en enfermant à l’intérieur de 
la diode une vapeur de césium ionisée au contact de la 
surface chaude émettrice. Nos expériences ont été 
menées en vue d'exploiter et de compléter ces diverses 
observations dans le cas de diodes présentant un espace 
interélectrode assez important. 


 tungstène) en la plaçant au foyer d’un miroir parabo- … 


Dans une première phase de Laboratoire, on utilise 
une lampe constituée d’un triple filament de tungstène, L'1 
entouré d’une cage de nickel, la distance du filament | 
à la cage est de l’ordre de 2 à 3 mm; ; l'ampoule qui 
contient le césium est placée dans un thermostat dont. 
on fait varier la température en maintenant constante | 
celle du filament qui est porté à 1 500 0C par effet 4 
Joule. de | 

Dans ces conditions, on observe que pour les plus 7. 
basses températures du thermostat le filament est 
positif par rapport au nickel. 5 | 

Lorsque la température du thermostat croit la 1 | 
puissance débitée passe par un maximum bien qu étant A 
faible et s’annule pour une température 7; (tempéra- 
ture d’inversion). : 

Au delà de 7;, la puissance augmente à nouveau 
mais le filament devient négatif. La valeur maximum 
par laquelle passe alors la puissance est beaucoup plus x 
importante que dans le cas précédent. 5 

$ 


Il existe, aux températures supérieures du thermos- 
tat, une seconde température d’inversion 73 qui 
redonne l'effet normal. E : 

Au cours de ces manipulations, nous avons observé « 
des rendements qui variaient de 8 à 12 %. : 

Ce phénomène complexe semble lié aux variations 
de températures de l’enceinte de l’ampoule, et par. 
conséquent à la pression de la vapeur de césium à. 
l’intérieur de celle-ci. Ce dernier paramètre agissant 
d’une part sur le nombre d’ions formés à la surface du …. 
tungstène et sur la nature même des surfaces, pouvant = 
en particulier, modifier complètement le travail d’ex-4 
traction de ces dernières. 4 

Dans une deuxième phase, nous avons tenté de réa- 
liser le chauffage de la surface émissive (une boule de 
lique et en l’entourant d’une grille de nickel. Une partie « 
de la lampe était calorifugée pour obtenir l'effet th a 
mostat. a 

Ces expériences en cours n’ont pas encore permis 4 
d'obtenir des rendements aussi importants que ceux 


ionization phénomena in gases, Uppsala, 1959-1 », » 


vérifiés en laboratoire. 1 
Lettre reçue le 9 Sctabrd st 1961. La 
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